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Preface 


The renewed interest in general relativity is due to 
several factors. For the theorist, quantization of Einstein’s 
gravitational field offers a challenging approach to a new 
theory of elementary particles and amelioration of the diffi¬ 
culties of quantum field theory. For the experimentalist the 
technological advances of the past two decades make possible 
some new gravitation experiments and more precise ways of 
doing the older ones. 

A book of this length on general relativity cannot be 
complete. I have attempted to give a fairly thorough intro¬ 
duction to the foundations of the theory, to the Riemannian 
geometry and tensor calculus ordinarily required in this field, 
to the conservation laws, and to the classical experiments. 

About a fourth of the book is mainly devoted to a treat¬ 
ment of the theoretical and experimental aspects of gravita¬ 
tional radiation. The last chapter gives a brief discussion of 
the deduction of the equations of motion, of unified field 
theories, of Friedman’s solution of the cosmological problem, 
and of the Hamiltonian formulation of general relativity. 

This tract is serving as a text for my relativity course in 
the Physics Department of the University of Maryland, 

I wish to express appreciation to Professor John A. 
Wheeler for the very stimulating year I was privileged to 
spend with him, to Professor James Anderson in whose 
lectures I saw these field equations for the first time, to 
Professor Peter G. Bergrnann for his many illuminating 
comments at the Stevens Institute of Technology colloquia, 
and to Dr. Charles Misner for his lucid remarks on a number 
of issues. Several sections of the book were substantially 
improved in consequence of discussions with a graduate 
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syeiit, i, George Hiiiiis. Messor C, J, Goebel read tk 
manuscript and gave valuable criticism, 

Tbe pbotograpb of Einstein whicb appears here was 
made available bv Mrs. T, 
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CHAPTER 1 


The EquiTalence Principle 


This assumption of exact physical equivalence makes it im¬ 
possible for us to speak of the absolute acceleration of the system 
of reference, just as the usual theory of relatimty forbids us to talk 
of the absolute velocity of a system. A. Einstein 

1.1 The Eotvos Experiment 

From the time of Newton it had been assumed that the 
ratio of the inertial mass to the weight (gravitational mass) 
of a body is the same for all substances. If we assume this and 
write the equations of motion for a body in the earth’s gravi¬ 
tational field, the mass cancels out and all freely falling 
bodies have the same acceleration. 

In 1890 Eotvos (1) performed an ingenious experiment 
designed to test the ratio of inertial mass to weight. Consider 
a mass on the earth’s surface (Fig. 1.1). There is a gravita- 



Fig. 1.1 

tional force G acting toward the earth’s center, and an 
inertial force I which is the centrifugal force associated with 
the earth’s rotation. The ratio of the two magnitudes, and 
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of the corresponding components of these forces, depends on 
the ratio of the gravitational to the inertial mass. Eotvos 
suspended two masses from a torsion balance, as shown in 
Fig. 1.2, at a latitude about midway between the equator 



and the pole. Suppose matters are arranged so that the 
balance is in equilibrium with the rod connecting the masses 
in the observer’s horizontal plane and pointing in the east- 
west direction. We can first conclude that the net torque 
component resulting from the vertical components of the 
resultant forces G -f- I on the two bodies is zero. If the 
ratio of inertial to gravitational mass is not the same for 
both, then the horizontal components of G -f- I will give rise 
to a torque which is canceled by an opposite torque of the 
suspension wire. If riow the entire apparatus is rotated 
through the angle tc , the bodies are interchanged and the sign 
of the torque associated with the horizontal components of 
G -j- I will reverse. The torque of the suspension wire, how¬ 
ever, remains the same. The result is that an angular deflec¬ 
tion of the rod and masses relative to the frame of the appara¬ 
tus will be observed if the ratio of inertial to gravitational 
mass is not the same for both bodies. 

Let the gravitational mass of one of the bodies be 
and let its inertial mass be Let be a unit vector from the 

body to the center of the earth, and let be a unit vector in 
the plane of the meridian, normal to the earth’s axis of 
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rotation. Let g, be the magnitude of the earth’s gravitational 
field. Then the gravitational force is given by 

= ( 1 . 1 ) 
Let (X be the earth s radius, let cd be its angular velocity, 
and let cp be the latitude. Then the inertial (centrifugal) 
force Ii is given by 

h = (rriiaco^ cos (p)i^ (1.2) 

Suppose the second body has gravitational mass and 
inertial mass 7712 - V^e compare the forces on it with those on 
the first body, by means of the torsion balance. Assume that 
and 2 3-re so chosen that the rod can be suspended in 
the center. Let the rod be represented by the vector b, and 
let the torque be denoted by T. We can write 

T = |x [G,-G,]+|x [1,-1,] (1.8) 


The resultant of the four forces must be in the direction 
of the thin wire which supports the rod, and is given by 

F = Gi+GgH-Ii+Ia (1*4) 

The component of the torque parallel to the supporting 
wire will tend to cause an observable rotation. Employing 

the preceding expressions enables us to write for the effective 
torque 


COS y [^i+m 2 ]i„,)-b X [Gi—G g+Ii—12] 

(1.5) 

In (1.5) we have omitted the centrifugal force in the denomi¬ 
nator since it is very small in comparison with the gravita¬ 
tional force. Evaluating (1.5) and making the substitutions, 

aj = Mj/mi; = Ma/mg (1.6) 

gives, for the effective torque, 

T , . CQS9P mim2(ai-aa)[b-Lxd 




(1.7) 
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Expression (1.7) vanishes if = ocg, and for ^ ag it 
will have a value which depends on the orientation of the rod 
b, with respect to the vector i, x which is normal to the 
meridian plane. It has its maximum when b points in the 
east-west direction. As we remarked earlier, the torsion 
balance is brought to equilibrium by turning it until the rod 
points in the east-west direction, in a plane tangent to the 
earth’s surface. Then the apparatus is rotated through the 
angle tt, reversing the sense of b; if ag there will be a 

torque which may then give a rotation of the rod relative to 
the frame which supports the balance. Eotvos observed no 
rotation and concluded that within one part in 10®, == ag 

for all the materials which were tested. This experiment has 
been repeated. (2, 3) The work of Southerns was done with 
pendulums and demonstrated the equality of a for radioac¬ 
tive materials. Professor R. H. Dicke (4) is now repeating the 
Eotvos experiment with greatly refined apparatus employing 
three bodies, and a threefold axis of symmetry, to minimize 
local disturbances. At this time his results agree with those 
of Eotvos, and the equality of a for certain substances is 
established to a few parts in 10^®. 

The Eotvos experiment enables certain conclusions to be 
drawn concerning the elementary particles. »The ratio of 
mass to weight for an electron plus a proton may be shown to 
be the same as for the neutron to one part in 10'^, and the 
reduction in mass of a nucleus resulting from nuclear binding 
forces can be shown to one part in 10® to be accompanied by a 
similar reduction in weight.^ With an accuracy of five parts in 
a thousand it can be concluded that the binding energy of the 
orbital electrons is accompanied also by a corresponding 
change in weight. 

Bondi (5) notes a possible distinction between mass which 
is acted upon and mass which is the source of a gravitational 
field. The mass which is acted upon he calls passive gravita¬ 
tional mass, and a mass which is a source is called active 
gravitational mass. The Eotvos experiment, in this view, de¬ 
termines the equality of the ratio of inertial and passive 
gravitational mass. 



THE EQUIVALENCE PRINCIPLE 


5 


1.2 Negative Mass 

Nothing in either Newtonian or relativistic gravitation 
theory precludes the existence of negative mass, but it is an 
empirical fact that it has never been observed. Both Newton¬ 
ian gravitation theory and general relativity indicate a quite 
different behavior for negative mass than for the correspond¬ 
ing situation in electrodynamics. If a small negative mass 
interacts with a large positive mass, again the (negative) 
mass cancels out on both sides of the equation of motion and 
the acceleration is still toward the positive mass. Thus a 
positive mass attracts all other masses, both positive and 
negative. A small negative mass would be expected to fall in 
the earth’s gravitational field. Similarly a negative mass re¬ 
pels all other masses, regardless of their sign. For a pair of 
bodies, one with positive mass and one with negative mass, 
with magnitudes about equal, we should expect the positive 
mass to attract the negative mass and the negative mass to 
repel the positive mass so that one chases the other! If the 
motion is confined to the line of centers the pair is expected 
to move with uniform acceleration. This problem has been 
discussed by Bondi (5). 

Schiff (6) has recently considered the possibility that the 
gravitational mass of an anti particle, the positron, might be 
nega-tive. His arguments are based on the renormalized 
quantum electrodynamics. The Coulomb field of an atomic 
nucleus produces a polarization of the vacuum. This effect, 
first calculated by Uehling (7) for hydrogen, produces a 
27 Me contribution to the Lamb shift of the 25 state in 
hydrogen. The virtual electron-positron pairs associated with 
vacuum polarization would be expected to contribute to the 
renormalized mass of atoms. We know from experiment 
that the inertial mass of the positron is positive. If the 
gravitational mass were negative, different atoms would be 
expected to have slightly different ratios of inertial to (pas¬ 
sive) gravitational mass. This follows because the relative 
contribution of virtual pairs to the mass would depend on the 
nuclear charge and its distribution. This would vary for 
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different atoms. For the case where the gravitational rest 
mass of a positron is assumed equal in magnitude and oppo¬ 
site in sign to that of a negative electron, but its kinetic 
energy is acted upon normally by a gravitational field, the 
difference between gravitational mass and renormalized 
inertial mass is finite and approximately equal (6) to 

{SmlSn)(Zj\S7Y r 

Here m is the electron mass, ju = mcjh, Z is the atomic num¬ 
ber, and F{(i) is the Fourier transform for the nuclear charge 
distribution, normalized to unit total charge. This expression 
has a ratio to the atomic mass of 2 x 10“’, and 

4.3 X 10“’ for aluminum, copper, and platinum, respective¬ 
ly. Since these numbers are larger than the uncertainties in 
the mass ratios determined by Eotvos, Schiff concludes that 
the possibility that the gravitational mass of the positron is 
negative is ruled out. 

It is likely that some experiment to see if anti neutrons 
fall in the earth’s gravitational field may be attempted. As 
we remarked earlier, if existing theories of gravitation are 
accepted they will be expected to fall, in any case. 

1.3 Equivalence of Different Frames of Reference 

The empirical fact that the two kinds of mass are equiv¬ 
alent did not fit anywhere in theoretical physics until 
Einstein (5) pointed out that it could be understood in terms 
of the equivalence of different frames of reference. Einstein 
postulated that an accelerated frame in a region free of gravi¬ 
tational fields is equivalent^ to a rest frame in a given infinite¬ 
simal region ^ of a time-independent gravitational field. 
This is called the principle of equivalence. From this postu¬ 
late the equivalence of gravitational and inertial mass can be 

t This implies that all observations made locally on a system in a uni¬ 
form static gravitational field will be the same as that on a system subjected 
to a uniform acceleration. 

§ The nonuniformity of a real gravitational field precludes replacing it 
by a single accelerated frame over a large region. 
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shown to be a consequence, by the following (8) argument 
(see Fig. 1.3). Here we have two systems Z and Z', whose 
2 axes coincide. The Z frame is assumed so be at rest in a 
gravitational field of strength g. The Z' frame is assumed to 


t 



cx 


(S 


Fig. 1.3 

have no gravitational field but it is being accelerated with 
acceleration g, in the positive z direction. Let Z' be at rest 
at i = 0 when a light pulse of energy is emitted from point 
a. This light is absorbed at point at which time the veloc¬ 
ity of Z' is gl/c. 

The light had energy and momentum EJc, The 
energy at point p can be obtained by use of the Lorentz trans- 
formation and is given by 

= Ea, COS ^p iEoc sin yj (1.8) 

where 

cos y) = (1 — sin yj = (— ivlc){l — 

Evaluating (1.8) gives 

Ef = £„[(c-f-w)/(c— u)]* E^(l+vlc) = Ecc+E^gllc^ (1.9) 

From the equivalence principle we assume that expres¬ 
sion (1.9) holds for the same process taking place in the Z 
frame. Imagine a mass M to be initially at a, in Z, then 
moved to Light of energy is emitted at a and absorbed 
by M at The total gravitational mass of M plus the ab¬ 
sorbed light is M'. Now lift M' back to a and re-emit light 
so that the mass at a is again M. There is no net change of 
energy during the process, so the change in energy going 
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from cx. to p can be set equal to the energy change on the 
return. 

Mgl -i- Efi — = M'gl (1.10) 

Making use of (1.9) then leads to 

M' - M = EJc^ (1.11) 

Expression (1.11) states that the increment in gravita¬ 
tional mass is the change in the inertial mass, and it therefore 
follows that the equivalence of mass and weight can be con¬ 
sidered as a consequence of the equivalence of an accelerated 
frame to a gravitational field. 

1.4 Gravitational Red Shift of Spectral Lines 

It follows also from the equivalence principle that we 
should expect a gravitational red shift of spectral lines. For 
consider again the emission of light at a by an atom in 
accelerated frame E', which is momentarily at rest. The light 
of frequency v is received by an observer at p who measures 
the frequency in units of his own proper time. The Doppler 
shift in the frequency at p gives rise to 

Vfi = Valic + v)l{c — v)]^ (1.12) 

and 

^ ’'a(l + (1*13) 

In the equivalent gravitational field, (1.13) again holds, and 

= vpjv^ = 1 + gljc^ (1.14) 

The quantity gi in (1.14) is the change in the gravita¬ 
tional potential, and we write the frequency shift as 

Av = v((p^~<p^)lc^ (1.15) 

In (1.15) 9 ?^ is the gravitational potential (a negative quan¬ 
tity) at the point where the light is received and 99 ^ is the 
gravitational potential at the point where the light is emitted. 
For light received on earth from a star, <p^ > 9 ?^. If M is the 
star’s mass, G is the gravitational constant, is the star’s 
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radius, and m and the mass and radius of the earth, 
respectively, (1.15) becomes 


Av 


vG 



(1.16) 


This predicts a red shift. In deriving (1.16) we have 
assumed that the frequency remains invariant in terms of 
the local proper time of the atom or molecule, even in a gravi¬ 
tational field. This assumption is really a statement about 
the atom. It would clearly not be valid for a pendulum clock 
but might be fulfilled at least approximately for a quartz 
oscillator clock,^ and more precisely for “atomic” clocks. 
In the absence of a complete quantum theory of atomic 
spectra which includes effects of all fields in interaction with 
a quantized gravitational field, it seems reasonable to sup¬ 
pose that (1.16) will be valid to a very good approximation. 
However, one might expect some effects due to the space 
derivatives of the gravitational field which would cause very 
small departures from (1.16) in some cases, even for an 
atomic or molecular system whose center of mass falls freelyJ 
Additional discussion of the red shift is given in Chapter 5. 


1.5 Further Remarks on the Equivalence Principle 

Since an accelerated frame is equivalent to a certain 
gravitational field it follows that we can annul a gravita¬ 
tional field by an appropriate acceleration. For example, a 
freely falling elevator in a gravitational field appears to be 
an inertial frame insofar as gravitational forces are con¬ 
cerned. A body would move within it as though no gravita¬ 
tional field were present, and no observation made on the 


t A quartz oscillator at rest on earth is compressed by its own weight, 
thus altering its dimensions slightly from the free-fall value it would have in 
a satellite. 

§ The fact that the atom is in free fall implies that there are minimal 
stresses resulting from the gravitational field. An atom at rest in an acceler¬ 
ated frame might be perturbed by such stresses. These might lead to some 
observable (and calculable) effects for experiments which are carried out 
in accelerated frames. See, however, C. W. Sherwin(18). 
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body could enable a distinction to be drawn between an 
inertial frame and the space inside the elevator. 

It is not so clear that this will still be true if the body 
within the elevator is electrically charged. There have been 
a number of investigations of the possible radiation from a 
uniformly accelerated point charge (9, 10, 11, 12). Bondi and 
Gold, and Fulton and Rohrlich predict that such a charge 
does indeed radiate. The radiation reaction (for uniform 
acceleration) is zero and the issue of energy conservation is 
complicated by the infinite self energy. It may be that 
when the internal structure of elementary particles is proper¬ 
ly taken into account, a charged particle will be found to 
radiate and have a non vanishing radiation reaction when 
falling in a uniform gravitational field. It would follow that by 
observing a charged and an uncharged body falling freely 
we can distinguish by local measurements whether we are in 
an inertial frame or falling freely in a gravitational field. 
The equivalence principle then becomes merely a guide for 
formulation of the equations of the gravitational field alone, 
and not a general law of nature. 


References 

1. R. V. Eotvos, Math. u. naturw, Ber. Ungarn 8, 65 (1890); Beibl. Ann. 
Physik 15, 688 (1891); Ann. Phys. 59, 354 (1896); Ann. Physik 68, 11 
(1922). 

2. L. Southerns, Proc. Roy. Soc. {London) 84, 325 (1910). 

3. P. Zeeman, Proc. Amsterdam Acad. 20, 542 (1917). 

4. R. H, Dicke, Revs. Modern Phys. 29, 3, 355 (1957). This does not 
discuss the new experiment hut contains Dicke’s ideas on the equiva¬ 
lence principle. 

A 5. H. Bondi, Revs. Modern Phys. 29, 3, 423 (1957). 

6. L. I. Schiff, Phys. Rev. Letters 1, 7, 254 (1958). 

7. E. A. Uehling, Phys. Rev. 48, 55 (1935); R. Serber, Phys. Rev. 48, 49 
(1935). 

8. A. Einstein, Ann. Physik 35 (1911). 

9. H. Bondi and T. Gold, Proc. Roy. Soc. {London) A229, 416-424 (1955). 

10. B. DeWitt and R. W. Brehme, Annals of Physics 9, 220, (1960). 

11. D. L. Drukey, Phys. Rev. 76, 543, (1949). 

12. T. Fulton and F. Rohrlich, Annals of Physics 9, 499 (1960). 

18. C.W. Sherwin, Phys. Rev. 120, 17(1960). 



CHAPTER 2 



The general laws of nature are to be expressed by equations 
which hold good for all systems of coordinates, that is, are co¬ 
variant with respect to any substitutions whatever {generally 
co-variant). A. Einstein 

2.1 The Idea of Covariance 

Minkowski made the important discovery that trans¬ 
forming from one inertial frame to another moving with 
relative velocity v corresponds (1) to the rotation of axes in 
a four-dimensional space-time coordinate system. The re¬ 
quirements of special relativity are indeed met in a most 
elegant fashion by writing the laws of physics as relations 
among four-dimensional vectors. While this procedure is not 
essential, it does add elements of beauty and simplicity. 

We employ coordinates labeled by superscripts. In the 
Minkowski space of special relativity we have the coordinates 
os, y, z, and ct, which we will denote x^, and x^. The 

important quantities are the relations between events. An 
event has no extension in space or time. It is a point in a 
four-dimensional space. The interval between two events a 
and b is denoted by the symbol and defined by 

If we transform from one Minkowski system of coor¬ 
dinates to another, the interval between the same events 
remains invariant and in terms of new coordinates x''^, x'^, 
x'^, x'^ is given by 

11 


( 2 . 2 ) 
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While the magnitudes of some quantities such as are 
invariant, others such as xj- ~~ do change under these 
coordinate transformations, which constitute a group. The 
form of the relation for s^b does not change, and such a 
relation is said to be covariant under the group of coordinate 
transformations (rotation of axes in Minkowski space). For 
two events separated by a differential interval we can write 

■ ds^ = - dx^^ - dx^^ - dx^^ (2.3) 


2.2 The Metric Tensor 

We shall employ a great many sums and shall follow the 
summation convention that an index which is repeated rep¬ 
resents a sum. This enables us to write (2.3) in the form 


—ds^ = gf^ydx/^dx^ 


(2.4) 


The elements of are represented by 

— 10 0 0 


g/iJ/ 


0 

0 

0 


+1 

0 

0 


0 

-j-l 

0 


0 

0 

+1 


(2.5) 


is called the metric tensor and the expression (2.5), which 
retains its form under Lorentz transformations, is called a 
Xorentz metric. 

In curvilinear coordinates, the metric tensor assumes a 
different form, but a coordinate transformation can trans¬ 
form it to the form (2.5), throughout the space, provided 
Euclidean geometry is valid there. 


2.3 The Metric Tensor in Curved Spaces and Accel¬ 
erated Frames 

The curved two-dimensional space on the surface of a 
sphere is described by the squared line element 


gf^ydx^^dx^ = r^dd^ •+- sin^ 6 dq)^; 


g/ti 


0 

0 sin^ 6 


( 2 . 6 ) 
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Triangles composed of arcs of great circles have the sum 
of their angles greater than n and less than Stc. This non- 
Euclidean, curved character is implied by the metric (2.6), 
through relations involving the derivatives of No coor¬ 
dinate transformation can reduce (2.6) everywhere to a 
diagonal form in which the diagonal elements are one. The 
concept of curvature will be discussed in more detail later. 

Suppose now that we have a plane triangle in an inertial 
frame. If we measure from a different inertial frame, the 
shape of the triangle may change, but the sum of its angles 
remains n and the geometry will not have changed. 

Consider a set of points which lie on a circle in an inertial 
frame. If we observe from a rotating frame of reference with 
center coinciding with the center of the circle and survey the 
same points, we find that the ratio of circumference to dia¬ 
meter now depends on our operational procedure for measur¬ 
ing length. To avoid, for the moment, the issue of clock 
synchronization within a rotating system, we assume that 
length is measured by bringing pairs of points into coincid¬ 
ence using the clocks of the inertial frame to define simul¬ 
taneity. The ratio of circumference to radius using measuring 
rods at rest within the rotating system will exceed 7 t, This 
follows because a rigid measuring rod would measure the 
same diameter as before, but when laid along the circum¬ 
ference it would be foreshortened by the Lorentz contraction. 
The metric tensor in the rotating frame now has to describe 
a non-Euclidean, curved space. 

In the inertial frame we may set up a system of fixed syn¬ 
chronized clocks throughout the space. In the rotating frame 
clocks at different radii will have different time measures. 
Identical clocks at different radii are similar to identical 
clocks in different inertial frames. They cannot be syn¬ 
chronized. Suppose that coordinate time is measured by 
pulses emitted by a clock at the center of the rotating frame. 
A clock at each point in the rotating frame ticks intervals for 
an observer at that point. The relation between these inter¬ 
vals and the coordinate time differences corresponding to 
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receipt of successive light pulses will be described by the goo 
component of This will now be a function of radius, and 
the metric (2.5) of special relativity does not apply.^ 

The equivalence of a gravitational field to an accelerated^ 
frame then implies that the special theory of relativity can¬ 
not be valid in an extended region (2) where gravitational 
fields are present. A curved-space metric is needed. To con¬ 
sider another example, suppose we make up a triangle in an 
inertial frame with sides which are light rays. The sum of the 
angles is n. If this is repeated in a gravitational field the 
sides of the triangle will become curved because of the action 
of gravitation on the energy of the photons. The sum of the 
angles of the triangle will now differ from n. In a covariant 
theory this is described by saying that the paths of light rays 

t An appropriate metric for the rotating frame is 

— ds^ = dr^-\-r^ d<p^~^dz^+2a)r^ dq)dt—(c^ — co^r^)dt^ 

^ A. few remarks on the clock paradox are in order. It was noted in 
Einstein’s first paper on the electrodynamics of moving bodies that if we 
have two identical clocks and keep one at rest in an inertial frame and then 
move the other in a closed path which returns to the position of the first 
clock, the two clocks will no longer agree. Consider now a pair of twins. 
One remains at rest in an inertial frame and the other sets off in a rocket and 
then returns. The traveler will on returning find that the stay-at-home twin is 
older than he. Darwin has pointed out (C. G. Darwin, Nature 180, 976 (1957)) 
that the entire problem may be understood within the framework of special 
relativity, for the acceleration times of the rocket may be short and the period 
of uniform motion extremely long. Then the result cannot depend on what 
happened during the short acceleration periods. If the velocity of the travel¬ 
er relative to the twin at rest in an inertial frame is v then 

where is the elapsed time for the twin at rest and is the elapsed time for 
the moving twin. This result follows immediately from special relativity, 
since we must carry out all calculations in the frame of the twin who remains 
fixed at all times in the inertial frame. The identical result may be obtained 
(see C. Moller, The Theory of Relativity, Oxford University Press, New York, 
1952) if we calculate in the frame of the moving-twin, using the formalism of 
general relativity, which is appropriate for frames that may undergo ac¬ 
celerations. There is, therefore, no paradox. 
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are always geodesics and that a curved-space metric is re¬ 
quired in a gravitational field. 

It has been pointed out (3, 4) that by suitably redefining, 
operationally, the measurement of length and time, the 
Lorentz metric may always be used. We could similarly 
insist that the earth’s surface is flat. By suitably defining 
the operations of measurement as a function of latitude and 
longitude, internal consistency would be achieved with 
Euclidean geometry. This is not the point of view which we 
adopt in this tract. 

2.4 General Covariance 

If we consider gravitational fields alone, the equivalence 
principle denies us the possibility of distinguishing, by local 
measurements, between an inertial frame and a freely falling 
system in a gravitational field. There is then no a priori 
reason to give special significance to inertial frames. Also it 
is not possible to set up the required system of synchronized 
clocks throughout a gravitational field. For these reasons 
Einstein was led to postulate that all systems of coordinates 
are equally good for the description of nature and that the 
laws of physics should have the same form in all. This is the 
principle of general covariance. 

If we adopt this principle the coordinates become noth¬ 
ing more than a bookkeeping system to label the events. 
The principle of general covariance has been a valuable guide 
in deducing correct equations. It leads us to avoid principles 
which seem simple only in certain coordinate systems, and to 
retain those which can be simply expressed in arbitrary 
systems of coordinates. It has been pointed out by Kretsch- 
mann (5) that any physical law can be written in a covariant 
form. The result is usually not simple. Kretschmann also 
pointed out that the principle of general covariance therefore 
has absolutely no necessary physical consequences. The re¬ 
quirements of simplicity of form plus covariance have none¬ 
theless been a valuable guide in deducing equations which in 
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a final analysis must stand or fall on the basis of comparison 
with experiment. 

The treatment of generally covariant equations in a 
curved space is facilitated by the formalism of the tensor 
calculus. 
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CHAPTER S 


Riemannian Geometry and Tensor Calculus 


The inner properties of surfaces are **most worthy of being 
diligently exploited by geometers.*' C. F. Gauss 

3.1 Some Ideas about Curvature 

A more general kind of geometry, in which g^j, is not 
necessarily reducible by a coordinate transformation every¬ 
where to the Lorentz metric, will now be discussed. 

Gauss considered the following question. Suppose we 
have a two-dimensional curved surface, inhabited by intelli¬ 
gent two-dimensional animals. Can they determine that 



their space is curved? Is it possible to determine the elements 
of curvature by means of measurements made within the 
surface alone? He found that this can indeed be done. First 
we proceed by labeling the points of the surface in any 
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regular but nonetheless arbitrary way. Two arbitrary 
families of curves, — constant and = constant are the 
coordinate system (Fig. 3.1). Direct measurement of length 
between points a and h gives ^22 5 similarly by measuring 
lengths ac and ad we obtain Gauss gave formulas 

which allow the curvature to be written in terms of the 
and their derivatives. 

Curvature is an intrinsic property, at any given point 
the same value is obtained in every coordinate system. We 
shall see how the idea of curvature can be extended to more 
than two dimensions. Einstein’s theory of gravitation relates 
the curvature of the space to the distribution of stress and 
energy. This follows in part the suggestion of Mach to the 
effect that the properties of the space-time continuum are 
determined by the distribution of energy. 

3.2 Transformation Laws for Different Kinds of 
Tensors 

Let us start our discussion with the assumption that we 
have M variables x^, x\ . . ., x^. A set of particular 
values of these variables is now regarded as a point in a hyper¬ 
space or manifold having M dimensions. The space is made 
up of all the points corresponding to the range of values which 
can be assumed by these variables. Suppose we employ a 
different labeling for the points x'\ x'^, x'\ . . ., x'^, such 
that 

^'cc _ . . ., X^) ( 3 . 1 ) 

We assume that derivatives exist, and write 

dx'<^ = —~ dx^ == —— dx^ (3.2) 

dx^ dx^ ^ ^ 

The coordinate differentials dx°^ are said to be the com¬ 
ponents of a contravariant vector.^ Similarly, any set of 
quantities F« are defined to be a contravariant vector if they 
obey the transformation law 

t A vector is a tensor of the first rank, a scalar is a tensor of rank zero. 
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dx'^ 

dx^ 




(3.3) 


Consider now quantities such as dcpjdx^, where 9 ? is some 
function of the variables x^, x^, . . x^: 

dcp _ d(p dx^ ^ ^ 

dx''^ dx^ dx'°^ (^•'^) 


The quantities d<pldx°^ are seen to obey a different trans¬ 
formation law from (3.3) and any set of quantities trans¬ 
forming according to 




a 


dx^ 

dx''=*‘ 



(3.5) 


are said to form a covariant ^ vector. 

Note that our definitions of covariant and contravariant 
vectors require the existence of derivatives on our manifold 
and not the existence of a metric. We are following the con¬ 
vention of denoting covariant vectors by subscripts and 
contravariant vectors by superscripts. 

The product of two contravariant vectors and 
will satisfy the transformation law 




dx'°^ 

dx^ 


dx'^ 

dx^ 


B^ 


(3.6) 


A set of quantities which obey the transformation 
law (3.6), are said to form a contravariant tensor of the 
second rank. Similarly a covariant second-rank tensor is one 
which obeys the transformation law 


T' 


flV 


dx^ dx^ rj. 
dx' dx'*' 




(3.7) 


A mixed tensor of any rank obeys the transformation 


t The word covariant has two quite different meanings. A covariant 
theory or equation has the same form in all systems of coordinates. The 
word covariant is used also to signify that a tensor obeys the transformation 
law (3.5). Thus a covariant equation may contain contravariant tensors as 
well as covariant ones, and other objects which are not tensors. 
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J is the Jacobian determinant \dx'<^ldx 'The superscript W 
is the power to which J is raised. T^\\ \ is said to be a tensor 
density of weight W. 

A function S, which transforms to S', such that S — S' 
at every point and in all coordinate systems is said to be an 
invariant, or scalar. 

In all cases when a quantity is given, its form provides a 
prescription for getting it in other coordinate systems. Thus 
if = A^B^, then in a different system of coordinates 
TV = 

The product A^B^ of a co variant and a contra variant 
vector transforms in the following way: 

= Sx ^ = A,B> (3.10) 


This product is therefore a scalar. It also follows that the 
inner product (summation over upper-lower index pairs) is a 
scalar for tensors of higher rank. These notions may be 
employed to test for tensor character. Let be an arbitrary 
contravariant tensor and let A ^ be a set of quantities which 
may or may not have tensor character. Then if the product 
Afj,Bf^is an invariant we can show that is a tensor, for 

A^B^== A'.B-^ (3.11) 


An index such as oc or /u, over which a sum is to be carried 
out, may be given any convenient letter. This is an impor¬ 
tant aid in manipulation. We may then write (3.11) as 
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It follows from (3.12) that A transforms according to 
A= A'y dec'll and is therefore a covariant vector. 

The Kronecker delta, written as d/, is a quantity which 
is unity H /u = v and zero it fz ^ v. If we write 


dcc^ _ dx'^ _ dx^ dx'^ 

dxf^ dx'^^^ dxf^ dx'^ dx/^ 


(3.13) 


it is evident that is a mixed tensor. 

If a tensor 5/**^“ = S it is said to be symmetric in the 
indices [j, and v. If A‘^^y — — A^°^y, the tensor is said to be 
antisymmetric (or skew-symmetric) in the indices a and 
By writing out the transformation laws it becomes clear that 
the symmetry properties of a tensor are retained under co¬ 
ordinate transformations if the pair of indices are both sub¬ 
scripts or both superscripts. In general the symmetry prop¬ 
erties are not retained if one index is a subscript and one a 
superscript. The symmetry properties are therefore meaning¬ 
ful only for the same kind of indices. 

We may write any tensor as 

Ap-^J = A- (3.14) 

and it follows from (3.14) that any tensor may be considered 
as the sum of a part which is symmetric and a part which is 
antisymmetric in a given pair of upper or lower indices. 

It follows from the transformation laws that if all com¬ 
ponents of a given tensor vanish in one coordinate system, 
then all components vanish in all coordinate systems. This 
fact is of great importance in theoretical physics. If a law is 
written in tensor form, for example by saying that one tensor 
equals another, the difference of the two tensors will vanish 
in all coordinate systems and the law has a validity independ¬ 
ent of the coordinates which may be employed. Similarly, 
if we establish a tensor equation in a special coordinate 
system, it is valid in general. 
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We have seen that the derivatives of a scalar form a co¬ 
variant vector. The derivatives of a covariant vector, on the 
other hand, do not have the required transformation law of 
tensors. For 

- A /A! A'] = ^—^ + A' ( 3 . 15 ) 

\ dxP “/ ^ 

A new kind of derivative, the covariant derivative, does 
permit us to form tensors from the differentiation of other 
tensors. In order to see how this comes about, we must con¬ 
sider the concept of parallel displacement of a vector. 

3.3 Parallel Displacement and Covariant Differentia¬ 
tion 

It is possible to extend the idea of parallel displacement 
of a vector to curved spaces in a consistent way. We assume 
that there always exists a reference frame such that the 
geometry is Euclidean (or Minkowskian) in the immediate 
vicinity of a point ^ P, and Cartesian coordinates are employ¬ 
ed in that neighborhood. In such a coordinate system the idea 
of an infinitesimal parallel displacement of a vector simply 
means that none of its components change. Also the scalar 
product of two vectors A and P, does not change 

under parallel displacement. For arbitrary coordinates we 
define the operation of infinitesimal parallel displacement of 
a vector from a point P to a neighboring point P' to be 
one which leaves the scalar product with an arbitrary vector 
Bp invariant. 

t We can always transform to a Lorentz metric at any given point where 
gjiP is well defined. Suppose the space lias n dimensions. The g^p at the point 
are a set of n{n +1 )/2 numbers. Write da;" = for convenience assume 

= m^ix. Substitute this in ds® and set g'^p = Sj^p. This gives n{n-\~l)/2 
equations for the same number of w’s which are then determined. There 
are other ways of doing this. It cannot in general be done over a region, 
because the m’s become functions rather than numbers, and the ecjuations 
da;“ == will not ordinarily be integrable to give the required coordinate 

transformation. 
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The idea of a parallel displacement along some given 
curve in a two-dimensional surface can be given an intuitive 
interpretation. Suppose the surface is developable. Then we 
can unroll it on to a plane and parallel-displace vectors in 
the plane. The surface is then rolled back and we have the 
required parallel-transported vector. If a given surface is 
not developable, we must first select a path for parallel 
transport, then erect a tangent plane at each point of the 
path. These tangent planes will envelope a developable sur¬ 
face. This new developable surface can then be unrolled and 
the operations of parallel transport and rerolling carried out. 
If the curve along which the parallel displacement is to be 
carried out happens to be a geodesic, it becomes a straight 
line when unrolled on to a plane. It is then clear that the 
angle between a geodesic and a vector remains unchanged in 
a parallel displacement. The angle between two vectors 
would also be expected to remain constant in a parallel dis¬ 
placement. We emphasize that the operation of parallel dis¬ 
placement from one point to another depends in general on 
the path which is specified between the points. 

In curvilinear coordinates in both curved and flat spaces, 
the components of a vector would be expected to change, 
under a parallel displacement. Let dA’' be this change, 
for an infinitesimal parallel displacement. dA’' should be a 
linear function of the coordinate differentials and the com¬ 
ponents A'', It can be written 

dA’' = —F^'ocfiA^dx^ (3.16) 

The quantities are coefficients whose properties are to 
be determined. First we show that symmetric in a 

and f}. Let yl*' be a coordinate differential dx’': 

d{dx’') == — F^^ffdx^ dx^ (3.17) 

We now return to the local Cartesian coordinate system 
by the transformations 
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x<^ = cc'^, . . .) 

x'«' = 0 ?^ . . .) 


(3.18) 


The primed coordinates are local Cartesian coordinates. 


dx^ = 


df^ 

dx'^ 


dx'fi 


(3.19) 


Under a parallel displacement, d{dx'^) == 0 by definition, 
so that from (3.19) we can write 


S{dx^) = 




dx'^ dx'y 


dx'^ dx'y 


dcp^ dcpy 


dx'^ dx'y dx“- dx^ 


dx°^ dx^ 


(3.20) 

Comparing this expression with (3.17) it follows that 



d^p d(p^ dcpy 
dx'^ dx'y dx°^ dx^ 


(3.21) 


The right side of (3.21) is clearly symmetric in the in¬ 
dices a and so must be symmetric in a and p also. 
Earlier we stated that the use of upper or lower indices is 
employed to give the transformation properties of tensors. 
Unfortunately the index notation is also employed for quan¬ 
tities such as which, as we shall see, do not obey the 

transformation laws of tensors. 

If the operation of parallel displacement from a given 
point to all other points in its neighborhood is defined for all 

vectors, the point is said to be affinely connected to its neigh¬ 
borhood. 

It is desirable now to employ the metric tensor 
which we defined earlier by 

~ds^ ~ g^^dxt^dx^ 

The tensor character of follows from the invariance of 
the squared “length” ds\ since all the coordinate transforma- 
ions which are employed here leave unchanged the “lengths” 

o curves. If the are arranged as a matrix and the inverse 

is taken, we obtain new quantities and 
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= a/ (3.22) 

Since the Kronecker delta is a mixed tensor, it follows that 
g''" is a contravariant tensor. The covariant and contra- 
variant metric tensors are useful for raising and lowering 
indices ^ to give new§ tensors; i.e., 

(3.23) 

We derive formulas for the in terms of the metric 
tensor. From the definition for a parallel displacement, 

d(A^A^) = ^ -J- dA^][A^ -j- SA/^] 

- (po“.) 

Carrying out these operations gives 

^ + g^.A/^SA- + = 0 (3.25) 


Making use of (3.16) to eliminate dA/‘ and dA" gives 

00 

= 0 ( 3 . 26 ) 

The symmetry of with respect to the lower indices 

allows permutation of the indices v and a to obtain 

= 0 (3.27) 

Similarly we write 

^ (^* 28 ) 

Solving (3.26), (3.27), and (3.28) gives, with the help of 
(3.22), 

t Note that in arbitrary coordinates, <5^^, is replaced by for 

= ^/ivf Only dji^ — 1 if /.t =b v and zero if ju:^v, 

in general. 

^ In general Bp/^ Btf-p. The left to right order of indices must be pre¬ 
served when raising and lowering them. 
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= is-'" 


-Sg 


VII 


+ 


igva Sg 


fiCC 


dx/^ 


dx^ 


(3.29) 


The three-index quantity 


r”' 

1 ‘ 
2 ■ 




q- 
dx'^ 


Sgya H 


flCL 


2 jxi^ J 


= r 


v,/ia 


■ +!iP' Phristoffel symbol of the first kind and is often 

is called the Christoffel symbol of the 

second kind and is often written as {^}. As we remarked 
earlier the are not components of a tensor. By starting 
with the differentiated transformation law for g^„ it is not too 
difficult to show that [/*«, r] transforms according o 

(3.30) 


[f.., vY = 1 ^. St, [^y. 


dx'^ dx'^ 


dx'^ dx'i^dx'^ 


Making use of the transformation law for then leads to the 
transformation law for b^ 

ry. * r\c\ 

(3.31) 


dx:> 8x1^ Sx’c „ , 8x'^ 8^x< 

~ flrV r)x/- 




A study of (3.31) indicates that it is always possible to 
choose a coordinate system such that all the Christoffel sym¬ 
bols vanish at a point. One way to do this is the following. 
Suppose the given point is labeled P and the Christoffel 
symbols do not vanish there. Carry out the coordinate trans¬ 
formation 

=^x^- ir“(P) + 4P“^5{p)(^^ -- 

Here the subscript (P) refers to the value at point P. 
Making use of (3.31) to directly calculate the new values of 
T'^^ shows that the new values vanish at point ^. This 
constitutes a proof that the Christoffel symbols may be trans¬ 
formed away because (3.32) can always be carried out wher¬ 
ever the Christoffel symbols are not infinite. A coordinate 
system in which the Christoffel symbols vanish at point P is 
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called a geodesic coordinate system: P is said to be the pole. 
It is also possible to transform away the Christoffel symbols 
along a curve (4). 

We return to our main objective, which is the formation 
of new tensors by differentiation. As we saw earlier, the 
partial derivative of a vector is not a tensor. A new kind of 
derivative is formed in the following way. Let a vector at a 
given point be Ai^. At a neighbouring point the vector is 
A/^ dAf^. The vector resulting from parallel transport to 
the neighboring point is A/* + SA/^. Subtract these quanti¬ 
ties to obtain 


dAf^ — SAf^ ~ 


dx- 


(3.33) 


(3.33) would be expected to be a vector since it is the differ¬ 
ence of two vectors at the same point; the quantity 


dAf^ 




then, is a mixed tensor called the covariant derivative oi A P' 
and written 

^ A 

^ (3.34) 

OX^ 

From d{A^Af^) = 0 it follows, using (3.16), that 

= r^f^pAocdx^ (3.16a) 

From this and a procedure similar to that of (3.33) and (3.34) 
we write the co variant derivative of as 

(3.35) 

The tensor character of (3.34) and (3.35) can be formally 
established by showing that they obey the required trans¬ 
formation laws, employing (3.31). The contravariant deriva¬ 
tive is formed by raising the index which denotes differen¬ 
tiation. 


(3.36) 
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The cova.ria.nt derivatives of a second-rank tensor 
may be obtained from the requirement that the 
remain invariant under a parallel displacement for arbitrary 
and BK Setting d{C^sA^B^) equal to zero gives dC^s, and 
the covariant derivatives may be written 

p = ^ - r-,, (8.37) 

Similarly, 

^ + r*„p cf- (3.38) 

OXP 

= + ( 3 . 39 ) 

The extension to higher-rank tensors is evident. 

Consider the covariant derivative of the metric tensor 
This is a tensor. In geodesic coordinates all components 
of this tensor must vanish. It therefore follows that 

^fiv',(r — ^ (3.40) 

in all other systems of coordinates. 

3.4 The Curvature Tensor 

Earlier we noted that if it is possible to carry out a 
coordinate transformation such that the metric tensor is a 
Lorentz metric everywhere, a region of space is said to be flat. 
A given vector may be parallel-displaced throughout such 
a space and give a constant vector field. Parallel displace¬ 
ment about a closed path results in a vector identical with 
the original one. In a curved space, parallel transport about 
a closed path does not, in general, yield the original vector. 
Consider, for example the surface of a sphere (Fig, 3.2) on 
which is a spherical triangle composed of geodesic curves. 
We start with the vector A and parallel transport it. It 
becomes B, then C, and finally D. Clearly A and D are not 
the same. We now proceed to calculate the change in com¬ 
ponents of a vector under a parallel displacement about an 
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Fig. 3.2 

infinitesimal closed path defined by four curves of the two- 
parameter set 


= (3.41) 

The path to be considered is shown cross-hatched in Fig. 3.3, 
in which the sides are at constant v and dt; at constant u[ 



Fig. 3.3 

The change in for the entire closed path is 

6nf^ — ~ ^ 


(3.42) 
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Considering the algebraic sum of the contributions from op¬ 
posite sides of the “parallelogram” and retaining terms only 
up to first order in Au and Av gives 

dni^ = {ri^a.^n^)Av ^ Au — ^ {r/^^yn«-)Au ^ Av 
dv ou ou cv 

(3.43) 

We carry out the operations indicated by (3.43), utiliz¬ 
ing the expression for the change in under parallel displace¬ 
ment. The result is then arranged in the form 

\ Pfi per __ pf, per 

dxy dx^ . . , , 

--— AuAv (3.44) 

cv ou 


Since only first-order terms have been retained it is 
more precise to write, from (3.44) 


lim 

Av—^0 


dnf^ 

AvAu 


8^_er^ 

dx^ 3xy 


dxy dx^ 

72“ --^ 

dv du 




(3.45) 


The left side of (3.45) is a vector, dx^jdu and dxyjdv 
are vectors because u and v are parameters. It follows, 
therefore, that the quantity defined by 




dxy 


+ r^ersF-^. 


-*• cry oLp 


(3.46) 


is a tensor. It is called the Riemann-Christoffel tensor, or 
simply the curvature tensor. For an infinitesimal closed path, 
then^, 

an/* = ^ ^ dudv ( 3 - 47 ) 

t This discussion implies that the result of two covariant differentiations 
depends, in general, on the order in which they are carried out. From (3.47) 
and (3.47a) it follows that 

;yS;a= — 
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Also for a parallel displacement d{n^n^) vanishes. This leads 
to 


dn 




fifiy ^a. 


Zx^ dxy 
du dv 


du dv 


(3.47a) 


From (3.47) and (3.47a) it follows that for any finite closed ^ 
path (requiring integration over u and i;), the vanishing of 
Rf^a/Sy guarantees that an arbitrary vector will not change 
under a parallel displacement. 

For a Lorentz metric in a given region the curvature 
tensor vanishes because of the constant g^y. If coordinate 
transformations are carried out in such a space, the trans¬ 
formed components of Rf^ocjsy will still vanish in consequence 
of its tensor character. It follows that a necessary condition 
for the space to be flat is the vanishing of On the 

other hand, if Rf^a./sy vanishes everywhere, then a Lorentz 
metric may be generated by parallel propagation of a 
(Minkowski) set of axes because this is now a unique opera¬ 
tion, independent of the path. The components of an arbi- 


t Since is a parallel-displaced vector, its value at any point within 
a given contour depends on the path taken. A meaning can be given to the 
integration of (3.47) or (3.47a) by carrying out the parametrization as in¬ 
dicated in Fig. 3.3. Suppose the initial point is O. Parallel-displace n“ from 
0 to 0^ along all the given curves. This generates a continuous vector field. 
The vanishing of (3.47) then guarantees no change in w“ as it is parallel- 
displaced around tlie (boundary) contour returning to 0. The entire matter 
can also be viewed using (3.10) alone. For a given contour, (3.10) may be 
written 


ds 






This is a set of c<mations to determine AM' as long as the tangent vector 
dxPjda is prescribed by specifying the path. If the change in AM is to vanish 
for all closed paths, the <piantity must be an exact differential, 

and this leads to 

dxy dxp ‘ 

Evaluating this with the condition that >1“ is being parallel-displaced then 
gives the requirement that the right side of (3.40) has to vanish. 
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trary vector will not change on parallel-displacement any¬ 
where. All the Christoffel symbols will vanish everywhere 
and all first derivatives of the metric will vanish everywhere. 
Since the metric was a Lorentz metric at the starting point it 
is then clear that it is a Lorentz metric everywhere. It fol¬ 
lows therefore that necessary and sufficient conditions for the 
space to be flat are the vanishing of all components of 
It is also clear that the vanishing of guarantees the 

integrability of the set of equations (3.16). 

Expression (3.46) tells us that Rf^a^y is antisymmetric in 
the indices ^ and y. The tensor R$afiy is given* by 






dx^ dx^ 8xy dx^ dxy 
+ gi,v{ri^apruy — 

From (3.48) it follows that 


Radfiy Rscc^y 


■R. 


R, 


^ocSyfi -^^fiyaS 

From (3.48) and (3.46) it may be verified that 

^^'fiyS + R^SBy + R^ySB = 0 


R 


ByS 

otByS 


+ ^xdfiy + ^ay6B ~ 


7 SB 


3.5 The Bianchi Identities 


^2 fj 

^ t>xy 

8a>^ dx^ 
(3.48) 


) 


(3.49) 


(3.50) 


An important differential identity may now be proved. 
Select at a point, a geodesic coordinate system; then the 
Christoffel symbols vanish and the covariant derivative of 
^^sBy becomes 




SBy'.v 


dx^dx" 8xr Sx” 


(3.51) 


t The rearrangement of (3.48) can be done by -writing ontgus^i/^aSvi 
then noting that ^ 


SM 


dxr 


— ig/iS )—rBo^p 


dxy 


^gjiS 

dxy 


|>y, 5] + [6y, fz} 


gftsRBoB ~ [a/5, and 
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From (3.51) it then follows that 

= 0 (3.52) 

The left side of (3.52) is a tensor which we have calculat¬ 
ed in a particular coordinate system and found all compo¬ 
nents to be zero. The result must therefore be valid in all 
other coordinate systems. The relations (3.52) are known as 
the Bianchi identities. 

A. mixed tensor of the second rank or higher can be 
summed over an upper and a lower index to give a tensor two 
ranks lower. This operation is called contraction. The tensor 
is formed by contraction of 

- sv. - ^ (,.m) 

and is called* the Ricci tensor. Its symmetry is evident. The 
scalar formed by contraction of R^^ according to 

gf^^R^y = R (3.54) 

is called the curvature scalar. § 

Multiplication of (3.52) by gives 

g^^lRy Sfiy-P 4 - Rysp0-,y + Rysyv,p\ = 0 (3.55) 

The covariant derivative of vanishes; also from (3.46) 
we have Equation (3.55) may then be 

written 

{Ry^ \dy^R).^^ = Gy^.y = 0 (3.56) 

(3.56) is the contracted Bianchi identity. The tensor 

i Contraction on. the first and second indices gives zero, since 
— and this vanishes because HiK^y is antisymmetric in i and k. 

Contraction on the first and fourth index gives the negative of (3.53). 

§ In a two-dimensional space the following relations are valid at any 
given point: 

^ ^ _ lim 

g ” 2 “ ~~ s-^~s 

Here g is the determinant of g^ip , rj and are the principal radii of curvature, 
S is the area of a small geodesic quadrilateral at the point, and Ad is the 
excess of the sum of its angles over four right angles. 
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Q V ^ R p—^dy^R is called the Einstein tensor. 


3.6 Geodesics 

Suppose we ask what is the equation of the curve defined 
by the requirement that each element of it is a parallel dis¬ 
placement of the preceding element. The equations of any 
curve are the one-parameter family of points 

= Ms) (3.57) 

The tangent vector is dx^jds, and the new tangent vector 
which results from a parallel displacement is 

r/> = ^ - rPaB ^ dxf (3.58) 


The new tangent vector T'p is also given by 


TV = ^ 4- A Z'—\ 
ds ds \ds } 


(3.59) 


Setting (3.58) equal to (3.59) gives the equations of the curve 
whose elements are parallel-displaced preceding ones, as 


d'^xP ^ dx^ dx^ 
ds^ ds ds 


(3.60) 


As might he expected, (3.60) represents the equations of 
the curves of extremal length, called geodesic curves. To 
show this we seek the relations which must be satisfied to 
give a stationary value to the integral 

jds — j V—gf^udxP'dx^ (3.61) 

We introduce a parameter k and write 


Let 


dxp‘ dx^ j. 


dxP dx" 


(3.62) 


dk dk 


(3.63) 
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The condition for a stationary value for (3.61) becomes^ 

d^Ldk = 0 (3.64) 

k may be chosen to be the length s along the geodesic. It is 
not the length for the other curves. 

Accordingly, (3.64) becomes 


The second term of the integrand may be written as the two 
terms 



and on integration the first of these contributes nothing 
because the variations vanish at the end points. The Euler- 
Lagrange equations for this problem are then 



We note that X = 1 along the geodesic. 
Making use of (3.63) gives 


d dL 


ds 



‘dg^^ dx^ 
ds ds 


_|_ 




d^x''' 


dL 


(3.66) 


dx/^ dc^ . 
dx^ ds ds 


and (3.66) becomes 

t Here <5 means a variation corresponding to change in coordinates and 
the vector dx^-Jdk as we go from one possible path to another; k merely labels 
corresponding points on different paths. 
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g 




KV 


ds^ 

gKV 


2 
d^ 




KV 


ds^ 


dx>^ 

1 

2 


+ 


dxf^ 

'^gKV 


dx^ dx^ 


ds 


dxf^ 


+ 


ds 

^gtcn 


dg 


fiv 


dx^ 


dx>^ 


dx^ dx^ 


ds ds 


0 (3.67) 


Multiplying (3.67) by and summing then leads to the 
geodesic equation 


d^x’^ Pa dxf^dx^ _ 

ds^ ds ds 


(3.68) 


(3.68) and (3.60) are identical. 

If we are dealing with the propagation of light, we have 
null geodesics. The deduction as given above cannot be ap¬ 
plied because ds vanishes throughout. Clearly the argument 
leading to (3.60) remains valid. It is therefore consistent to 
define the geodesic equation as the appropriate one for the 
null case, with the important proviso that the variable s is 
a parameter along the geodesic other than the length. 


3.7 Some Useful Galculational Aids 

We complete this section by giving a number of useful 
aids in manipulation of tensor quantities. First consider the 
CO variant divergence 

Ai‘.„ = ^ -h (3.69) 

From (3.29) we obtain 

= 4 g'*” ^ (3-70) 

(3.70) can be written in terms of the determinant of 
which we denote by g. The rule for expansion of a determi¬ 
nant leads to the relation 

= Zl/*" (8.71) 

Sg/.v 

In (3.71) is the cofactor of the element g^^. From 
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the rule for obtaining the inverse of a determinant, and from 
the definition of (3.71) may be written 

(3.72) 

dg = gg/‘’'dg^^ = —gg^^dg/‘'’ (3.72a) 

The expression on the far right of (3.72a) follows from 
d(.g^yg'‘’’) = 0. We then have 




ggf^v 




tsfiv 


dx°^ 




^g/iv 

dx^ 


(3,73) 


The use of (3.73) enables us to write (3.70) in the forms 


ll \g 


dgfiv 1 ^g 


(IV 


dx^ 


2g dx°^ dx°^ 


Expression (3.74) is now employed 
variant divergence (3.69) as 


(In V—g) (3-74) 
to write the co- 


For a contravariant tensor of the second rank, (3.38) and 
(3.74) give 

= :^g i (3.76) 

and for a mixed tensor, (3.39) leads to 

c/,, = (C»^V-g) - (3.77) 


For an antisymmetric tensor F“>^, the last term of (3.76) 
vanishes and the covariant divergence is 


Pafi 




(F^^W-g) 


(3.78) 


\/—g dx^ 

For a symmetric tensor 5“^, rearrangement of the last term 
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of (8.77) leads to 


oir-,0 


‘\/—g dx^ 


{SJV-g) 5^- (3.79) 


The Levi Civita symbol is often useful in perform¬ 
ing calculations. It is defined by the statement that e = 0 
if any two of the four indices are identical. It changes sign 
on interchange of any pair of indices, and == 1* A con¬ 
venient way of expressing the expansion of the determinant 


K of the tensor Kjj,y is 


Ke 


(xfiyS 


^ lkXja^ la ^lid 


(3.80) 


The transformation law 

dx'^ dx'‘^ 


K 


liv 


dxf^ dx^ 


K' 


afi 


may now be employed in (8.80) and the result rearranged to 
show that the determinants (K) and {K') are related by 


©)]’^ 


(3.81) 


Employing (8.80), (3.81), and the transformation law for 
then gives 

dx‘' dx'^ dx^ dxf^ 


<x,Py 8 



dx'‘^ dx'^ dx'y dx'^ 


(8.82) 


(3.82) tells us that is a tensor density of weight — 1 . 

Also it can be similarly shown that is a tensor density^ 
of weight + 1 . Note that this treatment requires the total 
number of indices on s to be equal to the number of dimen¬ 
sions in the space if s is to be a tensor density in the given 
space. We now let and from (3.81) conclude that 

the determinant of is related to that of by 

(-g')^ = /(J)(-g)i (3.83) 

t The metric tensor cannot be used to raise or lower indices ot e since 
Safiyd a-nd have different weights. 
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A four-dimensional volume element d^x' is related to 


d^x according to the Jacobi rule 


d^x 


(3.84) 

Combining (3.83) and 

(3.84) gives 



g d^x = V— g' d^x' 

(3.85) 


the quantity —gd^x is seen to be an invariant. 


3.8 Length Measurements 

We have been dealing with coordinates and intervals. 
In the special theory of relativity we deal also with lengths. 
A rod at rest in a given Lorentz frame has a well-defined 
length. The world lines of the ends of the rod are parallel to 
the time axis. In a different Lorentz frame with coordinates 
x'°^, the increment of length is given by 

{Ax'^)^ + {Ax'^f + {Ax'^)^ ( 3 . 86 ) 

Ax''^ refers to the difference of x'^ coordinates at the two 
points where the hyperplane x'^ — constant intersects the 
world lines of the ends of the rod (see Fig. 3.4). If we have 



Fig. 3.4 
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a pair of events we proceed differently and project each 
space-time point normally onto the hyperplane a?'® = con¬ 
stant, again using (8.86). 

In a curved space we can define length in the same way. 
Let the coordinates be a?®, a?^, a?^, and a?^, with metric 
Introduce at a given point a local Lorentz frame with coor¬ 
dinates cc'^, x'^, x'^. Let a?'o be parallel to x^, so that 

^ dx'^ ^ ^ 

—— == Oi - 0 

^a?® dx^ 

The index i refers to a space coordinate. For a pair of events 
we define 

dP — dx'^dx'^ == -— -— dxf^dx^ (8.87) 

dx^ dx^ ^ ’ 

The transformation law for the metric tensor is written 


^ _ Bx'^ Bx'^ Bx'^ dx'^ 

Bx^ Bx^ Bx^ Bx^ 

^ Bx'^ Bx'^ 

8.7i» ’ Sos- - 

dx'^ Bx'^ 
dx^ dx^ 

(3.88) 

Use of (3.88) in (3.87) then 

gives 



dl^ = U, - ^*=1 
\ ^00 / 

dx^dx^ = 

dx^dx^ 

(3.89) 

with defined by 





_ i0j^Oh\ 

^00 / 


(3.90) 


can be shown^ to be the reciprocal of the three-by-three 
matrix The determinant of the g^^ (four-by-four) is relat¬ 
ed to the (three-by-three) determinant 


t Write 

== < 5 ^^ g^^g.o-hg^^goo = 0 
then eliminate to obtain = <5^*. 
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of yjj,, which we denote by y, according^ to 

g = gooy (3.91) 

For a pair of events we may similarly define the time 
differential by its value in the locally Lorentz frame as 

= dl^ — g^^dx/^dx^ 

Employing (3.89), and extracting the square root gives 
dx'o = (—— {—gm)~hoid«!^ 

3.9 Determination of the Metric Tensor 

If we have a bendable (but not stretchable) rod and a 
“natural” clock which ticks at invariant intervals, these 
may be used to obtain the g^^ at a given point by direct 
measurement, using suitably chosen pairs of events, with 
known coordinates. No rigid rods exist, but ordinary ones 
may be used if the departure from complete rigidity is cal¬ 
culable. 

An alternative procedure makes use of a locally Lorentz 
frame and knowledge of the coordinates of events in both the 
original frame and the Lorentz frame. The squared lengths 
may be measured with rods in the Lorentz frame. No correc¬ 
tion for lack of rigidity is needed since the rods are in free fall 
and are not stressed. Measurements on six pairs of events 
then allows us to calculate using (3.89). 

Suppose the original time coordinate is established by 
the receipt of radio signals, and we have a caesium beam 
clock at the given point where the are desired. Compari¬ 
son with the time signals gives g^Q. To determine goi we 
arrange radio receivers at two points which differ in their 

t First set 


^00 

0 

0 

0 

^10 

Vxx 

yiz 

yia 

g^o 

yzi 

yzz 

yza 

gzo 

yai 

yza 

733 


Add (giolgoo) times the first column to the second column. Similar operations 
Involving the third and fourth columns then lead to (3.91). 
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coordinate. We then arrange to transmit a pulse of radiation 
from one receiver to the other, measuring the times of arrival 
in terms of the time coordinate signals. 

= 0 = (3-92) 

from (3.89) we have 

{dZi)2 = _ M\ [dx^Y (3.93) 

\ goo / 

(3.93) is then used to eliminate gii from (3.92), giving 
„ + + + ( 3 . 94 , 

Everything in (3.94) is known except goi, which is therefore 
determined. Similarly go 3 determined. Knowl¬ 

edge of these plus' the fact that the were measured then 
gives all the remaining components of 
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CHAPTER 4 


Field Equations of General Relativity 

and Electromagnetism 


Gravitation occupies an exceptional position with regard 
to other forces, particularly the electromagnetic forces, since the 
ten functions representing the gravitational field at the same time 
def ine the metrical properties of the space measured. A. Einstein 

4.1 The Gravitational Field Equations 

In Chapter 2 we discussed a rotating frame of reference. 
In consequence of the acceleration a Riemannian metric was 
needed. A nonuniform gravitational field is equivalent, 
within each small region, to an appropriate accelerated frame. 
A Riemannian metric would therefore be expected also to 
give a description of a gravitational field. The square of the 
line element is given by 

— ds'^ = (4.1) 

We identify the metric tensor as the gravitational 
field. This identification, which has followed from the equiv¬ 
alence principle, is perhaps the most important new feature 
of general relativity. The task now is to formulate the differ¬ 
ential equations which relate to the distribution of matter 
energy. The starting point is Newton’s law of gravitation, 
described by the Poisson equation for the gravitational po¬ 
tential cp as 

\J^(p ^ 4^7tGpM (^* 2 ) 

Here G is the constant of gravitation, pM f^e mass per unit 
■volume. 

The left side of (4.2) may be made at least Lorentz-in- 
variant by writing 

□ 9 ^ = 4i7tGpM (^■^) 
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(4.3) is similar to the set of equations for the four poten¬ 
tial in electrodynamics. There, the electric charge density 
is one component of a four vector in consequence of the in¬ 
variance of total electric charge. Mass, however, is not an 
invariant, so the right side of (4.3) is not a component of a 
four vector. It is a component of the second-rank stress- 
energy tensor already known from the special theory of rela¬ 
tivity. This second-rank tensor character leads to a more 
complex theory than electrodynamics. Also the particles of 
zero rest mass associated with the gravitational field, the 
gravitons, will have spin two, again in consequence of the 
second-rank tensor character of the gravitational ‘‘poten¬ 
tials, ^[iv- 

We recall the conservation laws. For a fluid in non- 
relativistic mechanics the conservation of mass is expressed 
by the equation of continuity 


_L_ c ^ = 0 


(4.4) 


is the velocity, and the summation over i includes the 
three space coordinates. (4.4) can be multiplied by dx^ dx^ 
and integrated over a given volume. The first three terms 
are readily transformed to an integral over the closed surface 
bounding the volume. The result then states that the time 
derivative of the total mass within the region of integration 
equals that carried in through the surface. 

The special relativity generalization of (4.4) is 

^ = 0 (4.5) 


T is the stress energy tensor. 

For the remainder of this tract we shall usually indicates 
the operation by a comma followed by a. Thus (4.5) 

becomes v = 0. 

In the theory of fields of special relativity a procedure i.s 
available for obtaining T from the Lagrangian density X, 
which is assumed to be a function of the field variables and 
their first derivatives The action function is the four- 
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volume integral of LIc, and the principle of stationary action 
states that for variations of the field variables whi< 2 ih vanish 
on the boundary, 

(5 J Ld^x = 0 ( 4 . 6 ) 

Here the symbol d is taken in the usual sense, and it is 
assumed that the field variables and are functions of a 
parameter. If the parameter is k, d means dk{dldk), with 
djBk meaning is constant. Thus the order of differentiation 
involving a coordinate and k may be interchanged. (4.6) can 
be expressed as 

r/dL . . dL 


I 


0 = <5 J Ld^x — 

d 


(I? - 


8q«. 


dq 




ar 

dxy dq 


d^x 4- 




d^x 


dL 


— 

/ 


(4.7) 


d^^x 


J dx^ \dq‘ 

The last term on the right of (4.7) is then written as a surface 
integral which is zero since are assumed to vanish on the 
boundary of the region. (4.7) has to hold for arbitrary varia¬ 
tions and the field equations 

d dL dL , , 

dxy'8^y~ ^ 

follow. We multiply (4.8) by and note that 

^ -J- * ^7% 

dx^ dq°^ dx^ dq°^^p 'dx^ 

also that = 9.^,p,y result is then rearranged to 

give 

\ r 

dfVL —q“,f 




0 


(4.9) 


From (4.9) and (4.5) it follows that 

dL 

T^y = ^ (4.10) 

,r 

From classical mechanics we recognize —as the energy 
density, since the correspond to velocities. A study of 
(4.9) and (4.10) indicates that the definition (4.10) does not 
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make T/ unique. Any quantity can be added to T 

provided « is antisymmetric in y and a. It has been 
customary to select the function 0 in such a way that the 
resulting tensor T (both indices raised or lowered) is sym¬ 
metric (1,2). The reason for this is the desirability of having- 
the angular momentum density, defined by 

Mafiy = T^pXy - T^yXp (4.11 ) 

satisfy a conservation law 

= 0 ( 4 . 12 ) 

Carrying out the differentiation indicated by the left side of 
(4,12) shows that (4.12) is indeed satisfied by a symmetric 

^a/?* 

For a fluid, is given in terms of the four velocity 

'I'IIV = (p + /lU V + ^iivP (4.13) 

Here p is the pressure and E is the total (mass) energy 
density, evaluated at each point in the rest frame of tlie 
local matter. Equation (4.13) is the four-dimensional gener¬ 
alization of the three-dimensional stress^ tensor which 

gives the force dF^ transmitted across the element of surface 
ds^ as 

dF^ = Tij.ds^' (4.14) 

(4.4) gave the conservation of mass; (4.5) includes the 
conservation of both total energy and momentum. The mo¬ 
mentum and energy are given by the space and time com¬ 
ponents of with 

■Pi = J Tf>dx^dx^dafi; P’0 = J dx^ dx^dafi 

Pq, transforms like a four vector, under Lorentz transfor¬ 
mations. 

t Historically this is where the stress energy tensor had its origins, 
in elasticity is symmetric in consequence of the balance of turning mo¬ 
ments on volume elements. 



GENERAL RELATIVITY AND ELECTROMAGNETISM 


47 


In electrodynamics^ the stress energy tensor is given in 
terms of the field tensor by 

r/= -^ (4.15) 


This follows from (4.10) by using the Lagrangian density 
for the Maxwell field which is given later in this chapter. The 
term 


1 

4b7t 






was added to make symmetric. 

These expressions were obtained from arguments ap¬ 
propriate for Lorentz frames. However (4.13) and (4.15) are 
tensor equations and are therefore valid in arbitrary coor¬ 
dinates. Any result obtained from (4.10) will either be in 
tensor form or can be readily modified (for example by chang¬ 
ing ordinary derivatives to covariant derivatives) so that it 
transforms as a tensor. A manifestly covariairt expression 
for T is given later in this chapter. 

We return now to (4.3) and note that the right side of 
it is identified as the zero-zero component of This sug¬ 

gests that the gravitational field equations be formulated 
as a second-rank tensor relation set equal to the stress 
energy tensor associated with the other fields. The left side 
of the tensor field equations should reduce to the D’Alember- 
tian operator in a certain level of approximation. The number 
of possibilities for this tensor is reduced by the requirements 
that it be formed from and contain no higher than second 
derivatives of It should reduce to the D’Alembertian for 
weak fields. 1\) guarantee this and to ol)tain dimensional 
consistency in a simple way it is reasonable to require that 
the equations at least be linear in tlie second derivatives of 


the 

Suppose then that a gi ven tensor does not contain higher 
derivatives of g^^;, than the second and it is linear in these. 


t VVe arc usinjr C.G.S. absolute units. 
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Introduce geodesic coordinates so that the Christoffel sym¬ 
bols vanish. From (3.4)8) we see that the second derivatives 
of gfiv can in these coordinates always be written as linear 
functions of the curvature tensor. If a second-rank tensor is 
needed, its most general form must therefore be 

“h ^tiv (4.16) 

Cl, Ca, and C 3 are constants. The first two terms of 
(4.16) are linear combinations of curvature tensor com¬ 
ponents. Since (4.16) is a tensor equation, it is valid in any 
system of coordinates. The requirements of general covari¬ 
ance, a second-rank tensor, linearity in the second derivatives 
^f plus absence of higher derivatives, have given us the 
form (4.16). 

An argument of Mach (3) and Hilbert ( 4 ) is useful at 
this point. Suppose the laws are in the form of a relation 

^fiv fiv (4.17) 

between symmetric tensors. Assume that a solution is 
obtained, which gives the as a function of the coordinates. 
Coordinate transformations can be carried out by introducing 
the four functions 

= Fp'{x) (4.18) 

These functions may be selected in such a way that = 
and the first derivatives of g^^ are everywhere equal to 
those of g'^p at some initial time or on a spacelike surface.^ 
Because of the covariance of (4.17), the transformed field 
equations will have exactly the same form as the original ones. 
Therefore since the equations do not have higher derivatives 
^f ^fiv than the second, it follows that g^^ = g'iiv everywhere. 
This would be in contradiction to 




flV 


dx^ dx°^ 


(4.19) 


It was therefore concluded by Hilbert that a four di- 


J Consider for example — x^, x'^ ~ ~ x^ x'^ = 

then g'l^v = gfiv everywhere at = 0, and all first space and time deriva¬ 
tives of are the same at £»® = 0. 
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mensional covariant expression such as (4.17) should not be 
ten independent equations, but that four identities would 
have to hold such that the set consists only of six independent 
equations. In this case the solutions would contain four 
arbitrary functions which become uniquely specified when 
the coordinate system is singled out in some noncovariant 
way. A clue to the four identities which have to be satisfied 
by (4.17) comes from the fact that with a Lorentz metric the 
stress energy tensor which forms the right side of (4.8) satis¬ 
fies the conservation law 


T 


V 

H ,v 


= 0 


(4.20) 


The logical extension of (4.20) for a generally covariant 
theory is 

T/,, = 0 (4.21) 

It was shown (Eq. (3.56)) that the Bianchi identities 

give 

(R/ — = 0 (4.22) 


These identities apply if the constants in (4.16) are suitably 
chosen. From these arguments it follows that the field equa¬ 
tions must be given by 


R/ — — Ad/ =:KT/ (4.23) 

Here A and K are constants. Experience (5) and logical 
simplicity have indicated that A may be set equal to zero. 
Einstein’s original (7) formulation included the statement 
that the motion of particles in a gravitational field is describ¬ 
ed by the co variant generalization of Newton’s laws, which 
gives the geodesic equation 




+ 


dx^ dxy 
ds ds 


= 0 


(4.24) 


Later it was shown that the equations of motion are in 
fact already contained within the field equations (4.23) and 
do not need to be postulated separately. According to (4.24) 
the force per unit mass on a body at rest is given by the three 
components —c^jT^oo* weak field approximation the 
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are very close to the Lorentz metric values, and for a 
time-independent metric 

By analogy with other fields, it is then natural to regard 
i^^gfiv as the gravitational potentials. For weak fields, 
(4.23) should reduce to the Poisson equation,^ (4.2). This 
can be accomplished by choosing the constant K in (4.23) 
equal to SjtGjc^ where G is the gravitational constant. 
G — 6.67 X 10 ~® cm^ gm~^ sec"^. The field equations then 
become 

IT-. StcG ^ 

Ri.v — \gf.vR = —2"^. (4.25) 

4.2 Variational Principle Deduction of the Field 
Equations 

We start with the action function /, given by 

I ~ 1q -F I p :==. j {^Lq -f- Lp)^/—g d^cc (4.26) 

Here Iq is the gravitational part and Ip is the part of the 
action function due to all other fields, Lq is the Lagrangian 
density for the gravitational field, and Lpis the Lagrangian 
density for all other fields. For JLg we select the curvature 
scalar R times a dimensional factor The principle of 

stationary action then gives 

{c^ll67cG)d j RV—g -f- dip = 0 (4.27) 

For the variation of the gravitational part of the action func¬ 
tion we have then 

8Ig = (c®/163tG) [JJ V'^'g-“'')d*a!] 

(4.28) 

We recall that 

Rfiv ~ ~ r<^p,a,v -h 

First let us choose a geodesic coordinate system; then we 
may (6) write 

J The weak-field treatment is carried through in Chapter 7. 
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^Rfiv /ia)',v (4.29) 

This is a tensor equation which is therefore valid in all coor¬ 
dinate systems. 

The variation of — g is needed. (3.72) gives 

d{~~g) = = ggfiv^gf^'' (4.30) 

Expressions (4.29) and (4.30) and the vanishing of the 
covariant derivative of the metric tensor enable us to express 
(4.28) in the form 

dIo = (<^ll6nG) r f [§/*>’ 

L-*' _ (4 34 \ 

+ V-g 

The first integral in (4.31) will now be shown to vanish 
(6). The quantity gP'^dR^y is a scalar; from (4.29) it follows 
then that the quantity is a vector. 

Expression (3.75) then gives 

■’ . ,_ (4.32) 

= /(V-g[g'“'<5/’V - 

Gauss’ theorem may he used to convert the right side 
of (4.32) to a surface integral which is zero in consequence 
of the vanishing of the variations on the boundary. (4.31) 
becomes 

sro = (y^) J K]5g^'' (4.33) 

This result was obtained witliout making entirely def¬ 
inite statements concerning what (in addition to g^^) are 
the variables. For the remainder of the action function it 
will be assumed that no higher than first derivatives of the 
are present. For the variation of I we have 

8{LFV—g) , SiLp V—g ) 



d^x (4.34) 
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The second term of the integrand in (4.34) may Tbo 
written 


d(L^V-g) 








{Lp V-g) 


,a 


LSg'"'. 


{LpV-g) 


a 


{4S.SB ) 


dgf^'’ 


,a 


The first term on the right of (4.35) can be transformed 
into a surface integral which again is zero because the varia¬ 
tions vanish on the boundary. This gives 


dl, 


=i/[ 


^gilV 


{Lp a /~^) 


( 


sg'“',« 


{Lp y'- 


g)) 1 dgi“’d*cc 

(4.3e ) 


We set the bracket of the integrand of (4.36) equal to 3. 
second-rank tensor density (T^j, will in a momeixt; 

be identified as the stress energy tensor.) 



d 

dg/^^ 


(Lp '\/—g) 




(4.26), (4.33), (4.36), and (4.37) then lead to 


(4.3r ) 



SnG 



(4.33 ) 


The left side of (4.38) satisfies the Bianchi identities;» 
so defined by (4.37) satisfies 

T /,^==0 ( 4 . B £>) 

and may therefore be identified as the stress energy tens oar ^ 
subject to the usual issues of uniqueness. The prescription 
(4.37) for calculating T is in some respects superior to 
(4.10) since a symmetric quantity always results. 
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4.3 Maxwell’s Equations 

The Lagrangian density for electrodynamics is 


L 




M 


a/? 






(4.40) 


Lp is the Lagrangian density of the charged particles. 
/« is the four-current density. Consider first a Lorentz metric. 
The field tensors are obtained from a four potential A 
with 

F fip = A A (4.41) 

The components of the four potential are considered to be 
the field variables in (4.40). The Maxwell equations are given 
in terms of the field tensors and the four-current density as 



(4.42) 


dcc^ 


(4.43) 


Substitution of (4.41) into (4.43) shows that (4.43) 
becomes an identity satisfied by any F obtained according 
to (4.41). For a field tensor obtained from a four potential, 
(4.42) is therefore the entire content of Maxwell’s equations. 
If we adopt the Lorentz gauge, given by 

A^^y = 0 (4.44) 


then substitution of (4.41) into (4.42) gives 






(4.45) 


as the four-potential formulation of electrodynamics. 

We wish now to generalize these results for arbitrary 
systems of coordinates. (4.40) is still valid and (4.42) becomes 

Ff^^.y =■ (47r/c)f-“ (4.46) 



54 


GENERAL RELATIVITY AND GRAVITATIONAL WAVES 


(4.41) becomes which by (3.35) and the 

symmetry of lower indices again reduces to (4.41). 

Since is an antisymmetric tensor we may use (3.78) to 
write 

-- (Ft^W—^),v == {4>7clc)j/^ (4.47) 

v—§ 

The appropriate generalization of (4.43) is 

0 (4.48) 

However, the quantity in the brackets is a second-rank tensor 
antisymmetric in the indices y and 6, so it follows, again 
from (3.78), that (4.48) reduces to (4.43) in arbitrary systems 
of coordinates. 

The expression (4.44) becomes 

A^.^ = 0 ( 4 . 49 ) 

A calculation, which is most easily performed in geodesic 
coordinates (noting that the Christoffel symbols, but not 
their derivatives, vanish), leads to the following result for 
changing the order of covariant differentiation of a vector: 

= T?" A«. 

ip] or cr; p acrp^ 

We apply this now to (4.46), which is first written in terms of 

Af^. 

— A^’^.^, = (47i/c)y/^ 

Lowering the index [z, writing in terms of 

and using (4.49) gives 

Ap-v'"' — FpaA°^ — {~4i7ilc)j^ (4.50) 

as the generalization of (4.45). 

4.4 Motion of a Charged Particle 

The action function for a charged particle with charge e 
is given by 

I = —me j ds + {eje) j A,, dx/‘ 


( 4 . 51 ) 
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The equations of motion result from setting the varia¬ 
tion of I equal to zero. For the first term of (4.51) we obtain 
from (3.62) through (3.68) 


■mcd ds 


J' 


me 


J( 


„ dx^^ dx^ 


ds^ 


ds ds 


^g^ydxyds (4.52) 


For the second term of (4.51), 

(elc)d j A^dxf^ = (eje) j dA^dxi^ -j- {eje) ^ A ^ddx/^ 

= (e/c) j ld{A^dxf^) — Sx^dA^dA^dx^] 


(4.53) 


The first term in the integrand on the far right vanishes when 
integrated because the variations are zero at the end points 
of the path. For the remaining terms we note that 


while 


dA. = °Aj: dxP 
dxP 

d 4 

SA = ^ Sx^' 

^ dx^ 


Making use of (4.41), combining (4.52) and (4.53), setting 
the result equal to zero, and raising one index on then 
yields 


d^xf^ 

ds^ 


4" 


dx°- dxA 
ds ds 


rJy^oc 

Fft^ 

ds 


(4.54) 
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CHAPTER 5 


Experimental Tests of General Relativity 


5.1 The Schwarzschild Solution 

Up to the present time, experimental tests of general 
relativity have consisted of the Eotvos experiment, the grav¬ 
itational red shift, the advance of the perihelion of planetary 
orbits, and the deflection of light by the sum. The last three 
of these effects are based on the solution of the field equations 
for a spherically symmetric field. We now obtain this solu¬ 
tion, following K. Schwarzschild (1). 

The starting point is the spherically symmetric squared- 
line element of flat space-time: 

ds^ ~ c^dt^ -}- dv'^ -|- v'^^dd^ -j- sin^ Q d(p^) (^*1) 

Now we introduce at the origin a mass with spherical 
symmetry. (5.1) must be modified but in a way that retains 
spherical symmetry and symmetry with respect to time re¬ 
versal. This leads to 

-ds^ = -/o(r', t)c^dt'^ -h h{r\ t)dr'^ 

-f ty^dd^ y sin 2 e dy) ^ ^ ^ 

We choose a new coordinate r such that f^{r\ ty^ — r^. We 
set /o = /i = e\ to obtain goo = — = e\ 

§33 = sin^ 6; then 

—ds^ = —e^c^dt^ -f e^dr^ -f r^{d6^ + sin^ 6 dcp^) (5.3) 

For the empty space surrounding the body, T == 0, 
and the field equations become 

Multiplying this by and contracting gives = 0. The 
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gravitational field equations for vacuum become 

Rfiv = 0 (5.4) 

To facilitate calculations we give the non vanishing com- 

= —sin 0 cos 0 

r\, = i 

or 

jr ®23 == cot 0 

jTigg = —r sin^ 0 e~^ 


For the vacuum case we use (5.4) or == 
set equal to zero. The expressions for G^^ are somewhat 
more simple and may also be used when a problem other than 
the vacuum with only gravitational fields is considered. The 
nonvanishing components of are given by straightfor¬ 
ward calculation as 


ponents of as 

dX 


= i 


12 


r\Q = I 

J.2 

ro 

10 

r\. 


dr 
c3t 

dx 


e 


1 
r 

X-~v 


cdt 

dv 

dr 

~re~^ 


(I — }l^\ — L 

Vr^ r dr) r^ 



fl! ^ 

rc dt 




(5.6) 
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Setting (5.6) equal to zero leads to the independent equations 

^ + 1 - = 0 (5.7) 

dr r r 


dr 



(5.8) 


^ = 0 (5.9) 

Bt 


The sum of the equations (5.7) and (5.8) gives 

— a-\- v)= 0 (5.10) 

dr 

(5.9) requires A to be independent of time, (5.10) then 
implies that any time dependence in v must arise from a 
function independent of r. (5.3) indicates that such a time 
dependence can be eliminated everywhere by a coordinate 
transformation involving only the time. This is equivalent 
to the statement that the assumption of spherical symmetry 
guarantees the possibility of a time-independent description 
of the geometry of the space. All time derivatives in (5.5) 
and (5.6) can be taken to be zero. The solution of (5.7) and 
(5.8) is 

e->^ = e” = \ (5.11) 

r 


The constant may be determined from the require¬ 
ment that Newton’s law of gravitation be approached, at 
large distances from the mass. 

From the geodesic equation it follows that the accelera¬ 
tion of a small test body at rest relative to the central 
mass M is 



Comparing this with the Newtonian value — GMjr^ gives 





(5.12) 
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The Schwarzschild metric (5.3) now takes the form 

(sin2 e d<p^ +d(92) -f. 


■ds^ 


— 






c^r 


(5.13) 

(5.13) exhibits^ the “Schwarzschild singularity” at 
r = 2GMjc^. This appears to limit the size of massive bodies 
to r > 2GMIc^. For an electron this is 13.2 x 10”®® cm, 
which is much too small to be related to present elementary 
particle theories. For a charged particle at the origin, the 
arguments which led to (5.3) are still valid. However, the 
energy associated with the electric field of the particle is 
distributed throughout space. The trace of the Maxwell 
stress tensor (4.15) vanishes. If we raise one index and then 
contract the field equations (4.25) the result is that R = 0 
for the combined Einstein-Maxwell field and we obtain 

r,, (5.14) 


for the case of energy sources which consist entirely of the 
electromagnetic field. The use of the metric (5.3), expressions 
(5.6) and (4.15), and the Maxwell equations (4.47) for spheri¬ 
cally symmetric § Aq then give in place of (5.13) for a particle 
with electric charge e 

O 


ds^ 


(• 


2GM . Ge^ 


H- 


^2 


dt^ 


r 

-f- r‘^(sin^ 6 dcp^ -f dO^) -f- 


dr^ 


(5.13a) 




2GM ,e^G\ 

■I" c4y2j 


C^T 


t H. P. Robertson has shown that a test body tfxkes a finite time to cross the 
Schwarzschild singularity. New coordinates remove the Schwarzschild singu¬ 
larity, but not the one at the origin, with singular Rieinann tensor invariants- 
Finklestein® has given an analytic extension of the Schwarzschild 
solution, with no singularities except at the origin. It does not possess time 
symmetry. He attributes this to the nonlinear character of general relativity. 

§ To evaluate the constant we may compare the invariant JF/**' at 
large distances with the flat-space value for a charged particle. 



60 


GENERAL RELATIVITY AND GRAVITATIONAL WAVES 


For a spherical uncharged body of finite extension the 
^5.13) is "valid only outside the body. A-n. interior 
solution for a fluid sphere can be written in the forna 


- A-BV ] 


c^dt^ 


-[- [sin^ 6 d<p^ dO^] 


(5.13b) 


1 — 




and by requiring this to properly join up with (5.13) at the 
boundary r = r-, we can obtain in terms of the density p, 

= {Sc^lSnpGy, A = (3/2) ]/ I 1 - : -B = i 

5.2 The Gravitational Red Shift 

We return now to the metric (5.13) to again deduce the 
gravitational red-shift formula. Suppose we have an oscilla¬ 
tor at the surface of a star, with interval between vibrations 
given by the invariant As, Imagine that a pulse of light is 
emitted at the peak of each vibration and let the time coor¬ 
dinate at any point in the space be established by the receipt 
of these pulses. First we note that the interval of coordinate 
time At, between successive pulses at a fixed point, is equal 
everywhere in the space to the value at the surface of the 
star. Each pulse is propagated along a null geodesic in the 
radial direction and (5.13) gives, for ds = 0, 

(5.15) 

egooir) 

Integrating (5.15) for the time of the arrival of the nth pulse 
at the radial distance r leads to 


t(r, n) — t(r„, n) = — j = R(r, Tq) (5.16) 

Similarly, 

i(r, n + 1) — t(ro, n + 1) = F(r, Tg) (5.17) 

Subtracting (5.16) from (5.17) gives 

l(r, n + 1) — t(r, n) = t(rg, n + 1) — t(r„, n) = At (5.18) 


*• gii(r')(ir' 

-0 cgoo(»'') 


F(r, ro) (5.16) 
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which is independent of r. Let the coordinate frequency be 

_ 1 _ c-\/“goo(distant oscillator) ___ 


I*^or the identical oscillator located at the position of an 
observer at radius r, the locally generated frequency is 
assumed to be given in terms of the invariant intervals Ashy 


= — V—goo (local oscillator) 


Comparing (5.19) and (5.20) results in 


(5.20) 


V. 


V, 




goo (distant oscillator) 


) 


goo (local oscillator) 

The frequency shift at radial coordinate r for 
which has its surface at radial coordinate Tq is 


(5.21) 
a body 


. GM / I l\ 

{5.21a) 

(toward the red), in agreement with (1.15), since r and are 
very nearly equal to the corresponding distances. The differ¬ 
ence between (5.21) and the equivalence-principle result 
(1.15) amounts to a fractional frequency shift of about 
(GM/cVo)2, roughly one part in lO^s for a terrestrial meas- 
■urement of the entire red-shift. 


Measurements on the companion of Sirius give a red- 
shift equivalent to the Doppler shift associated with a veloc¬ 
ity of 19 km/sec, while the calculated equivalent value is 
20 km/sec. For the binary star 40 Eridani, measurements on 
the white dwarf give 21 ± 4 km/sec, which is to be compared 
with the calculated value 17 Jz ^ km/sec. Line-shift measure¬ 
ments on the sun give different results at different points on 
the surface. This is not understood. 

Experimental work is in progress to measure the red 
shift by other methods, involving the earth’s field. One type 
of experiment compares precision frequency standards of the 
“atomic clock” (4) type on earth and at a different gravita¬ 
tional potential, using rocket-launched satellites (3). A dif- 
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ferent method makes use of the discovery of Mossbauer (5) 
that some of the gamma radiation (in the energy range 
ev) emitted by long-lived (greater than 10"® sec) “isomeric” 
states of nuclei is recoil-free. The recoil momentum is taken 
up by the solid as a whole, and there is no significant Doppler^ 


t The recoil momentum may excite the lattice vibrations or may be 
transferred as linear momentum to the entire crystal. In the former case 
there will be a Doppler breadth because the emitting atom recoils. Tlie 
criterion that the lattice vibrations not be excited can be stated in tlie 
following terms. Let the gamma-ray recoil momentum be p — fiwylc. Let 
the mass of the atom be m. Then we may expect no appreciable Doppler 
effect if (p^Jm) where is the Debye temperature and k is Bolte- 

mann’s constant. This result can be understood by consideration of the 
effect of the radiated momentum transfer. Suppose the lattice vibratioir 
wavefunction is y,- before emission of the gamma ray. •ipj is an eigenfuncdion 
of the Hamiltonian operator. If the gamma-ray emission transfers mo¬ 
mentum p in some way to the crystal, then the wavefunction immediately 
after emission is Here is the harmonic oscillator coordinate and 

this wavefunction is the one which has an expectation value for momentum 
p units greater than the initial one. After emission we may write 

a^ipr — e^'PaQ!‘lfi-y)^ 

ttr = Sr{+iPs\q%/fi' 

The term Ptlq’lrj gives rise to phonon emission or absorption while the terra. 
drj does not. To obtain the relative probability of emission without phonon 
excitation we need to evaluate 


modes 

k\ri are the harmonic oscillator matrix elements. 


or 


•Nph 


Nmoji ' ' NmcDi 

where is the number of phonons, N is the number of atoms, mV is t-bo 
mass of the crystal, and coj is the lattice vibration frequency. 


modes mcOj 


(■Nph +1) or 


h 

ma>i 


Nph 


Since energy is at least approximately conserved, to will be a high lattice 
vibration frequency, of the order of the Debye freciuency (o^ -= /cOp/fi, also 


V 


ph 1- Using these relations and requiring 

modes 


then gives (p^jm) k6j^. 
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shift, or broadening. The long lifetime is a consequence of 
the fact that the angular momentum of the nuclear ground 
state is different from that of the excited state. The radiation 
is either magnetic dipole or of a higher multipole order. The 
ratio of the frequency to line width is unusually great. Sug¬ 
gestions that this radiation be employed for terrestrial red- 
shift measurements were first published by Pound and 
Rebka (6), and independently put forth by Wilkinson, by 
Boyle, by Devons, and by Cranshaw (7). 

The properties of the 14,700-ev radiation of the lO-^-sec 
Fe^’, obtained from the radioactive decay of 270-day Co®’ 
have been reported. Pound and Rebka observed the radia¬ 
tion with a scintillation spectrometer and arranged to move 
the source toward and away from the spectrometer by ce¬ 
menting it to a moving coil transducer. In this way they ob¬ 
served that the radiation has a half width Av 2 x 10® 
cycles/sec. Since the red shift is proportional to frequency, 
a shift of roughly one part in 10^^ would be expected for a 
change in height of --^50 m in the earth’s gravitational field. 
For the Fe®’ radiation this amounts to about 3 x 10^ cycles/ 
sec, or roughly 1 percent of the line width already achieved. 
This is, therefore, a feasible experiment.^ Pound and Rebka 


t Results of such an experiment were recently reported by Cranshaw, 
Schiffer, and Whitehead {Phys. Rev. Letters 4, 103 (I960)). They employed 
and a total height difference of 12.5 m. A red shift 0.96±0.45 times the 
predicted value was observed. A more precise result was achieved later by 
Pound and Rebka, Phys. Rev. Letters 4, 837, (1960). They obtained a red 
shift 1.05 ± 0.10 times the predicted value. Anomalous frequency shifts as¬ 
sociated with particular sources and absorbers, and a temperature dependent 
frequency shift had to be taken into account. The latter effect is a time dilata¬ 
tion associated with the temperature dependent velocity. The angular fre¬ 
quency shift is 



ZcvykT 

I COy 


Here T is the absolute temperature, and the approximation on the far right 
is valid at high temperature.s (where e<iuipartition occurs). 

This effect was minimized by controlling source and absorber tempera¬ 
tures. The other anomaly was corrected for by exchanging positions of source 
and absorber. 
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note that with a source of limited strength, statistical fluctu¬ 
ations decrease the precision, owing to decreased counting 
rates, as the distance between source and absorber is in¬ 
creased. This just about compensates the increased red shift 
as the height separation between source and observer is 
changed. Larger vertical distances do not therefore give 
much increase in precision, unless some means is provided to 
concentrate the radiation. 


5.3 Effects on Planetary Orbits 

(5.5) may be employed to write the geodesic equations^ 
for either a massive particle or a light ray. The differential 
equation for the 6 coordinate is 


dH 

ds^ 


2 dr dd 
r ds ds 


— cos 6 sin 6 



0 


(5.22) 


For the particle with nonvanishing rest mass, s is the inter¬ 
val along the geodesic. The problem is simplified if the plane 
containing the initial velocity and the central mass is taken 
as 6 — 7 c/ 2 . dd/ds — 0 initially, and from (5.22) d'^d/ds^ — 0. 
The entire orbit will then lie in the equatorial plane. The 
remaining geodesic equations are 


d^r .dX/drY ~x i i ( * 

df 

=0 (5.23) 

dsj 

d^(p ^ 2 dr d(p ^ 

ds^ r ds ds 

(5.24) 

d^t dv dr dt 
ds^ dr ds ds 

(5.25) 

Equations (5.24) and (5.25) are readily integrated and give 

*2 _ T( 

d^~ ^ 

(5.26) 

dt 

==K 

ds “ 

(5.27) 


For these calculations the assumption is made tliat the coordinate 
system has axes fixed with reference to the stars. 
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and are constants of the orbit. Instead of inteffratinff 
(5.23) directly we obtain a first integral of it by writing the 
Schwarzschild metrie (5.13), for 6 = jr/2, in the forni 





(5.28) 


(5,28) is tli 6 rel<a»tion, between 9 .ny set of (quantities dT/dSy 
dcpjds, and dtjds in the space. If, however, we restrict dcold^t 
and dt/ds to the values given by (5.26) and (5.27), we can 
eliminate and dt and thus obtain an ordinary differential 
equation in terms of the variables r and 9 ? for the orbit given 
by 








— 1 


(5.29) 


We substitute W = 1 /r and differentiate (5.29) with respect 
to 93 , obtaining with the help of ( 5 . 11 ) 


d^W 

d(p^ 


-f- IT = 


GM 


ZGM 




(5.30) 


(5.30) can be solved in successive orders by an approximation 
procedure. Let p = 2GMIc^ be the Schwarzschild radius and 
let 


pW = pIF(i) + 4- . . . (5.31) 

where is a first-order and is a second-order quan¬ 

tity. Substitution of (5.31) into (5.30) and setting terms of 
corresponding order equal to zero separately yields the equa¬ 
tions 




-l_ pp'd) 


rf 2 q^( 2 ) 


d<p^ 

The solution of (5.32) is 


__ 11 ^( 2 ) 


GM 

3GM 

~c^ 




(5.32) 

(5.33) 


GM 

Ti/d) = (1 4- e sin (93 ~ 930 )) (5.34) 

; c 

Here e is the eccentricity of an elliptical orbit described by 
(5.34). We orient the coordinate system so that 93 ^ = 0 . Let 
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bt the leiig^th of the major axis of the ellipse (Fig. 4.1): 

An ex,.rej«ion for the constant {K,)^ can be written, using 
( j,t o), 111 ernis of the eccentricity and semimajor axis, as 

aGM 


{K^y 


(1 - e^) 


(5.36) 



Employing (5.36) and. (5.34) enables us to write the solution 
of (5.33), and the entire solution (to second order) of (5.30) 
is found to be 


W 


1 -|“ e sin (p 


SGMe 


a(X e^) 

GMb^ 
c^a{'l — 

In (5.37) the term 




c^a{\ — e^) 
GM 


(p cos Cp 


COS^ ^ - 1 - _ 

«2) c2a(l - ^2) 


SGM a (p cos g) 


(3 + 8^) 


(5.37) 


c2a2(l -£2)2 

will cause an advance of the perihelion. The second and 
terms can be expressed in the form 


e 


a{l . 8 ^) 

8 

a(l e^) 


Sin (p — 

1 “h 


SGM (p cos (p' 
c*^a(l — e^) _ 
9G^M^(p^ - 
c*a^l^ — 


i 


( 


sm 9 ? -f tan”^ — 


(5.38) 
SGM(p 


c^a{l-e^) 


) 
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(5.38) implies that for each revolution the perihelion will 
advance an amount A given by 

^ _ QtiGM 

~ 7^ a(X - (3-S9) 

The dependence of the perihelion advance on the eccen¬ 
tricity £ and semimajor axis a is evident. The quantity c in 
(5.89) is the velocity of propagation of gravitational inter¬ 
actions. Thus the agreement of (5.39) with observed values 
for the planet Mercury (astronomical observations give 
42.6 ± 0.9 see per century and expression (5.39) predicts 
43.0 sec per century) indicates that general relativity provides 
a good description of gravitation and that gravitational 
interactions propagate with the speed of light. 

It is also of interest to note that the failure of a planetary 
orbit to be exactly periodic provides an especially sensitive 
test ot the departures from an inverse-square law, for slow- 
moving bodies (v < c). The force on a body at rest is not 
the Newtonian value but rather 



GM 

^2 


GM- 

c^r 


(5.40) 


with r given in terms of tlie distance I from the surface of the 
central body as 



5.4 Deflection of Light 

To discuss the deflection of light in a gravitational field 
we must again solve the geodesic equations given by (5.22), 
(5.23), (5.24), and (5.25). The integrals (5.26) and (5.27) are 
still valid. Instead of (5.28) we obtain a corresponding rela¬ 
tion lor null intervals. First we set the Schwarzschild squar¬ 
ed-line element equal to zero, then we divide by 



S is now a parameter other than interval. Employing (5.26) 
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and (5.27) in (5.41) and letting W — 1/r leads to the equation 

It is easier to deal with the equation obtained by differ¬ 
entiating (5.42), which is 

d^W , BW^GM 

_+TF = _^ (5.43) 

(5.43) can be solved, to second order, using the method 
employed for (5.30); the result is 

1 cos 9 ? GM /...ox A A\ 

^ = 7- =-7- 9’) 

^ ^0 ^ '0 

In (5.44) the integration constant Vq is the distance of closest 
approach for the unperturbed light path (see Fig. 4.2). 





Fig. 4.2 

For r ^ 00 we obtain two values of <p which are solutions of 

GM . 


cos cp — — 


Tt 


(1 -1- sin^ <p) 


(5.45) 


Let the values be 


7^2 = 


7c dcp 

~2 ~¥ 

n d(p 


(5.46) 



experimental tests of general relativity 69 

Solving (5.45) gives for the total deflection 

d(p — ^GMjc^rQ (5.47) 

For light from a star just grazing the sun’s edge, (5.47) 
gives <p " 1.75 sec. 

A summary of the astronomical data on the red shift and 
deflection of light has been given by Trumpler (2). For the 
deflection as observed during eclipses, the individual meas¬ 
ured values have a probable error of roughly 10 percent and 
are usually within 10 percent of the value given by (5.47). 
The average for eleven eclipses agrees with (5.47) to roughly 
one part in 500. 

5.5 Concluding Remarks 

The four experiments discussed thus far constitute the 
entire experimental verification of general relativity. This 
is in marked contrast with the abundance of experimental 
data which support the quantum theory. It has been a 
great challenge to conceive other experiments which are 
feasible. The problem is to some extent a technological one. 
Techniques are usually not available to observe the very 
small effects which the very weak interactions of the theory 
predict. In later chapters additional experiments which were 
recently proposed will be discussed. 
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CHAPTER 6 


The Conservation Laws 


6.1 The Canonical Stress Energy Pseudotensor 

In special relativity the relation = 0, which is 

satisfied by the stress energy tensor, is integrated over the 
three space coordinates to get the conservation laws for 
energy and momentum. For we have 

^ r 

J ^0® dx^dx^dx^ 

udT^ dTg^ ar„3\ r 

+15^ + ~^) = - j To'dS, 

On the far right the volume integral has been transformed to 
a surface integral with dS^ s. vector element of area of the 
boundary. This result relates the change of energy in a given 
volume to the flow of energy momentum into it. 

The situation is different if a Lorentz metric is not em¬ 
ployed. The stress energy tensor now satisfies = 0. 

Since is symmetric we may use (3.79) to write 

( 6 . 2 ) 

If the metric is time-dependent then for /< = 0, the second 
term of (6.2) represents a time-dependent exchange of energy 
and momentum between the gravitational field in each vol¬ 
ume element, and the other fields. In order to be able al¬ 
ways to earry out the operations which led to the integration 
of (6.1), we search for a quantity t/ such that 

(—g)“*(VV—g),-= — 
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(6.3) 
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Then (6.2) becomes 


[(T/ +V)\/~^Xv = 0 (6.4) 


With (6.4), an integration similar to (6.1) may now in general 
be carried out which would be expected to yield conservation 
laws (1-5) including the gravitational field. Expression (6.3) 
does not uniquely define Also it follows from the right 
side of (6.3) that is not a tensor. We shall see that be¬ 
haves like a tensor under certain restricted groups of trans¬ 
formations. It is called the stress energy pseudotensor. One 
way of obtaining an expression for tf/ is to follow essentially 
the procedure which led to (4.10). The problem is different 
from that ordinarily encountered in field theory because the 
gravitational action function 


/(? = (c®/167rG) J Ry^ —g d^x (b-5) 

contains second space and time derivatives of the field vari- 
ables However, the part containing second derivatives 
may be written as a surface integral. We have, in fact, 

RV-i = V-g = 

The terms involving second derivatives are the first two on 
the right of (6.6). These may be rearranged in the following 


way: 

(g'‘''-r“^r.« - g'“'J"%...) V-g = g).» 

- V~ g).. - V- g), a+-r“^a(g'“’ V- g),. 

If we now use (3.73) 

3 g 


(6.7) 


9 ®“ 


gg'^g.K.a, g^'';«=0; 


rf‘yo,g'^''—r%^gi‘y; 


^g^'^ g IK, ft 

The last two terms of (6.7) may be written 
~R%K(g'‘W-g).<K + R‘‘^(g'‘’'V-g).K 

= v-g 


( 6 . 8 ) 
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Using (6.6), (6.7), and (6.8), we may express (6.5) in the 

form 

C , _ ( 6 - 9 ) 

Thus we see that in (6.9) the first integral on the right 
is only a function of g/**' and its first derivatives; the second 
integral can readily be converted into a surface integral using 
Gauss’s theorem. For variations of the action function which 
vanish on the boundary, only the first integral on the right 
side of (6.9) will be needed. We may now define 



and a total Lagrangian density by 

^ -|- Lp ( 6 . 11 ) 

with Lp the Lagrangian density of fields other than the 
gravitational field. Since (6.10) contains only the part of the 
gravitational Lagrangian without second derivatives, and 
Lp also contains only the field variables and their first deri¬ 
vatives, we may apply the standard Lagrangian formalism 
of field theory (section 4.1) to the action integral 

I = (6.12) 

It therefore follows thd.t the field e<^ua,tions &re 

— _ S(.Sf v-g) 

dx‘^ ^ Sg^*^ ^ (6.13) 

Using (6.11) and (4.37) we may write 

c* 1 _ dg/“' ~ dx^ Sg^ J 

(^..-Jg..i?)V-g = (5^)r,„v-g 


(6.14) 
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and using (4.10) we may define a total canonical stress energy 
pseudotensor by 


V v-g 


( 5 /^ V^g —gpfi^ ^ 


g(-gv-g) 1 


with a gravitational part 


(6.15) 


with 







G 


V— g — 




g( -^gV—g) ~ 


(6.16) 


(VV-g),v = ( V v-g), V = O (6.17) 

The expressions (6.13) and (6.14) are especially interesting 
because they exhibit Einstein’s field equations in the form 
of Lagrange’s equations. 

It is more instructive to derive (6.15) in a different way 
making use of invariance of under certain coordinate trans¬ 
formations. We recall that the conservation laws of a theory 
are associated with certain invariance properties. For exam¬ 
ple, in mechanics energy conservation is related to invariance 
of the Hamiltonian under time translation, angular-momen¬ 
tum conservation is related to invariance of the Hamiltonian 
under rotations. 

Let us carry out the infinitesimal coordinate transforma¬ 
tion* 

x'P' = H- SxP (6.18) 


making use of the transformation law for gives for the 
change corresponding to (6.18), 


dg^P' 


gKOC 


d B 

= g''‘/‘{x')—g'=/‘{x) (6.19) 


t In this chapter the <5 symbol always means what is expressed by (6.18) 
and (6.19). Some authors employ the symbol <5 for the different operation 

'SgK/i = g'Kfi{x) ~ = dxKifi + (5aj/t;/c 
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^(^-20) 




dse dse ^ 


This follows because Lf is a scalar, so dLp = 0. 
T\Tfl.Ving use of (6.19) and (6.20) leads to 


( 6 . 21 ) 


6£e. 




g/CV _|_ 2 




6 ,« 


dse 


(< 5 a:/‘),„ 


( 6 . 22 ) 


We may write 

3:Sf 




1 a(.sfv-g) 


5-. = [(^V-g)(-gH] = 


16.231 


Here we again made use of Ag = —gg^„ dg'"’'. We note that 
Y'—g does not contain g’‘i‘,f- Employing (6.23) in (6.22) gives 


(6.S?)V-g=[2 

a(.s?v-g) 




^gKll 




fi 




g'^. 




(d^^).v + 2 ^^''^g'^-- - - g"-(gfl;^)..,« (6.24) 


c^is invariant under linear coordinate transformations; for 
these, the last term on the right of (6.24) is zero, and (5«^ 
must vanish. The remaining terms vanish, giving 


2- (-^y g) g,„ 2 - - ■- g''',^ 




ag« 




8(^V-g) 
dg*'’... 


(6.23) 


0K/i 


0 


(6.25) does not contain SxM'; it is an identity. It then follows 
that for a general coordinate transformation, (6.24) leads to 
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^(-y/ — 

id^)V-g=2 ( 6 . 26 ) 

For any coordinate transformation which changes noth¬ 
ing on the boundary of the region of integration, the sur¬ 
face terms are zero and we write 


+ ( 6 . 27 ) 

The left side of (6.27) vanishes because the integral is a 
scalar, so for the right side 

== {dSfs/^gd*x +j^5{-V^d*x) = 0 ( 6 . 28 ) 

The second integral on the right of (6.28) is zero because 

V —g is a scalar. The vanishing of the first integral on the 

right together with (6.26) leads to 

's/ — g 


/ 


Sg''*. a 

We may write (6.29) in the form 


( 6 . 29 ) 




d^V- 




^g-1 ) 

- ,v/ 


d^x = 0 
(6.30) 

The first term in (6.30) can be converted into a surface in¬ 
tegral which again vanishes for with derivatives which 
approach zero sufficiently fast on the boundary. The re¬ 
maining integral in (6.30) may again be transformed into 
a part containing a surface integral which vanishes and 
a remaining integral which for arbitrary Sxf^ gi'ves 

rd^V~g 


L 4 ’"*, a 


<fKV 

O 


0 


(6.31) 


J, V, a 


(6.31) gives us the conserved quantity 

/d^V 


T/V' 


•g 


(r/ + v)V-g = -2 




-Jg,A ( 6 . 82 ) 

,a /,« 




76 


GENERAL RELATIVITY AND GRAVITATIONAL WAVES 


written as an ordinary divergence. (6.32) is very useful for 
calculations. 

We may now again obtain the forms (6.15) and (6.16) by 
multiplying the field equations (6.18) by and writing the 
first term as the difference of two terms to obtain 




dscy-g 

J 


a 


8g’‘>‘. a 


g 


KV 


g 


KV 


d^V- 


■g 




0 (6.33) 


Employing (6.25), differentiating with respect to and 
making use of (6.31) leads directly to the expressions (6.15), 
and (6.16) for //. 

For convenience in carrying out calculations we give 
expressions for 


—g 
r 


and 


d^G\/—g 


dgf^^ 


First we evaluate 


The use of (3.73) and (3.74) gives 


(6.34) 


SgP’^.y 


igpiy 


(6.35) 


Noting that g^f.^ — 0 = setting 

P = and employing (3.74) leads to 


3(g'“'-r»„.) 


+ (5<,»a^r) + Jg-rgp^ (6.36) 


8g<“'.y 

With the use of these relations it follows that 


a(g'‘-i'’%.rv) 

ag'’". y 


-igi‘'’rr^,gp. 


+ r^^l-iid^z-d^y + d^i^d^y) + igi^yg^^'] 


Also, the method used to obtain (6.36) gives 


(6.37) 
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v-g) 

dgp^^y ptr S 

From (6.10), (6.37), and (6.38) it follows that 
d(.^Q\/ g) 

== 1 — J. ^ -i- y .... — 00.7 rp „\a 

pt 


(6.38) 


dgp^^ 


-- {- + 4(g/‘>'rr^„ - g‘^-'r»af)gp. 

+ V—g\ 

The use of the field equation (6.14) allows us to write 


(6.39) 


d^Q's/—g 
dgp(r 


(Rp<r igpcrR) (i6^g V — g) + 


V — § 

. -1.7 

(6.40) 


Evaluation of (6.40) is facilitated^ by the fact that the 
result cannot contain derivatives such as since ^'sZ—g 

does not contain such derivatives, so these terms will cancel 
each other on the right of (6.40). We obtain 


/167rG\ d^G ‘\/—^ 


act 


dgf'^ 

H- r^p^ + gp. 

StzG 

gpcr^G “f" a/ § 


(6.41) 


We return to a consideration of the properties of the 
canonical stress energy pseudotensor r/ of (6.15). Following 
(6.1) we write 

3^5 J V V'— g = — J V g dSi (6.42) 

For an isolated system, for whieh the approach the 
Lorentz metric sufficiently fast, and for which T is local¬ 
ized, the surface integral on the right of (6.42) vanishes and 


t This follows Moller; it is more simple than direct calculation. The 
variables are gP^ and gpv may be expressed in terms of but gpv.a. 

can only be expressed in terms of and 
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the integrals 

P>,=^!r/V^gd^cc 


(6.4,3) 


are independent of time, (t is regular everywhere.) It fol¬ 
lows from (6.82) that systems in which 

8{^V-8) 

0KV —- 1 - 

^ a 

has the same value over a closed surface will have the same 
P^. Under Lorentz transformations transforms as a four- 
tensor. Pn will transform like the energy-momentum four 
vector of special relativity (13). 

We now calculate —Pofor a Schwarzschild particle and. 
show that it is Me, To do this we first carry out the coor¬ 
dinate transformation 


=/( 


GM\^ 


l[(- 


2GM7-\i 


+ r 


GM- 


(6.44) 


The coordinates r', 6, 93 , and t are called isotropic coordinates- 
The Schwarzschild metric (5.13) becomes 



If we now introduce 


+ r'^dd^ -1“ r'^ sin^ 6 dep^) 



(6.45 ) 


-{-y^ + = r '2 

X ~ r' cos q) sin d 
y = r' sin (p sin Q 
z — r' cos 6 

and assume that r' is very great, ( 6 . 45 ) may be written 
~ds^^ (1 + {dx^-\-dy^~\-dz^) — ^l ~ cW (0 45 ) 
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We note that all quantities are regular everywhere since a 
Schwarzschild interior solution is valid inside the body 
producing the field. 

The use of (6.46), (6.39), (6.43), and (6.32) then leads 
to 

= Me = (Pi = P^ = P3 == 0) 

^ (6.47) 

Pp = 


On the other hand, if we had used the Schwarzschild 
metric (5.13) directly we would have obtained 


Po=“-oo (6.48) 

Also if a flat-space metric 

— ds^ — dr^ -|- r^(d6^ sin^ 6 dep^) — c^dt^ (6.49) 

is employed an infinite result is again obtained for the 
energy. These results are not surprising because is not 
a tensor. 

The conservation laws (6.15) and (6.16) are seen to give 
energy-momentum four vectors under quite restricted con¬ 
ditions. The Lorentz metric must be approached sufficiently 
rapidly at spatial infinity and the special metric (6.46) has to 
be used. No consistent interpretation of as an energy 
density is possible since only the volume integral is meaning¬ 
ful. Coordinate transformations affecting the interior of a 
region will change the values of within it, but not the 
volume integral as calculated in a coordinate system for 
which (6.46) is appropriate. The energy is simply the coef¬ 
ficient of 2G/r'c^in (6.46), so all systems which have the same 
constants appearing in (6.46) will have the same energy 
regardless of the interior description. 


6.2 Other Conservation Laws 

A conservation law such as (6.17) states that something 
is conserved; in a given system of coordinates the conserved 
quantity may not be related to the energy density in a simple 
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way. The expression (6.16) is only one of a class of pseudo- 
tensor densities which satisfy conservation laws. Addition of 
any quantity with a vanishing ordinary divergence to 
will not affect (6.17). This arbitrariness in the 
specification of has been employed by Landau and Lif- 
shitz (6) to obtain a conserved symmetric pseudotensor 
and by M0ller (8) to obtain a quantity whose zero-zero 
component transforms like a scalar under coordinate trans¬ 
formations not involving the time. In order to obtain these 
results it is convenient to introduce (7) a quantity U^ 
which is antisymmetric in a and v. The circumflex symbol 
is often used to indicate this property. The canonical stress 
energy pseudotensor may be written in terms of U as (6.50). 

T/ ^/~-g = a- (6.50 ) 

satisfies the conservation law {r/\/—g)^y — 0 inconsequence 
of the antisymmetry of 17/^. We may search for an expres¬ 
sion^ for U by studying (6.32) and noting that it must 
differ from 


cr 


gccv 


by some quantity whose ordinary divergence vanishes. If 
the quantity is written the vanishing of its divergence 

is guaranteed. We therefore write 


~ 2 ( 6 . 51 ) 

(6.39) enables us to evaluate 

+ ^ V-g (6.52) 


^ The are sometimes called the superpotentials. 
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Remembering that 

r-fia = g.^/2g, = 0 and , 

rearranging terms then allows us to write (6.52) in the form 






actv 


i( V^Vg/5- — v^g a/g-'^),^ 


16 jrG 


(6.5S) 


Comparing (6.53) and (6.51) we see that 

"" i GnG (§■'“- g”'“g''^)].^ (6.54) 

^ ^ ihc ^ S"^ - S''”]). f (6.55) 


These quantities transform like tensor densities under linear 
(affine) transformations, and are called affine tensor densities. 

For formulation of angular-momentum conservation 
laws it is convenient (but not necessary) (4) to have a quan¬ 
tity corresponding to which is symmetric in the indices a 
and p. Such a quantity can be obtained from the in the 
following way. Consider 

+ ( -g)”i l-g{g”>‘g‘''‘ - g'^g"'’)],^, J (6.56) 

The last term in (6.56) is symmetric in ^ and v, but the first 
is not. If tlie factor (—g)“^ in front of (6.54) is removed, the 
modified quantity 

A/*- = l-gig'^g-P - g''>‘g”^)lp (^^) (6.57) 

will have a divergence symmetric in v and ju. We use to 
indicate the Landau-Lifshitz gravitational pseudotensor and 
iC/**' for their combined pseudotensor density. Landau and 
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Lifshitz employ (6.57) to write 

J^nv (j'/iv ^ = h^KP p (6.58) 

Kp'^ then satisfies a conservation law in consequence of the 
antisymmetry of in v and p. Integration yields a quantity 
given by 

= L r (_g) [j/.o ^ l^o^dx^dx^dx^ (6.59) 

C J 

Under linear transformations transforms as a four- 
vector density rather than as a four vector. Making use of 
the field equations and the definition of Landau and 
Lifshitz give the following formula for 

+ g’^gy‘‘{rp^prpy^+rpy^rpfp - rp^^rp^y - r^^yPp^p) 

+ g^rg-p{rp^,rpyp - PPfyr^^p)} 

Again Kp^ has the undesirable properties that to obtain the 
required result a special coordinate system is required in 
which the isotropic metric is approached at large distances. 
Again no consistent interpretation of as localized energy 
density is possible. 

To remedy these defects, M0ller noted that we may add 
another third-rank affine tensor density V to the right side 
of (6.51) without affecting the validity of the conservation 
laws, as long as is antisymmetric in v and a. M0ller sug¬ 
gested that we define a function X/r^ by 

(6.6i) 

Then a stress energy pseudotensor A/ is given by 

^/V-g = [(L/ + t/)V~g] = (b.62) 

Here i/ is the gravitational pseudotensor corresponding to 
The additional affine tensor density V is now sought 
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which will give the following properties. For physical 
systems for which it is possible to use the “Lorentz” metric 
(6.46) at large distances 

— g(Px (6.63) 

should behave like a four vector under all coordinate 
transformations which do not involve the time coordinate. 
This results in Aq^ being a scalar and Aq^ being a three vector 
under arbitrary spatial transformations. In order to find a 
V satisfying these requirements M0ller investigated the 
transformation properties of —g under arbitrary in¬ 

finitesimal transformations which do not involve the time 
coordinate. This indicates in what way 
satisfy the requirements of a scalar density under spatial 
coordinate transformations. Combinations of and their 
first derivatives, which for (6,46) lead to which fall off 

faster than 


.2 ’ 


for r —> oo 


must be considered. In this manner M0ller was led to pro¬ 
pose the function 

== 4- (6.64) 

Using (6.54), expressions for the functions and 

may be calculated. 

(6-65) 

In order to calculate A^^ we will need 

Xo- = ^ 3 ^ 1 " - (goA « - go., ^) g-" (6-66) 

We study now the behavior of (6.66) and Aq^ under the 
group of coordinate transformations which do not involve 
anywhere the time coordinate, but are otherwise arbi¬ 
trary. goa transforms like a four vector; recalling that for any 
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vector B — By. ^ — ^v,(i then tells us that the 

bracket in (6.66) transforms like a tensor. It follows then 
from (3.78) that —g transforms like a scalar 

density, xo''^ transforms like an antisymmetric tensor density, 
so it follows that Aq^ transforms like a three vector. Thus we 
see that M0ller’s tensor, (6.62) and (6.65), does have the 
transformation properties which were sought. 

Now we shall consider briefly the association of conserv¬ 
ed quantities with infinitesimal coordinate transformations, 
following Peter G. Bergmann (9). Consider an infinitesimal 
coordinate transformation 


x'f*- = xi* bxi^ (6.67) 

The quantity Mp defined by 

Mp = ( 6 . 68 ) 

satisfies Mp^^ = 0. Since (3<r« is arbitrary, it is clear that 
(6.68) represents an infinite number of such laws. If we 

choose (50?* to be a set of constants, the Einstein formulation 

is obtained. If we choose — g, where is a 

set of constants, the Landau-Lifshitz expressions result. 
Angular-momentum conservation laws may be obtained by 
selecting 

— (g<^^xy — (6.69) 

with Jpy again a set of constants. It is evident that any vec¬ 
tor field may be employed in (6.68). 

Komar (10) has shown how to construct a set of conser¬ 
vation laws in tensor form. Starting with M0ller’s result he 
writes 


Dp^-g = {dx-u^e + 

~ "I 

= gj.« 

Addition of the curl field 

g,.] V-gj.. 


(6.70) 


WP-s/—g= 


(6.71) 
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to (6.70) then gives 

~ (6.72) 

(6.72) is conserved from (3.75), and it is in tensor form. 
A study of (6.70) and (6.71) shows that for (6.72) 

gives the result of M0ller. 


6.3 Further Remarks on the Conservation Laws 


It is clear from the foregoing that in situations where 
the matter-energy is sufficiently localized that a ‘‘Lorentz” 
metric (6.46) may be employed at large distances, the total 
energy is well defined and calculable using several of the 
expressions given here. The pseudotensor proposed by 
M0ller does not require the metric (6.46), and the zero-zero 
component is the same in all coordinate systems which have 
the same time coordinates. The Landau-Lifshitz result is a 
symmetric quantity. 

Many other stress energy pseudotensors may be con¬ 
structed, and Dirac has noted the greater possible utility of 
particular ones in solving the equations of motion. Other 
expressions suggested by radiation problems are given later in 
this tract. In the chapter on radiation we shall note some 
other difficulties which lead to the conclusion that the prob¬ 
lem of energy in general relativity has not been solved in a 
completely satisfactory way. 

Peter G. Bergmann makes a distinction between con¬ 
servation laws which are identities and those which are valid 


if the field equations are satisfied. The identities such as the 


Bianchi identities, or = 0 are called “strong” laws. 

Relations which are valid only if the field equations are satis¬ 
fied are called “weak” laws. 


It is interesting that a fourth-rank tensor, a completely 
symmetric quantity, has been discovered by Bel (11). To 


construct this wc need first the concept of a dual tensor. If a 
given tensor is A its dual may be constructed by writing 


(6.73) 
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here is again the Levi Civita tensor density. The dual 
of the Riemann tensor is 


^(Xp yS 's/ §1^ 

(6.74) 

Bel’s tensor is 

'^a.pyS ^<x[ivS y ~\~ ^ ^a/ivd^^p^^y 

(6.75) 

and satisfies the relations 


(6.76) 

= 0 

(6.77) 


provided satisfies the field equations. These properties 
were established by I. Robinson in as yet unpublished work. 

In addition is symmetric in all indices. The “con¬ 
servation” law (6.76) does not lend itself to construction of 
the three-dimensional integrals used earlier. Bel’s tensor is 
being studied in order to utilize its properties in the problems 
of gravitational radiation. 
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CHAPTER 7 



7.1 Weak-Field Solutions 

One of the central issues in the general theory of relativ¬ 
ity has always been the question of gravitational radiation. 
Gravitational waves have thus far not been observed. Until 
recently (18) no solutions of the exact gravitational field 
equations were known which represent spherical gravitation¬ 
al waves. For these reasons a great deal of theoretical work 
has appeared on this subject during the past four decades. 
We shall see later that some experimental work now appears 
feasible. Some theoretical issues have been resolved in recent 
years, and it has been possible for a number of physicists to 
conclude that general relativity really does predict the exist¬ 
ence of gravitational waves. 

In 1916 Einstein (1) studied the weak-field solutions of 
the field equations 

SjtG 

~ ""^ 4 “ Tfiv (^* 1 ) 

obtained by assuming 

giMv == V + (7.2) 

is the Lorcntz metric and h^y is a first-order quantity. 

and h are defined by 

V = (7.3) 

h = (7.4) 

The Ricci tensor is given in terms of the Christoff el 
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symbols as 

The use of (7.2), (7.3), and (7.4) enables us to write to 
first order as 


The last three terms of (7.6) can be arranged in the form 
h^fiv — (7.7) 

(7.7) can be made to vanish by choosing the coordinate 
conditions such that 

(V - iVX/? = 0 (7.8) 

jR^j, will now consist only of the first term of (7.6), and 
the field equations may be written as 


= — T,, (7.9) 

Care must be exercised in the use of (7.9), since all compo¬ 
nents on the right will not ordinarily be of the same order. 
We now define 

9^/ = h/ — Id/h (7.10) 

The coordinate conditions (7.8) are written as 


9 /, I. = O / .ii; 

Raising the index v in (7.9) and employing (7.10) gives 

l^jiG 


U9/ 


T 


(7.12) 


The solutions of (7.12) are well known from electrodynamics, 
as 

C(T/)ret^rdeddV 


(p/{r, t) 


P 


(7.13) 


Here ]r r [ is the distance from source point to field 
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point. As we noted earlier, the relations (7.12) can be ob¬ 
tained directly from Newton’s law of gravitation, in which 
is replaced by □; the conservation law ^ 0 then 

yields (7.11). 

Some years ago Pauli and Fierz (3) considered the 
question, what relativistic wave equations would be appro¬ 
priate for particles of zero rest mass and spin two. Now in a 
relativistic theory for spin S, 2(2S -f 1) components are 
needed, so a second-rank tensor is required. They were led to 
propose the equations 

Uw/ = 0 (' 7 '- 14 ) 

with supplementary conditions 

= 0 (7.15) 

Since these equations are the same as (7.11) and (7.12) 
for vacuum it is expected that gravitons will have spin two. 
Since the only available direction is the direction of motion, it 
follows that the spin angular momentum must be oriented 
in the direction of motion. Since the gravitational forces 
have infinite range, it follows that the rest mass of the 
gravitons must be zero.* 

We consider now the conditions under which the ex¬ 
pression (7.11) may be satisfied. Suppose we carry out an 
infinitesimal coordinate transformation 



= 03“ -f |“(a!) 

(7.16) 

To a first 

approximation this can be written 



03“ = X'°^ — ^“(03') 

(7.17) 

The transformed metric tensor is then given by 



^ IX.V ^jJLV Boev (1 Bfifl V 

(7.18) 

Similarly 

^ fiv fi 

(7.19) 


t If the forces are transmitted by particles with rest mass w, a Yukawa 
potential <p ^ (e-wicr/ft/y) jy appropriate. This has infinite range only if 
m 0. 
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(7-20) 

<pV“« = ?'/*“«— <5'"" f A, (’'•21) 

If the coordinate conditions (7.11) are not satisfied, we 
have 

(7.22) 

(7.21) states that an infinitesimal coordinate transfor¬ 
mation may be carried out to make 

(7.23) 

provided 

□ ('^•24) 

On the other hand, if 9 ?^“,« ~ then we may carry out 
arbitrary coordinate transformations which maintain this 
condition, provided that the functions satisfy the wave 
equation, according to (7.24). 

We return to consideration of (7.12). In free space these 
equations represent waves with cp/ given by 

□ <^/ = 0 (7.25) 

It is of interest to inquire how many independent components 
of or (p^^, are needed to describe a plane-wave solution of 
(7.25). We will find it convenient to discuss first the case of 
a locally plane disturbance of arbitrary strength, then to 
apply the result to the weak-field case. 

7.2 Riemann Tensor for a Wave of Arbitrary Strength 
Which Is Locally Plane 

Suppose we have a gravitational wave which has van¬ 
ishing space derivatives in directions normal to the direction 
of propagation. Such a wave is therefore locally plane- In 
this case we shall show that at any given point all compo¬ 
nents of the Riemann tensor may be written in terms of the 
derivatives of the three components of which describe 
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intervals in the plane transverse to the direction of propaga¬ 
tion. Introduce a geodesic coordinate system so that the 
Christoffel symbols vanish at the given space-time point and 
the metric has the Lorentz values there; the Riemann 
tensor becomes 


^aPyd 4" §,ay, fiS — ^ fid, ay) (7'.26) 

Let the direction be the direction of propagation; then 
the partial derivatives of and with respect to and 
vanish. A study of (7.26) indicates that we may write the 
components of R^fiys in the three groups 


-^2020 = 

i^22,00 

i?122o ~ 

^^22,10 

■^1212 — 

4^22,11 

^2030 ~ 

i^23,00 

-^1230 ~ 

i^^23,10 

■^1213 = 

fe23,ll 

-^3030 = 

i^33,00 

•^1320 “ 
^1330 = 

4^23,10 

fe33,10 

■^1313 = 

fesaii 



-^1223 = ^ 

•^2323 

= 0 


i^l323 — ^ 

^2320 

= 0 


■^1023 ~ ^ 

^2330 

= 0 

^1030 — 

“fel3,OoH~fe30,10 

■^1020 ~ 

~tel2,00 + ^20.10 

^1310 “ 

“ 1^30,11+ fel3,10 

^1210 ~ 

“~fe20,ll+fei2.10 


-^1010 — ~feii,oo + ^ 10 , 10 “ i^oo,11 


(7.27) 


(7.28) 


(7.29) 


In vacuum = 0; at the given point the use of the 
Lorentz metric in conjunction with (7.27), (7.28), and (7.29) 
gives 


•^12 — 

•^1020 

= 0 


^13 — 

■^1030 = ^ 


i^20 ~ 

i^l210 

= 0 


II 

■^1310 ~ ^ 


^11 = 

•^1010 

"“■^1313 

4^1212 — ^ 

•^10 = 

-^1220 "f ^1330 = ^ 

(7.30) 

1^22 

■^2020 

7? 

^1212 ~ 

= 0 


i'^2030 -^1213 === ^ 


J?33 == 

-^3030 


= 0 

^00 ~ 

-^3030 + -^2020 H“ ^1010 

= 0 
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The above equations for the vanishing of R^^^ 

and Rqq imply the vanishing of the first four values of 
given in (7.29). The equations for the vanishing of R-^j^, 

J^ 33 , and i^oo combined to solve for the last expression, 

-^loioj in (7.29), which turns out to he zero. It follows there¬ 
fore that the only nonvanishing components of the Riemann 
tensor are the ten components (7.27), and these are given 
entirely in terms of g 22 » and ggg. The expressions (7.80) 
for R 22 , Rz 3 , and Rqq also lead to 

^ 22.00 + ^ 33.00 = 0 (7.31) 

(7.31) states that the required second derivatives of the 
two components and ggg are not independent. In the 
weak-field approximation, (7.26) is valid everywhere. At all 
points the Riemann tensor can therefore be written in terms 
of A 22 J ^ 33 , and 7^23, subject to (7.31). It follows then that 
coordinate transformations may be carried out to transform 
away” all components of except the “transverse” ones, 
everywhere. This can be done in a relatively straightfor¬ 
ward way, using infinitesimal coordinate transformations 
satisfying (7.24). In the next chapter we shall show that 
relative displacements of closely spaced free particles in the 
path of the wave occur only in directions transverse to the 
direction of propagation of the wave. 

7.3 Approximate Evaluation of Source Volume In¬ 
tegrals 

We return to a consideration of the weak-field solution 
integral expressions (7.13). These components may all be 
written in terms of in a certain level of approximation. 
The stress energy tensor satisfies, to a first approximation, 
the conservation laws 

^0k,7c ^ 00,0 ^ (7.82) 

^ 30,0 ~ (7.33) 

(7.33) can be multiplied by cc^ and integrated (by parts) over 
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all space, neglecting surface terms which vanish at infinity. 
This is written in the symmetrized form 

J (7.34) 

(7.32) can be multiplied by integrated over all space, 

again neglecting surface terms which vanish at infinity; the 
result is 


[ J = — J [Tiaa^+ T,qX''\ d?x (7.35) 

Expressions (7.34) and (7.35) enable us to write 

J Tijd^x— J TqqX'^x^ (7.36) 

The use of (7.36) facilitates certain calculations in which 
the distribution of mass and energy is more readily dealt 
with than the distribution of stress. The relation (7.36) has 
neglected retardation effects and is not valid where a source 
large compared to a wavelength is involved. 

For a source in a region which is small compared to a 
wavelength, (7.36) may be used to calculate all the space 
index components of in the lowest nonvanishing approx¬ 
imation. The use of the coordinate conditions then allows 
determination of the remaining components. In some cases 
this is a simpler procedure than direct use of (7.13). 

For a particle it is necessary to evaluate the integrals in 
the limit as the volume tends to zero, taking account of the 
fact that retardation effects imply that different volume 
elements are taken at different times. It is more direct to 
start with the result for a particle at rest, then to generalize 
it for a moving particle, as in electrodynamics. A particle of 
mass w which is at rest has 

4iG7n 

9^oo(r) = ^2|i. — r'l 

In the linear approximation we may infer that the generali- 
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zation for a moving particle must be 

, , 4^GmV\U'^ 

<pA^> t) = _ r») 

where now U\ is the four velocity of the particle, m is the 
rest mass, and 

— (r'i _ c{t' t); t' == t — \r — r'|/c 
(7.13a) may be written in the form 



7.4 Weak-Field Approximations for Energy Flux and 
Energy Density 

The covariant divergence of the stress energy tensor 
vanishes, in consequence of the field equations and the 
Bianchi identities. This may be expressed using (3.79), in the 
form 

V-g).. - igcA = 0 (7.37) 

The stress energy pseudotensor density ti/\/ — g is re- 
lated to the stress energy tensor density by 

+ {t/ V~^),V — 


0 


(7.38) 
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tf/ may be calculated using the results of Chapter 6, but 
for weak fields it is more direct to compare (7.37) and (7.38) 
and write 


{t/ '\/~§>),v — —V — i ^ (7.39) 

The form of (7.39) suggests that we multiply (7.12) by ^ 
to obtain 


□ 9a/f = 






(7.40) 


Utilizing (7.10) and (7.39) then leads to 

B2nG 

(V V-g),v 

We may write 

= [A'"’, , -hhj.h''+ Ih, h-p 5 /], „ 

Comparing (7.42) and (7.41) gives 


(7.42) 


VV-g = V”-AA->'+i5/(4fe ^ 

(7.48) 

Contraction of (7.10) results in 

= —K”- (7.44) 

Using (7.44) and 7.10 gives us the alternate form 




(7.45) 


7.5 Linear Mass Quadrupole Oscillator 

Consider now a linear mass quadrupole oscillator at the 
origin of a rectangular system of coordinates. The motion of 
the masses is along the axis. Let 1 be the peak value of the 
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time-dependent part of the moment of inertia. We can use 
(7.36) to write the solutions of (7.12) in the form 

2a>^GI cjo 

V 33 = ^ cos — [*» — ?•] (7.46) 


Oscillations in the direction can only give rise to 
the components 9933 , 9930 , and 9930 . The coordinate conditions 
90 ^^^ = 0 can be employed to calculate 993 ^ and 99^3 if we know 
that the coordinate conditions have been met. If the energy 
density in the vicinity of the sources is not so great that the 
space is appreciable curved, the Lorentz metric conservation 
law = 0 will be valid, to a good approximation. If the 

sources are in a small region so that retardation effects are 
unimportant, it then follows from (7.12) that the coordinate 
conditions are in fact being met. We may then write 

0 = 0 (7.47) 

9 >o ®,3 + 9>o‘’.o = 0 (7.48) 

Let the observer be situated along a radial line making 
an angle d with the axis. (7.46), (7.48) and (7.47) give 

9oo = e 9933 ; 9933 = —cos e 9933 (7.49) 

We could have calculated 99 Q 0 and 9933 directly from (7.12) 
or (7.13a); this would have required more careful considera¬ 
tion of retardation effects, and the energy in parts of the 
oscillator which furnish restoring forces. 

The radiated power can be calculated in the following 
way. We need 

— JroO-s/—(7.50) 


By use of the conservation laws, (7.50) may be written 

c 8 r ,— __ _ _ r ,— 


IV—g — c J 


gTo'dS, 


t^<Px -k c 


J V—gt^‘dSi 


(7.51) 
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The second and fourth integrals are over a closed (two- 
dimensional) surface. The source is assumed to be in the 
vicinity of the origin, and if our bounding surface is far from 
the source, the second integral in (7.51) will vanish because 
To' is zero outside the source. The third integral will be 
periodic in time while the last integral will give the steady 
decrease of energy of the source. Making use of (7.46), (7.45), 
and (7.49) leads (15,17) to the formula for the time-averaged 
radiated power Pi 


GPco^ 

15c^ 


(7.52) 


7.6 Radiation from a Spinning Rod 

Historically the first gravitational radiation problem to 
be considered was the radiation from a spinning rod (1, 2). 
This will now be ealculated. Assume the rod has a mass per 
unit length <y and spins in the xy plane with the angular 
velocity w. At t=0 the rod lies along the x axis. Then again 
using (7.36) and (7.13), we obtain 


<Pii 


c^7' 


2 

ar^ cos'^- dr- 

c 

retarded 


2co 

cos— {x^—r) (7.53) 


C^T 


2G 


7*22 


7i2 


c‘^r 




err" - dr 


4tGI ^co‘ 


c^r 


retarded 


2G 


|2 r 


9 o na I ■- 

c^r 8x^ J \ C 

retarded 


. / OJX^ 


\ I 

cos - 

/ \ c ; 


dr- 


2a) 

cos —(cT® —r) (7.54) 


‘^GI^ , 2ct) 

sin—(a?’^ 




(7.55) 


In (7.53), (7.54), and (7.55) is the moment of inertia 
of the rod. The coordinate conditions (7.11) are again em¬ 
ployed to calculate and and the radiated power is 

calculated using (7.45); the result is 

^ S2GIJoj^ 


5c® 


(7.56) 
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7.7 Further Remarks on the Weak-Field Solutions 

It is of interest to note that (7.56) and (7.52) can be 
obtained if we assume only that gravitational effects are 
propagated with the speed of light and follow the kind of 
heuristic argument* which has been used to calculate the 
radiation from an accelerated charge. 

Expression (7.36) states that the lowest-order multipole 
radiation is quadrupole radiation. This follows also from the 
fact that for an isolated oscillating system the dipole moment 
vanishes in consequence of conservation of linear momentum. 
For example, suppose we have a small mass which is coupled 
to a large mass by a spring (Fig. 7.1) and the system is 



Figure 7.1 

oscillating. If and Xm are the displacements of the small 
mass and the large mass, respectively, we have from momen¬ 
tum conservation 

mx^ + Mxm = 0 (7.57) 

and 

mx^ + Mxm = 0 (7.58) 

The left side of (7.58) would be expected to give the 
dipole contribution to the (accelerated mass) radiation, but 
it vanishes in consequence of (7.57). 

Bonnor (17) has calculated the loss of mass energy of a 
quadrupole oscillator in the second approximation and found 
it equal to the radiated energy given by (7.52) in the first 
approximation. 

t See, for example, Richtmyer and Kennard, Introduction to Modern 
Physicsf McGraw-Hill, New York, Fourth Edition, 1947, p. 62. 
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7.8 Exact Cylindrical Wave Solutions 

Einstein and Rosen (4, 5) obtained some exact cylin¬ 
drical wave solutions of the equations of general relativity. 
Their metric is 

— ds^ = — c^dt^) H- p^e~^'^d(p^ + e^'^dz^ (7.59) 

The functions y and 'ip are dependent upon p and t alone. The 
metric defined by (7.59) can be employed in the free-space 
equations 

= 0 (7.60) 

to give the following (exact) equations which have to be satis¬ 
fied by yj and y: 


1 

v.p.p + — y>.p — p 'P.t.t = 0 

(7.61) 

y.p = pli.v,p)^ + (i/c^)(v,«)^] 

(7.62) 

y,t = ^p'v.p'pjp 

(7.63) 


(7.61) is a linear equation, the cylindrical coordinate 
wave equation. A solution of (7.61) representing outward 
traveling waves is 

y) = A IJoicop/c) cos cot + Nq{copIc) sin cot] (7.64) 
Making use of (7.62) and (7.66) then gives 
y = (A^copj2c){jQ{coplc)jQ{o)plc) A-Nq{cop/c)Nq{copIc) 

+ {coplc)[{Jo{<oplc)y + Uoicoplc))^ + {N,{a>plc)y 

(-^0 (7.65) 

+ [J(,{oiplc)Joio>plc)—No(<oplc)Nf,'((oplc)] cos2cot 

+ [Jo{(^ple)N„'(a>plc) + No(coplc)Jo'((oplc)] sin Hmt) 

— (2j7i)A^(ot 

The prime means differentiation with respect to {(oplc).^ The 
last term of (7.65) increases linearly with time and at first it 
was thought that this represented the cumulative effect on the 
metric of continuous radiation of energy. Rosen argued (7) 
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against this interpretation on the grounds that if the energy 
is being radiated by a cylindrical system near the origin, 
then the loss of energy would make it impossible for the 
function y, given by (7.64), to remain periodic in time. The 
occurrence of the function Nq makes (7.65) singular on the 
axis. Marder remarks that it is not possible to join the free- 
space solution (7.65) to the solution for a region near the 
origin containing mass in a physically acceptable way since 
the sign of the density will ultimately change. He considers 
the difficulties associated with (7.65) as due to the fact that 
the field has been turned on for an infinitely long time. 

A more interesting case is that of a pulse of radiation. 
The linearity of (7.61) enables this to be treated, employing 
the Fourier integral. The following discussion is a summary 
of a paper by Weber and Wheeler (16). For the Fourier 
transform of yj we choose 

y}{co) = Jq{(opIc)Ic (7.66) 

This has the virtue that the integrations required to con¬ 
struct yj(t) and y(t) can readily be carried out. 

f OO 

^ Jo cos cot dcojc 

^ (7.67) 

= ^[([« — + p2)-^ +([a + ictY + P^)~^] 

Integration of (7.62) and (7.63) then gives 
y = p^[(a—ictY p^]-^ 

A study of (7.67) shows that for negative values of i, the 
‘pulse is imploding, onto the axis. It explodes back out 
again for positive values of t. It is a “time symmetric” solution 
of Einstein’s field equations. 

In order to calculate the “energy density” we employ 
the expression (6.82): 
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Direct calculation using (7.59) and (6.39) leads to 


^(^oV-g) 


= 0 


500 S 


^^22 , oq ^^33 


~ ^ V dx^)\ 

With these (7.69) becomes, for rQ^^Z—g 


-\-g^^ 


dx^ I \167i:G 


(7.70) 


V—g 


, d / 

To“A/-g = Yp\~^ + 


l) =0 (7.71) 


This calculation, (7.71), is valid for any metric of the form 
(7.59) whether matter is present or not. The everywhere- 
vanishing energy density, (7.17), is confirmed by use of the 
Landau-Lifshitz forms (8) (6.58) and (6.60) for The use 
of M0ller’s tensor, (6.62) and (6.65), gives a somewhat similar 
result. Direct calculation gives for the metric (7.59) 


z^oV—g 


4}7zG 




While this expression does not vanish identically, the inte¬ 
grated value for the energy per unit length. 


A^’>V—gdpd<f> = ^ 



oo 

(7.73) 

0 


does vanish for (7.67) and (7.68), (7.72) is positive at some 
places and negative at others for the pulse-wave solution. 
We note that the pulse, (7.67) and (7,68), represents at / = 0 
a space which is asymptotically flat at large p but which is 
not “Euclidean,” since according to (7.68), y at 
large distances. The metric is therefore “conical” at large p, 
for ^ = 0. 

Rosen has recently (20) calculated Tq® and for the 
cylindrical solutions in a “cartesian” coordinate system and 
finds them to he finite and reasonable. His result is not in 
conflict with that given by M0ller’s value for the energy 
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density since the passage from the cylindrical to Rosen’s 
cartesian coordinate system presumably requires a trans¬ 
formation involving the time. The M 0 ller value is only in¬ 
variant under spatial coordinate transformations. In the cy¬ 
lindrical coordinate system there is no guarantee that the 
“energy” will always be positive definite. If negative values 
are allowed, then the vanishing of Tq® may mean only that 
absorption of energy from the wave makes the wave energy 
negative. A lower bound for the energy would guarantee 
stabilitv. 

7.9 Interaction of a Particle with Cylindrical Gravi¬ 
tational Waves 

Calculation of the components of the Riemann tensor for 
the metric (7.59) shows that the components do not all vanish 
for the cylindrical waves, so the waves must be real. In this 
section we show that the cylindrical waves do carry energy, 
by analyzing the motion of a particle initially at rest which 
interacts with the cylindrical pulse. We write the geodesic 
equation for p, and integrate it in successive orders, regarding 
y; as a first-order quantity and y as a second-order one. Let 

P — P(o) P(i) + P(2) ■+•••• (7.74) 

From the metric (7.59), 

C^=l + vV2 + . . . (■^)72 + . . . (7.75) 

These relations lead to the first-order equation 

dt^ ^ dp (7.76) 

= --BcV(o){[(a-M<)"+P(o)T* + [(a+ici)2+p<o,^r*} 

and the second-order equation 

^V(2) __ ^P(n (dy\ , 


(7.77) 
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In obtaining these equations we have made use of 

Integration of (7.76) gives, for the boundary condition 
dpjdt = 0 at i = —• oo, 

d/ 0 (i) icB a — ict _ cl + ict 

dt P(o) _[(a— p(Q)^]^ {_(cL~\-ictyp(o)^]^. 

(7.78) 

The change in p from t = — ooto35 = 0 may be written 

z)p=|‘’ J‘ (^jd«'di = ^[(/5(„,^ + a")*-«]->2B (7.79) 

for large p. The distance from the axis, at i = 0, for large p, 
is 

^ jp 0 H** ^ 

Igii^dp = ey-'^dp ^Po^ 2B-2B In (2p(o)/a) (7.80) 

The results (7.78), (7.79), and (7.80) indicate that to first 
order, the particle, initially at rest, starts to move when the 
wave arrives; the change in the coordinate p at the instant of 
maximum implosion of the wave is -j- 2J5, but at that instant 
the change in distance is 2 B — 2 B In {2p^Q)|a). As the wave 
explodes back out again these steps are reversibly retraced 
in this approximation so that at i = + oo fbe particle ends up 
where it started and at rest. However, one could imagine the 
particle to be coupled to a second, distant particle by a 
mechanism giving irreversibility.^ The relative motion 
would then lead to absorption of energy from the wave even 
in the first approximation. 

To consider the second approximation we need to inte¬ 
grate (7.77). Carrying out the indicated operations leads to 

t In the next chapter 'we show that energy may always be absorbed from 
a system for which the Fourier transform of does not vanish, in a local 
Lorentz frame propagated along a geodesic. 
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t.^nty^wo integrals. We calculate only the velocity at 

. „ three integrals can readily be shown 

to vanish, by contour integral arguments, and the result is 


/ f ill >-|-“0o 

4i3’^ 


P(0) 


J oo 

-CO 


c^dt 


+ 


T 


cH^dt 


Pio) J-oo {[(“+*c<)^+/)(„)2][(a—irf)24.p,||,a]}| 

J_co [(a + icty + P((„2]* 

These integrals may be written in terms of 

p1Tl2 

F^j 

Jo 


de 


[1 — ^ sin2 e]i ’ 


P(0) 


as 


H~ P(o) 


(7.82) 


(: 


cdtj 


+ ^2 


B^F 

dF 


Lao) +• />(0)^)^ + 


-4_ -^(0) , ^Pi0)(5^^+P(0)^) 

I / . ‘"5 1“ ^ = --- 




lp(0>ia^+Pio)^)i ■■(a2+7,”2)i' 

gp,o>(5a^+6aV,„,2+p,„,4) 

(“ +PioA)i 


12a2 


(«^+P(o)“)5 

P( 0) ^ )P( 0) 


d^F 

dz^ 


“f- 


3^(0)(«^+/)(0)^)~^ 4-B2 

gp(0)(5a^4-6a2p(Q^24-p^^^4)^^2__^^ ^^ 2 ) 


/>(0)(^"-f P(Q)")^(a"~y0(0)2)2 

“f" /^(o)^)"^ 


(“ +P(0)^)* ^ _ 


(7.83) 
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This lengthy expression reduces, for the limits small and 
P(o) lai'ge t~^co to: 

dp^2) 27zB^cp(Q) 


y .«ij (’■•85) 

(7.84) and (7.85) indicate that the particle is left with a 
residual velocity^ relative to the axis after the wave has left 
and the interior space is again flat. The axis has a real signif¬ 
icance in this problem because a test particle could be placed 
there; according to (7.79) it would not move at any time 
thereafter. Observing a particle at radial coordinate p there¬ 
fore amounts to observing the relative motion of a pair of 
particles. 

Conclusions similar to these have been reached inde¬ 
pendently by L. Harder (18). He extends the pulse solution 
to the interior of a mass cylinder which is the radiator. After 
emission of the pulse he finds that the mass of the cylinder 
has decreased. 


dt 

( 

a 

dp(2) 

2 B^c 

v -. 

dt 





(7.84) 


7.10 Exact Plane-Wave Solutions^^ 

Plane gravitational waves have been considered by a 
number of authors. Taub (6) and McVittie (11) showed that 

t This velocity produces infinite displacement, over an infinite time. 
This does not affect the valitlity of our series expansion, which is in powers 
of Bja. 

it These waves arc not i>lane, since for propagation in the x direction 
the departures from flatness dei)end on y and z. Bondi calls them plane 
because they have as much symmetry as plane electromagnetic waves. 
These metrics admit a group of motions with five parameters. 

An illuminating discussion of these solutions has been given by Bonnor 
(24) and is summarized below. Gonsider the metric 

-d.v“ = Ad^^-\-Bdy^-\-Cdz^~-2bdydz~~Ddrp 
Here A, B, C, b, and D arc functions of alone. 


(footnote confd) 
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there were no unpolarized plane waves. Robinson and, later, 
Bondi were able to show (13, 14) that the field equations 
permit ‘‘plane” wave zones of finite extent between two flat 
space-time regions. They start with a metric originally given 
by Rosen: 

ds^ = — de) — u^{e^^dri^ -h (7.86) 

with u = r — and Q are functions of u alone: = 

This metric satisfies the vacuum field equations 
^fiv == If represents, in general, a curved space. If the 


Choose coordinates for which A ~ D; then let 
obtaining 

—— ~dxdT-]rBdy^-{-Cdz^ — 2bdydz 

Here B, C, and b are functions of x alone. 

For this metric the vacuum field equations are 

dx^ . 2x \dx/ \dxl dx dx ^ 

and X — BC—b^. The other vacuum field equations are satisfied without 
further restriction on B, C, and b. We have then three unknown functions 
with only one relation which they satisfy. Two of the three may be arbitrarily 
chosen. Take b == 0 and let B = C — q^; then the field equation for 
Rji becomes 

1 d^p 1 d^q 

p dx^ q dx^ ^ 

Take as a particular solution of this equation p ~ sin nx, q — sinh nx, where 
n is a real constant. This solution appears singular on the hypersurfaces: 
X = mnjn, a; == :£ co, where m is an integer. 

These singularities may be removed by carrying out the following trans¬ 
formation 

03 ® = x+ny^ sin nx cos nx-\-nz^ sinh nx cosh nx 

X^ ~ X 

— y sin nx 
x^ — z sinh nx 

and when this is carried out the metric becomes 

— ds^ = — dx^ dx^-{•dx^^-j-dx^^ — n^dx^^(x^^ — 

The singularities at a? == mnin have disappeared. The singularities at infinity 
may also be removed, by choosing a new origin. 
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relation 

Up,u,u -h 2/5,, = J3 (7.87) 

is satisfied, all components of R^fiys vanish and the space is 
flat. By carrying out the coordinate transformation 

u = T—^ = ct—XI = ct-{-x — 

rjue^ = ?/, u^e~^ = z 

for u :> 0, the metric (7.86) becomes 

ds^ — c^dt^ — dx^ — dy^ — dz^ 

+ 2/5,,[( 2 /%— zdz) {cdt — dx) — ( 2 /^— z^) {cdt — dx)^(ct —a?)"^] 

— —00 ''^){cdt—dxY (7.88) 

For /5,, = 0 except for a finite positive range of u, there 
is a curved region of space between two flat regions. The 
amplitude of the waves is determined by /5, which is an 
arbitrary function of u. 

A second wave type is described by the metric 

ds^ = e^^{dx'^ — d^'^) — (t—|^)^ [{cosh 2/5)(d97^-|-dC^) 

+ sinh2/5cos 26(dr]^ — dC^) — 2 sinh 2^5 sin 2ddyd^'] 
here again Q, /5, and 0 depend on u = r — with 

2i2,, = (r - ^) -f 0 ,,2 sinh2 (2^)] (7.90) 

Bondi notes that at the boundaries the conditions of 
Lichnerowicz on the continuity of and its first derivatives 
can be met. (These conditions are given later in this 
chapter.) 

7.11 Initial-Value Formulation of the Radiation* 
Problem 

A different approach to gravitational waves is to note 
that the field equations may be regarded as setting up require¬ 
ments on the initial values** of the fields,* then predicting 
what the future evolution will be. A study of the acceptable 

* Such investigations were first carried out by Lichnerowicz (10) and 
Foures-Bruhat (27). 

** This section is a partial summary of the thesis of D. Brill, Princeton 
University, 1959. 
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initial values then gives information concerning possible 
wave solutions (22, 23). 

The corresponding situation for the Maxwell equations 
is well known. For no charges or currents the equations 


V*F=0; = 0 (7.91) 

constitute conditions on E and H at say ^ = 0. The remain¬ 
ing equations 


Vx F = 


i 

c dt ’ 




then give the time evolution, preserving (7.91) at all stages. 

For the general relativity equations we may choose the 
spacelike surface a:® = 0. The equations contain second-time 
derivatives, so the initial values will be the and their 
first-time derivatives. By use of the contracted Bianchi 
identities, we write for the Einstein tensor 


( 7 . 93 ) 

(7.94) 

No higher than second-time derivatives occur on the 
right side of (7.93) and (7.94), so it follows that and G/ 
cannot contain higher than first-time derivatives of the 
The stress energy tensor will not contain derivatives of 
It follows, therefore, that the four equations 

R/^o _ ^g/^OR ^ jr^o ( 7 - 95 ) 

c 

or 




(7.96) 


may be taken as initial value equations, and that the six 
equations 

~ r« 


(7.97) 
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or 

SjtCr 

Ri, - igijR = T,, (7.98) 

are the equations of time evolution. 

We can now show that if the set (7.95) or (7.96) is satis¬ 
fied at the initial time, then any solution generated by inte¬ 
gration of (7.97) or (7-98) will continue to satisfy the initial 
value equations (7.95) or (7.96). Let 


X/iv = Xfcv _ 


SnG 


Xjti' 


(7.99) 


Since the covariant divergence of T vanishes, we may 
write, from (7.94), 

» V ^ 

For fi = k, (7.100) becomes 

= + n^,X‘‘’‘-\ 

and for fj, = 0 (7.100) gives 

= - -h H- 

Consider first expression (7.101). At ^ = 0, Xf^^ vanishes 
everywhere so the space derivatives X^\^ also vanish. There¬ 
fore we may conclude that X^^^q vanishes initially. Now 
consider (7.102); since Xf^^' and X^^ vanish everywhere at 
t = Oy X^\i also vanishes; therefore A*®®,o vanishes initially. 
Now differentiate (7.101) partially with i*espect to time. It 
can then be shown that initially vanishes. If we differ¬ 

entiate (7.102) partially with respect to time, then the vanish¬ 
ing of X^^Q everywhere initially guarantees the vanishing 
of which then guarantees the vanishing of 

the initial time. By continuing this process it then becomes 
clear that X^^ and all its time derivatives vanish initially. 
The Taylor expansion of therefore guarantees that X^^ 
vanishes for all time. 


(7.100) 

(7.101) 

(7.102) 
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7.12 Time-Symmetric Solution with Positive Energy 


A tensor T is called time-symmetric if 

x^) = x^) (7.103) 

where n is the number of times that zero occurs amongst 
a^yd . - . 

We shall see that introduction of this element of sym¬ 
metry substantially simplifies the problem of solving the set 
(7.95) and (7.96). A time-symmetric solution should give 
the same value for ds^ for a pair of events with the same coor¬ 
dinate differentials, at and — cc^, so = 0; also all 

tensors in the problem must have the same kind of symmetry. 
It follows, therefore, that the equations J?®* — = 

are automatically satisfied since g®*, and 
vanish. The initial-value relations reduce to the one equation 


7? 0 
-^0 




/StzGX „ ^ 

hr) 


(7.104) 


(7.104) may be written in the forms 

+ R/) = KRio^" - 


= -iRm‘ 


ik 


h 


8:7rG\ 

4 


) 


T 0 

J. 0 


(7.105) 


In a time-orthogonal system of coordinates reduces to 
the curvature scalar computed from the three-dimen¬ 
sional metric. This is because Since i 

and k are space indices and g^^ = 0, it follows that only 
sums of space derivatives of the space index g^^ enter into 
The initial-value equations become 


(3) R + (^^) r„<> = 0 (7.106) 

In vacuum, for “pure” gravitational waves, (7.106) is 

= 0 (7.107) 

Consider the solutions of (7.107). These allowed initial 
values should show the essential characteristics of the com¬ 
plete solution. Brill obtains information about the solutions 
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in the following way. Let the required metric be satis¬ 
fying = 0. Let ds'^ be obtained from another metric 

ds^, according to 

ds'^ = e^^ds^ (7.108) 


The unprimed is called the base metric. P is a function 
of the coordinates. The transformation (7.108) is a new kind 
of operation, for unlike the coordinate transformations used 
thus far, intervals do change. At a particular point all differ¬ 
ential intervals are changed by a factor e^. It follows that 
the shape of a small triangle or any small figure will not 
change. Transformations (7.108) are therefore called con¬ 
formal. It can be shown that the value of P' in terms of (the 
base metric) R may be written (25) 

R' = e-^^[R + (w-l)(2V2P (^_2)(VP)2)] (7.109) 


In (7.109), n is the number of dimensions of the space 

V^P = (VP)^ = g-P,P, 

In our case we are concerned with R' for the three space, 
so n = 3. It is convenient to make the substitution 


P = ^((n-2)/2)P 


(7.110) 


(7.109) and (7.108) become, for R' — 0, 

ds'2 PW. -h ^ ^ 0 (7.111) 

A careful study of (7.111) by Brill has led him to the 
conclusion that solutions of. (7.111) exist with P well behaved 
and not vanishing anywhere,^ corresponding to a metric 
with the asymptotic behavior 

( O \ 

1 + j (dr^ + r^Q^) + goocW (7.112) 

The function W may be expanded asymptotically in a 
series A + Bir-\ where i = 1, 2, 3, . . . A study of (7.111) 


t The second equation of (7.111) has the form of the curved space 
Schroedinger equation, for which such solutions are known. 
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shows that the first nonvanishing term in the expansion of R 
goes as r (7.112) provides an unambiguous way to specifv 
the rnass energy of any wave for which the expansion in 
negative powers of r is appropriate, so that (7.112) applies. 
Brill finds that this “Schwarzschild” mass is a positive 
efinite quantity for fields having axial symmetry*. To 
show this, the second equation in (7.111) is divided by S' 
and integrated to obtain ^ 




Since 



(7.113) 




1 -f- 


J 


2Gm 
c^r 
Vyj 

W 


and 


yi 1 _}_ 


d^x 


Gm 

2cV 


dS —— 27zGmlc^ 


(7.1U) 


From (7.113) and (7.114) we have 
27tGm f/Vv^\2 p 

0^ (7.115) 

writer symmetric case the metric may be 

^ dz^) p^dcp^] ~ (7.116) 

1 /^ 2 . asyrnptotically at least as fast as 

’p ~J u on the z axis, 

inteial f o®, ‘^'•^'^ket of (7.116)) the 

Ls1?rnnrr If ^P^ere centered at tie 

p.m) LThf^Sr*’’' '* 


m 


inG 




d^x 


(7.117) 


field approilStio^ withmit''th **** positive definite, in the weak 

PP oximanon. without the assumption of axial symmetry. 
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Since ^has no zeros anywhere, the “Schwarzschildian” 
mass given by (7.117) is well defined and is a positive definite 
quantity for a localized disturbance with axial symmetry 
which implodes for negative t and explodes for positive t. 
The conservation laws guarantee the invariance of m before 
and after the moment of time symmetry. Furthermore it is 
apparent from (7.117) and (7.112) that the mass is zero only 
if the space is flat. 

7.13 Conditions on the Differentiability and Continu¬ 
ity of Manifolds 

Until comparatively recently it was customary to re¬ 
quire that the entire space-time continuum be covered by one 
nonsingular coordinate system. This has turned out to be 
unnecessarily restrictive. More than one coordinate system 
may be employed, provided that certain conditions can be 
met in the region where they overlap. A region in which a 
particular coordinate system is employed is sometimes called 
a coordinate patch. 

Lichnerowicz (10) has carefully studied the conditions 
which have to be met in order to ensure that the gravitational 
field equations have unique solutions for the case = 0. 
He requires that at the intersection of two coordinate patches 
the local coordinates of a point in one of the coordinate 
systems be fourfold-differentiable with respect to the coor¬ 
dinates in the other system and have nonvanishing Jacobian. 
The first and second partial derivatives must be continuous; 
the third and fourth partial derivatives must be at least 
piecewise-continuous. The metric tensor must be continuous, 
and is required to have continuous first partial derivatives; 
the second and third partial derivatives must be at least 
piecewise-continuous. These conditions on the metric tensor 
are required everywhere. 

7.14 Six-Dimensional Treatment of Gravitational 
Radiation 

Pirani (12) has given an interesting discussion of gravi- 
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tational radiation. This is based on the work of Petrov (9) 
and the generalization of ideas familiar in the treatment of 
electromagnetic radiation. For the vacuum case, 

Rjiv = 0 (7.118) 

First consider a gravitational wave characterized by a dis¬ 
continuity of the Riemann tensor across the wavefront. 
The Lichnerowicz continuity conditions on the metric tensor 
and its partial derivatives ensure uniqueness. Suppose we 
have a given point in a coordinate system such that is a 
Lorentz metric at the point, the Christoffel symbols vanish, 
and the wavefront discontinuity there is moving in the 
direction. Let the new coordinates ^ and C be defined by 

^ = (^0 — a7i)/^2 (7.119) 

t == (ir® H- x^)|^/2 (7.120) 

The line element at the point becomes 

ds^ — 2didC — dx^^ — dx^^ ( 7 . 121 ) 

The surface of discontinuity is d^ = 0. Then and its 
first derivatives must be continuous, so it follows that 
niust also be continuous. However, 


( 7 . 122 ) 

will not necessarily be continuous. A study of the expression 
for the covariant Rieniann tensor Rcc^ys in the manner of 
ection 7.2 shows that in this special coordinate system the 
discontinuity in the Riemann tensor can be written in terms 

and alone. Only two variables, 

O’ and (p, with 




(7.123) 


are needed to describe it. The terms in (p may be made to 
vanish by an appropriate choice of axes in the 23 plane. 
A six-dimensional formalism proves convenient for the 
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classification of the discontinuities in the Riemann tensor and 
for other purposes. Coordinates are again chosen such that 
the six-dimensional space has a pseudo-Euclidean metric at 
the point, given by — (I 5 I 5 1? —Ij ——!)• The 
physical components of a given tensor at a point are defined 
to be those measured by an observer in a locally Lorentz 
frame. If are the physical components of an antisym¬ 

metric tensor, the corresponding six vector is obtained by 
relabeling the suffixes in accordance with the rule 


ajS 23 31 12 10 20 30 

^ 1 2 3 4 5 6 


(7.124) 


The physical components of the Riemann tensor Ra/sys 
related to the symmetric six tensor Rab> ^-Iso according to the 
above prescription, taking pairs a/S and yd. The discontinuity 
of the Riemann tensor may then be written 


ARab-= 


0 0 

0 —a 

0 — g) 

0 0 

0 — 9 ? 

0 (T 


0 0 

—(p 0 

(T 0 

0 0 

cr 0 

(p 0 


0 0 

—(p a 

a (p 

0 0 

a (p 

cp —a 


(7.125) 


This result followed from the Lichnerowicz conditions. 
Comparison of (7.125) with Section 7.2 shows what was to 
be expected, namely, that if we have a wave which is a dis¬ 
continuity over a plane propagating in the direction, the 
nonvanishing components of the discontinuity in R^pys 
just those of (7.27). 

Some aspects of energy flow in electrodynamics are now 
reviewed for the purpose of generalizing them for the gravita¬ 
tional field. An observer at rest in a Lorentz frame has a 
Poynting vector where is the electromagnetic stress 
tensor. Let a Poynting-like vector be given by Pp. Then a 
covariant expression for Pp , valid for an observer with four 
velocity Up may be written in the form 


(7.126) 
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Observers who move with the wave velocity observe no 
energy flux, so for them (7.126) vanishes. This leads to 


(7.127) 


(7.127) states that an observer who measures no energy flux 
has velocity U'^, which is an eigenvector of the electromagnetic 
stress tensor T If the electromagnetic field is a radiation 
field, the field tensor satisfies the relations 0 and 

F^^ Fr^ = 0. Such a field is called a null field. A study 
of the electromagnetic stress tensor indicates that if the field 
is a null field, has one null eigenvector corresponding 
to the eigenvalue zero. A radiation field, in electrodynamics, 
is, therefore, a field whose corresponding stress energy tensor 
has a vanishing eigenvalue. Since all four velocity vectors 
are unit vectors, it follows that the Poynting vector cannot 
be transformed away. 

To consider somewhat similar ideas for the gravitational 
case, the eigenbivectors of the Riemann tensor are de¬ 
fined by 


R.„P^ = AP. 


(7.128) 


is antisymmetric and is written as aj?t, — 

The unit vectors and bj, define a two dimensional space. A 
second two space which is completely orthogonal to the first 
is also defined by P^^. The term completely orthogonal 
means that every vector in the first two space is orthogonal 
to every vector in the second two space. The pairs of two 
spaces associated with different eigenbivectors may intersect 
each other to define at each intersection point a four vector 
tangent to the intersection curve. These four vectors are 
called Riemann principal vectors. If the Riemann principal 
vector is a null vector, the gravitational field is said, by 
Pirani, to be a radiation field. 

Petrov (9) has shown that by suitable orientation of the 
coordinate system in which the metric has the Lorentz form 
at a given point, the Riemann tensor may be reduced to a 
canonical form of one of the following kinds: 
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Type I: 

oci 0 0 0 

0 a2 0 0/^2 

_ 0 0 a3 0 0 

J<AB = 0 0 — ai 0 

0 ^^2 0 0 -a2 

0 0 /53 0 0 

Type II: 

— 2a 0 0 —2/3 0 0 

0 a —a O 0 /5 <7 

0 0 a-f-cr 0 a P 

— 2/3 0 O 2a 0 0 

0 /3 cr 0 O’—a 0 

0 O' /3 0 0 —a —a 

Type III: 

0 —a 0 0 0a 

-cr 0 0 0 0 0 

0 0 0 <T 0 0 

0 0 cr 0 O' 0 

0 0 0 o 0 0 

o 0 0 0 0 0 

The a and p are scalar invariants of the Riemann tensor. 
The value of o depends on the orientation of axes in the 10 
plane. ^ 

A study of these types shows that type I has one time¬ 
like, and the three spacelike, principal vectors; type II has 
one null principal vector and two spacelike ones. Type III 
has only one principal vector and it is null. Rirani’s criterion 
is, therefore, that gravitational radiation is present if the 
Riemann tensor is of type II (null) or III, but not if it is 
type I. 

t 3 X 3 matrices can also be vised, in. consequence of the obvious sym- 
rneti-y. If tlic iqiper left 3 X 3 is called M and the upper right 3 X 3 is called 
N, the matrices M H- iN give the essential features. 



(7.131) 



(7.130) 
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If at each point in a space in which the metric has the 
Schwarzschild value we introduce a Lorentz frame with axes 
parallel to the coordinate axes, the Riemann tensor turns out 
to be of type I above with 

—iai = aa = as = GMjc^r^-, /?* = 0 

For the Einstein-Rosen metric a similar procedure leads to 
type II with 

O' = Iv.p.t + + W.tV.p — ^f.pY.t'lh 

It must be emphasized that the components of 
must be calculated in a locally Lorentz frame at any given 
point, in order to make use of the classification scheme of 
Petrov. 

A special coordinate system which has a Lorentz metric 
in the vicinity of one point, vanishing first derivatives of 
and specified second derivatives of is the “normal” 
coordinate system. At the pole the following conditions are 
met: 

xf^ — 0 


(7.132) 


= g^„, p == 0 

Spv, pir ^{^p/cvtr ” 1 “ ^pvpcr') 

Pirani has calculated the average value of the canonical 
stress energy pseudotensor density 

t/V-g = V^eV-g - g 


pa-, ft 


8g 


per, V 


about a point in a frame described by (7.132). The average is 
defined by 


W—g t/y = (47rr4)-i r d^S (7.133) 

lim r—j -0 


This is done to provide support for the definition of radiation 
given here, t/ can be written in terms of so it vanishes 

at the pole of any geodesic coordinate system. The average 
defined above will not vanish since the second derivatives of 
can contribute. Unfortunately this average does not 
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have the dimensions of energy density. It is an energy like 
object which gives (section 8.2) a measure of the energy 
which may be abstracted by an observer with four velocity 
for whom (7.132) is appropriate. The average is calculated 
to be 

(7.134) 

For an observer at rest with type I Riemann tensor this 
gives (summations over k) 

w-gt/y = (^) [2{i<x,^(vv + 2v^/ 

and for type II, 

<V--gt/> = (*) [(x2(-42d/4-16V V + 22a/(5/ 

+ 22(5/V) -i- 4oc(T(d/^2>'—d/V) (7.136) 

+ Sc^d/ + d/){d,’'-d,>')-] (^) 

7.15 Other Petrov Glass II Wave Solutions 

The tetrad (sometimes called quadruped, tetrapod, vier- 
bein, or four nuple) is a set of four unit vectors which may be 
introduced at any |3oint to define a locally Lorentz frame.^ 
Three of the vectors are spacelike and normal to each other 

t The transformation to a locally Ijorentz frame with coordinates a?'* 
may be obtained by writing 


\16uzG/ 






doc.'^ 


A(a)‘'A(/?)v ~ 


for 
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and to the time direction which specifies the fourth one. The 
tetrad vectors are usually written where fi is the vector 

index and (a) labels a particular one. 

Peres has obtained another class of wave solutions of 
R[iv — 0, which have less symmetry than the plane and 
cylindrical waves. His metric is 

— ds'^==dcc^-\-dy^-\-dz^-\-2F{x^ y, z-\-r){dz -f- dr)^ — dr^ (7.137) 

jR-TT ^,xx H“ ^, vv ^ 

harmonic function. 


The field equations become —R^^. ~ 
and are therefore satisfied if jp is a 
A tetrad of orthonormal vectors is introduced by 

(1 + F, 0, 0, -F) 

(0, cos a, sin a, 0) 


^(0)^ 

^( 1 )^ 


(7.138) 


(0, —sin a, cos a, 0) 

(F, 0, 0, l-F) 

with tan 2a = F^/^.a!®* The components of Rapys which 
do not vanish are (in a frame defined by the reference tetrad) 

O' = Rf^Xvl — ^/i2,v<2. — {F,xx^ “b ^,xv^)^ (7.139) 

for y and v equal to 3 and 0 only. This class of metrics there¬ 
fore falls in Petrov’s class II. 

For 

F = {x^ ~ y^) sin {z r); o = 2 sin {z -f r) (7.140) 

we have a plane wave. A wave packet solution is 

F = a;2/(a?2-[-^2)-2ei:*'"~(*+T)2]-2^ for jsj + t| < 6 

F = 0, for I 2 : + t| > 6 

We now show that all wave solutions of R^^ = 0 of 
physical interest approach type II. Any physically interesting 
wave solution is likely to be generated by localized distribu¬ 
tions of matter, and will therefore be plane at large distances 
(not cosmologically large distances). It follows that our 
analysis of Section 7.2 describes it. We emphasize that our 
definition of a locally plane wave is one which has == 


(7.141) 
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where i and 7 are the two tetrad directions normal 
to the direction of propagation of the wave. The ten com¬ 
ponents of expression (7.27) are arranged in the canonical 
form following the prescription (7.124) with the result 




0 

0 

0 

0 

0 

0 




0 

fe33,ll 

i^^23,ll 

0 

i^23,10 

^^33,10 


^AB ■ 

1 

0 

0 

i^23,ll 

0 

fe22,ll 

0 

0 

0 

fe22,10 

0 

fe23.10 

0 

(7.142) 



0 

te23,10 

fe22,10 

0 

^22,00 

'1^23,00 




0 

i^^33,10 

^■^23,10 

0 

^23,00 

i^^33,00 


^23,10 

and 

fe.oo may be eliminated by suitably orienting the 


reference tetrad. A study of (7.142) shows that it is of type 
II with a = — 0 . It follows therefore that all physically 

interesting gravitational waves approach Petrov type II 
(null) with a = = 0 , at large distances. Since (7.142) is 

valid without approximation, we have shown rigorously 
that the Riemann tensor is Petrov type II (null) at all points 
where a gravitational wave is locally plane (31). 

Robinson and Trautman (30) have given some exact 
spherical solutions representing waves, with a Riemann 
tensor which again approaches Petrov type II (null) at large 
distances. Their metric is 

ds^ = 2dpd(x {K — 2Hp — 2mlp)do^ 

~ + \A'n 4- 

Here m is a function of o alone, p and q are functions of 
<7, ^ and 7], 

H = -I- p(p-^q),g„ — 

K is the Gaussian curvature (chapter 3) of the surface 
p = 1, <7 = constant. 

K = V\{lnp)^^g -f (^np),^J 
For this metric, == 0 reduces to 

+ ^,vv ~ “1" ^~ 4p“^(m 0- — BRm) 
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The Riemann tensor may be written 

^[ivafi P ^ ^(ivafi H~ P P fivccfi 

D, III, N are tensors of type I (degenerate), type III, and 
type II (null), respectively. They are covariantly constant 
on any ray of constant cr, rj 

The solutions are degenerate type I if m 0 and K is 
independent of rj and are reducible to m = 1, p = cosh 
g = 0. p is a real or purely imaginary constant. If is real 
and 0, the Schwarzschild solution results. 

If {K gY 4 - ^ — 0, the solutions 

correspond to type II non null or type III, with type III 
resulting when m — 0. If m = 0 and K is independent of 
1^, 7] the solutions are type II null or flat, with 



the condition for flatness. 

Solutions periodic in o can be constructed. For source 
velocity <c c, there is at least one singularity on any wave- 
front where fields do not vanish. 
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CHAPTER 8 


Detection and Generation of Gravitational Waves* 


8.1 Detection 

Suppose that we have a system of masses which may 
interact with each other. We start with the action principle 


^/ = 0 

The action function I may be written 

(8.1) 

. 

J = —me 

1 % 

ds + W 

(8.2) 


Here m is the rest mass and W is the part of the action func¬ 
tion associated with nongravitational forces arising from the 
motion of the mass relative to other masses with which it 
interacts. The principle of virtual work suggests that for 
dW vfe assume a function given by 

dW = - ^ {Fds (8.8) 

(8.3) identifies as the four force. The Euler-Lagrange 
equations resulting from (8.1) are arranged using the method 
of Section 3.6 to obtain 


ds^ 


dx°- dx^ FP' 
d$ ~ mc^ 


(8.4) 


t The results of this section appeared in the author’s Gravity Research 
Foundation prize essays, April, 1958, and April, 1959. This chapter follows 
the paper Phys. Rev. 117, 306, (1960). We thank the Physical Review 
for permitting its publication here. 


124 


GENERATION OF GRAVITATIONAL WAVES 


125 


(8.4) may be expressed in terms of the four velocity 
27/* = dxf^jds as 


d 

Js 




dU/^ Ff^ 
Ss mc^ 


(8.5) 


The symbol djds means the covariant derivative with 
respect to s. 

Following essentially the methods of Synge and Schild, 
we introduce a parameter v such that the world line of each 
mass element corresponds to a given value of v. Taking the 
CO variant derivative of (8.5) with respect to v gives 


< 3 ^ 27 /* dFf^jm 


dvds 

Direct calculation leads to 

^ 227 /* ( 3227 /* 


c^dv 


( 8 . 6 ) 


<5^5^; dvds 




dxy 

dv 


Employing this commutation law for covariant differ¬ 
entiation enables us to express (8.6) in the form 




dvds 


c^dv 


(8.7) 


In (8.7) dxyjdv is a unit vector normal to the world lines, and 
the four velocity 27/* is a unit vector tangent to the world lines. 
The vector nr, defined by 


nr 


dxy 

dv 


dv 


( 8 . 8 ) 


is an infinitesimal vector joining points with the same value 
of s on neighboring world lines with values of v differing by 
dv. The covariant derivative of dxyjdv with respect to s can 
be written in the forms 


d /dxy\ d (dxy 




- 




ds \dv / <3u \ ds 


(»r\ 

w) 


dUy 

dv 


(8.9) 
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Employing (8.6), (8.7), (8.8), and (8.9) then gives 

+ = ( 8 . 10 ) 

If the particles do not interact, the right side of (8.10) 
is zero and 

4- R^ac^yU-n^Ur = 0 (8.10a) 

For (8.10a) the particles move along geodesics. (8.10a) 
is called the equation of geodesic deviation. 

8.2 Mass Ouadrupole Detector 

In order to discuss the detector* of Fig. 8.1, we imagine 



Figure 8.1 


the two world lines are those of the two masses, 
given by 

ny = ry ^y 

with ry defined by 

dry 

= 0, (for all s) 
ry -> ny 


Let ny be 
( 8 . 11 ) 


( 8 . 12 ) 


in the limit of large internal damping and all components 


t An arrangement somewhat similar to this was independently sug¬ 
gested by H. Bondi at the Royaumont conference in June, 1959. 
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of Rf^oifid = (8.10) becomes 


IL 

mc^ 


(8.13) 


In (8.13) we have denoted by //* the differences in (non- 
gravitational) forces at the two masses. For />“ we assume 
a restoring force —and a damping force — cD^g^{d$^lds). 
hi^a, and are tensors associated with the spring. (8.13) 
then becomes 


^ I 

ds^ cm ds mc^ 


U-Uy[r^ + 


(8.14) 


We now let time run in the direction of the tangent to 
the world line of the center of mass. The oscillator is in free 
fall. We use coordinates in which the Christoffel symbols 
vanish and write (8.14) in the approximate form^ (assuming 
i < r), 


dt^ m dt m 




(8.15) 


In (8.15) we see that the driving force for the harmonic 
oscillator is the Riemann tensor. Measurement of displace¬ 
ment amplitude or power absorbed enables one to calculate 
certain components of the Riemann tensor.** 


t If the particles arc free, the second and third terms of (8.15) vanish. 
In Section 7.2 wc showed that a locally plane gravitational wave traveling in 
the direction has to ssero. It therefore follows from (8.15) that 

if the vector whi<!h joins two closely spaced free particles is in the direction 
of propagation of the wave, no relative disi)lacements occur. For free par¬ 
ticles the maximum relative displacements occur if tlie particles are in the 
transverse plane. This result has appeared in the work of Firani. 

Measurement of the liiemann tensor by comparing accelerations of 
free test particles has been considered l)y F. A. Firani.2.3 The results of 
this chapter indicate that interacting partiedes must be used, in practice. 
The correspondence V)etween voltage in a piezoelectric crystal and 
which is discussed in Section 8-2, may provide a basis for consideration 
of measurement problems in quantized general relativity. A very small 
volume of (trystal is needed in consequence of the fact that the wavelength 
of the acoustic resonance vibrations is five orders smaller than the wave¬ 
length of the gravitational waves. 
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From (8.15) we may conclude that energy may always 
be absorbed from any gravitational field for which the Fou¬ 
rier transform of RiOio does not vanish in a locally Lorentz 
frame which is propagated along a geodesic line. This geo¬ 
desic is then the world line of the center of mass of our har¬ 
monic oscillator. 

Suppose now that sinusoidal (weak-field-approximation) 
gravitational waves are incident. An orthogonal comoving 
coordinate system is employed, with the oscillator oriented 
in the direction of the axis, and !)>“« are imagined to 
have one component only, k^^ — k and ~ Taking the 
Fourier transform of (8.15) leads to 


~ {co^m - icoDdi‘i — kdi-j,) 


(8.16) 


(8.16) is a maximum at resonance, —co^m -f- = 0. 

The total dissipation D == + Din where is the ex¬ 

ternal dissipation and Din is the internal dissipation associat¬ 
ed with irreversible processes within the antenna. The power 
which can be delivered to auxiliary apparatus with Dgx is 




m2c^(D/^0a0^“)^^ex 


(8.17) 


2(Dgx “h Din)^ 

(8.17) is a maximum when Dg^ = Din, s-nd the maxi¬ 
mum power Pm is given by 



8Din 


(8,18) 


The sinusoidal gravitational waves are now assumed to 
be radiated by a linear mass quadrupole oscillator. The 
transformation laws indicate that to a good approximation, 
R^oocoi as seen in a frame fixed in the center of mass of the 
radiator, is the same as that seen in a frame fixed in the 
center of mass of the detector, for small velocities. Using the 
solution for the linear mass quadrupole oscillator of Section 
7.5, the mean-squared value of Rf^Qao'^"' is calculated and aver¬ 
aged over all possible orientations of the receiver. Let 4r 
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the radiated power per unit area averaged over a sphere, for 
the linear mass quadrupole oscillator. The total radiated 
power P is given by 


P ~ 


\57tc^ 


(8.19) 


where is the amplitude of the quadrupole moment. (8.19) 
and the known expressions for the fields then give, for the 
mean-squared value of in a direction normal to the 

axis of the radiator 




\r\^G 


( 8 . 20 ) 


In (8.20) p is the propagation vector of the gravitational 
wave. Employing (8.18) and (8.20) gives 


|^|2 Q 
2cDi„ 


( 8 . 21 ) 


The influence of the internal dissipation will now be con¬ 
sidered. First we assume that no irreversible processes take 
place within the antenna itself and that is due entirely to 
radiation damping of the detector. The known solution for 
a linear mass quadrupole oscillator enables us to calculate the 
radiation resistance of the detector An as 


\r\^ 

__ 


( 8 . 22 ) 


(8.21) and (8.22) give 


15P 


M (radiatioa damping only) 


16 jr 


‘'or 


(8.23) 


The implication of (8.23) is that the average absorption cross 
section for a detector which is damped only by its own 
reradiation is 


15_^ 

1671 


(8.24) 
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We see from (8.24) that under these conditions the 
average absorption cross section is roughly a wavelength 
squared, and is independent of the constant of gravitation. 
Unfortunately the condition that the internal damping be 
only due to radiation cannot be attained in practice because 
other irreversible phenomena within the antenna are many 
orders greater than the radiation damping. In order to make 
this clear, we calculate the quality factor, denoted by the 
symbol (J, which is defined by 


CO (maximum stored energy) 

^ power dissipated 

The Q associated with radiation damping, denoted by Qr, is 



15c^ 

2Gco^m \r ^ 


(8.25) 


For an antenna at co = 27r X 10"^, a reasonable value of 
mr^ = 10 g cm^ and (25) gives Qr ^ 10®^. A practical an¬ 
tenna might be expected to have Si Q 10®. 

We therefore must deal with systems limited by internal 
damping orders larger than gravitational radiation damping, 
and under these conditions the average absorbed power will 
not be independent of the kind of antenna. For an antenna 
orientation arranged for maximum response, 


[{R'^occor^n 


|r|^ 


(8.26) 


(8.26) and (8.18) lead to power absorbed P^, given by 

^ 157tGm^P^\r\^ ^ 157tGmQin^^\T\^ ^ 

ScDir. 


8(oc 


(8.27) 


In (8.27) (Jin is the Q associated with internal irrever¬ 
sible processes, (Jjn = The cross section 5, implied 

by (8.27), is 

5 _ \r\^ 

Scoc 


(8.28) 
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For a continuous spectrum the absorbed power is 

7i^Gmc~'^P^\r\^toj.{coQ) (8.29) 


In (8.29) 4r(<^o) is the power spectrum of 4^- in the vicinity 
of the resonant frequency cdq. 

In order to further discuss these results we must consider 
the excitation of a continuous medium by a gravitational 
wave. This is necessary in order to be able to account for the 
interaction of the mass of the spring with the wave and to 
account for the effects of the finite velocity of propagation of 
the elastic forces of the spring. 


8.3 Interaction of a Crystal with a Gravitational Wave 


The starting point for our discussion is expression (8.10). 
The infinitesimal vector is from a reference point in the 
crystal to a neighboring point. The mass m is imagined to 
belong to an infinitesimal volume surrounding the neighbor¬ 
ing point. On the right side of (8.10) we must now include 
both elastic forces and dissipative forces. We write as 

(8.30) 

r/* is defined by the conditions 


— = 0, (for all s) 
r<“ —> 


(8.31) 


in the limit of large internal damping and flat space. We 
may now express Eq. (8.10) in the form 

(520 ^0 d20 

(8.32) 

Here we are assuming an isotropic crystal. The second 
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term accounts for internal damping and the third term ac¬ 
counts for the elastic forces. B and are normalized to 
unit mass density. Again is a unit vector tangent to the 
world lines. Since r/* can he arbitrarily specified, it follows 
that 


s^o so s^o 

(8.83) 

In (8.33) the fourth term is clearly symmetric in the 
indices v and fj.. The last term in (8.33) may ordinarily be 
dropped because it is many orders smaller than the fourth 
one. For the strain tensor which is the symmetric part of 
we may therefore write 


6^ 


(IV (sym) 

'ds^ 


66 


6 ^ 6 , 


I ^ ^^(iv (sym) I Q ^/iv (sym) 

ds 6x^ dx^ 




(8.33a) 


We now consider a special case of (8.33a), namely ex¬ 
citation of longitudinal acoustic waves. An approximate 
form, for waves in the direction of an orthogonal coor¬ 
dinate system (with the time direction tangent to the world 
line of the observer) is 

d^e d^e ee 

^ d{x^Y‘ ~ (8.34) 


In (8.34) Pm is the density, y is an appropriate modulus, and 
b is a damping constant. We assume that its origin 

in incident sinusoidal gravitational waves so that 


(8.35) 

In (8.35) the index k runs from 1 to 3. Let v, be the 
sound velocity \^lylpM'l> be the wavelength of sound, k^ == 
27r/Ag, cf.—hl2pMVs, andThen to a good approxima¬ 
tion the solution of (8.34) is 

6 = [Ay cosh yx^ — /ca -2 j 


(8.36) 
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The displacement of points relative to the center of mass of 
the crystal, which we denote by may be calculated by 
integration of (8.36) with respect to as 


A sinh yx^ 


f 




(1 






(8.36a) 


Making use of the boundary condition that 6 vanishes at 
the ends leads to 


f cos P-J, 


a>‘^y{icd sin k^l + cos k^l) 


(8.37) 


In (8.37) I is half the length of the crystal. The first 
term of (8.36) gives the contribution of the acoustic waves* 
and the second term gives the strains which would be set up 
if there were no internal forces at all. (8.37) must be modified 
if the crystal is piezoelectric. 

(8.37) has maxima when k^l is an odd multiple of 7r/2; how¬ 
ever it is clear from the denominator that the largest maxi¬ 
mum is the first one for which the total length is half an 
acoustic wavelength. The system composed of the two 
masses and spring (Fig. 8.1) must be described by an equa¬ 
tion such as (8.34), when the spacing of the masses approaches 
half an acoustic wavelength. It is clear that the largest 
value we can expect from (8.28) will occur when r is half an 


t A study of (8.30a) indicates that the apparatus located on the earth’s 
surface will behave as though it were in free fall. For consider the effect 
of gravitational waves on the internal motions of the earth, applying (8.36a) 
to the normal modes of the earth itself. The first term on the right represents 
the acoustic waves, the remaining terms represent the motion if there were 
no interactions. The ratio of the terms is, using (8.37), 

^iCO^A A, Vs 

”~7™ hxl “ c(a0 

The ratio of the velocity of sound to the velocity of light is about 10“'^, 
much smaller tlian the product of length and absorption coefficient for 
sound a in the earth. This argument is only valid for gravitational wave¬ 
lengths comparable with or shorter than the earth’s radius. 
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acoustic wavelength in the spring. This is an important 
limitation because the velocity of acoustic waves is about 
five orders smaller than the velocity of light, so the cross 
sections implied by (8.28) are ten orders smaller than would 
be the case if the elastic forces of the spring were propagated 
with the velocity of light. Such a limitation could be over¬ 
come in a number of ways. One might employ restoring 
forces transmitted by electric and magnetic fields, with the 
velocity of light. The piezoelectric effect may be employed, 
in which case the polarization charges on the crystal faces 

give rise to some stress components which do not change 
sign every half acoustic wavelength. 

In a piezoelectric crystal a strain results in an electric 
polarization P^. given by 

Here is the piezoelectric stress tensor. The electric 

polarization gives rise to an electric field over the crystal. 
Its integrated value may give a terminal voltage large enough 
to be observed with a low-noise radio receiver. Measurement 
of this voltage measures components of the Riemann 

tensor if a crystal with suitable constants is employed. 

The system of stresses in the crystal is modified in a 
significant way if it is piezoelectric. Additional terms in¬ 
volving the piezoelectric constants need to be added to 
Eq. (8.33). We consider a very simple example. Suppose a 
single longitudinal mode is excited, with sound velocity in 
the direction. Let the thickness in the direction be 
small and assume that the crystal faces normal to the x^ 
direction are plated with a conductor. The piezoelectric 
relations (5) are 

-T = 0Yo 4- DH 14^711 E = DjK + Hd (8.38) 

In (8.38) T is the stress, K is the dielectric constant, B 
is the strain, is the elastic modulus, E is the electric field 
intensity, and D is the electric displacement. Both D and E 
are assumed to have components in the x^ direction only. 
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H is the piezoelectric constant relating open-circuit voltage 
to strain. A study of (8.38) and the equations of motion of 
mass elements of the crystal indicates that a wave equation 
similar to (8.34) results with 

y = Yo- mKl4>7t 


Since the crystal surface normal to the direction is 
plated with a conductor, dEjdx^ — 0. At the free ends of the 
crystal T = 0. If the crystal is coupled to an external im¬ 
pedance Z we may write 



Edx'^—(4i7t)~^Z 



Ddx^ dx^ 


These boundary conditions and the wave equation (8.34) 
then lead to the result 

Q = [A^y cosh yx^ — (8.36a) 

where y, p, and / are as defined earlier. This has the same 
form as (8.36), but now the constant A-^ is given in terms of 
the length in the x'^ direction and lengths and in direc¬ 
tions x^ and £C^, and the “clamped” capacitance C as 



-\-ia)CZ){Y^^— H^K|4>7t)cosP;Jl|2+iH^K^l^o)Z sin ^iZi/2 ~ 
L8^/~( ylj2){1-^icoCZ) ^ia)ZH^K% sinh yl^l2_ 

(8.39) 

and the voltage which, appears at the crystal terminals when 
coupled to an impedance Z is 

^ = (i) (r-^iSz) ~ 

The electrical network theorems now permit straightfor¬ 
ward ealculation of the power which can be delivered by the 
detector to a radio receiver. For a crystal with constants 
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similar to ordinary crystalline quartz on which sinusoidal . 
gravitational waves are incident, the power which can be 
transferred is roughly 

Pa ergsjsec (8.41) 

In (8.41) CO is again the angular frequency and is the 
incident gravitational power flow in ergs per square centi¬ 
meter per second. is the Q of the crystal and associated 
electric circuit. A cubic meter of crystal at co 10® gives a 
cross section for absorption-- lO-^^ cm^. This is for perfect 
crystals driving amplifiers without positive feedback. Im¬ 
perfect crystals may have cross sections three orders smaller. 
For a continuous spectrum with a power spectrum function 
^or(ft>)> the power absorbed is 

Pa ^ 10”^^ F4j.(coo) ergs/sec (8.42) 

(8.41) and (8.42) provide a basis for discussion of sensi¬ 
tivity. In microwave spectroscopy it has been found that all 
spurious effects other than random fluctuations can be re¬ 
cognized. A similar assumption will be made here. The ran¬ 
dom fluctuations are partly thermal in origin, partly the 
result of spontaneous emission processes. For synchronous 
detection of sinusoidal waves, the power output of the de¬ 
tector must exceed the noise power (6) Pj^^i given by 

p _ Nhco 

k IS Boltzmann’s constant, T is the gravitational antenna 
ernperature, N is the noise factor of the receiver (expected 
o e ess t^n 25 and more than 1), and is the averaging 
time. A different expression is required if radiation with a 
continuous spectrum is being studied. In this case the power 
delivered by the detector must exceed 


/ TZ^CO \ i 

NHco 


'^oilkT _ 
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Experiments are being planned to search for interstellar* 
gravitational radiation** using methods described here. For 
the first method the earth itself is the block of material 
constituting the antenna. The normal modes of the earth 
(about 1 cycle per hour) are excited by incident gravitational 
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Figure 8.2 


waves. This procedure is limited by the relatively low Q of 
the earth and the high noise background of its core. The 
apparatus of Fig. 8.2 is employed in the second method, in 
which the strains induced in the crysal are employed. Search 


i J. A. Wheeler has noted (onziifeme Conseil de I’Institut International de 
Physique Solvay, La Structure et revolution de Vuniverse, Editions Stoops, 
Brussels, 1958, p. 112) that the density of gravitational radiation could be as 
high as 10-29 ]()~28 g/cin3 lo® ergs/cm^ sec) and still be consistent with 

present information al)out the rate of expansion of the universe. He and 
M. Schwarzschild (private communication) have subsequently noted that if 
this radiation were set free by the same process which caused the inhomo¬ 
geneous collection of matter into galaxies, it would be characterized at that 
time, and therefore also now by the same scale of lengths, of the order of lO^^ 
cm today (10® years vibration period). 


®/(typical\ ^ PM ^ 

--- I -- 

dx / 6*2 


0.2 X 10-®« crn-2 


typical b.5 X 10~2® cin~^ X 10^^ cm 10“^ 

This would appear to be not too small, but too slow to measure, by these 
methods. 

tt Experimental work along tliis line is now being carried out by David 
Zipoy and Robert L. Forward, in collaboration with the author. The piezo¬ 
electric effect gives enhanced sensitivity when a mass which is many acoustic 
wavelengths on a side is used. At low fretiuencies this is not important be¬ 
cause a mass which is one-half acoustic wavelength long is already quite large 
and may not be obtainal)le as a single crystal. Excitation of resonant 
acoustic vibrations in a large block of metal, by the gravitational wave, is 
being consi<iered, along with the arrangements of Figs. 8.2 and 8.3. 
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at frequencies cycles/sec is planned. The earth rotates 

the apparatus. If radiation is incident from some given direc¬ 
tion it may be observed from the diurnal change in amplifier 
noise output. The arrangement of Fig. 8.3 should not require 



Figure 8.3 

rotation. If radiation is incident it will cause correlated out¬ 
puts. All sources of internal fluctuations will be uncorrelated. 
Low-noise amplifiers such as masers (7) may be employed. 

The apparatus must not respond to earth vibrations. 

Dr. Zipoy has proposed ingenious acoustic filters with am¬ 
plifiers to sense, then drive the mounts to cancel the effects. 

8.4 Rotations Induced by Gravitational Radiation 

Dirac (private communication) has suggested that as¬ 
tronomical anomalies might be correlated with effects of 
gravitational radiation. To discuss this and to consider 
detection by observing rotations we return to expression 
(8.5). Let a group of masses be situated near the space 
origin of a coordinate system and let the infinitesimal vector 
rf^ be the position vector of one of the masses. Let 
the Levi Civita tensor density. Multiply (8.5) by to 

obtain 

^ d d , .1 

(8,43) 
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Here again is a unit vector tangent to the world lines 
and in the second term of (8.43) we have used the identity 
~ the world line of the origin be a path for 

which the Christoffel symbols vanish. Then it follows that 
(8.43) can be written 

^ d dx°^ d ^dx^ 

SiK" UrU^r<^ (8.44) 

In these coordinates = dr^^ysjdx’^, so (8.44) becomes 

^ doc^ d dcc^ 

~ ^ ^jiapK R^ycjsUyU^r'^r^ (8.45) 

If we now use (8.43) and (8.45) and sum over all masses we 
obtain 


2 

masses 


d 

ds 


dx^^ 

rfs" 

masses 



masses 


mc^ 


(8.46) 


(8.46) is a generalization of the relation between torque and 
change of angular momentum. If there are no nongravita- 
tional forces acting and if we take the time direction tangent 
to the world line of the center of mass, (8.46) becomes 


2 ^ ^ficc^K J . “ 2 ^/4a/?/c (8.47) 

masses masses 


We may apply (8.47) to t‘he calculation of the irregular 
fluctuations in the period of rotation of the earth caused by 
incident gravitational radiation with a continuous spectrum. 
Under these conditions a straightforward calculation leads to 
the result 


(ZIL) 


2 

ms 


£2 


^57zG 


^or 


(8.48) 


Here [ALY^^ is the mean-square fluctuation in the 
earth’s angular momentum, L is the angular momentum of 
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rotation, and is the total gravitational wave flux in ergs 
per square centimeter per second, assuming its Fourier trans¬ 
form is concentrated near zero frequency. If we arbitrarily 
assume that all the earth’s rotational anomalies are due to 
incident gravitational waves, ior is calculated to be 5 X 10® 
ergs/cm^ sec. It is clear from this that the earth’s rotation is 
not a useful detector unless the size of the anomaly can be 
reduced. The other astronomical anomalies lead to larger 
figures. 


8.5 Generation of Gravitational Waves 


It would be very desirable to be able to generate gravi¬ 
tational waves with sufficient energy to be detected in the 
laboratory. A number of important experiments could be 
done. 

For a spinning rod, (7.56) gave the formula for the 
radiated power Pr as 

Pr = 1.73 X 10“®® CO® ergs/sec (8.49) 

For a given moment of inertia the angular frequency co 

can be increased until the rod ultimately breaks. If we write 
the maximum value of a> in terms of the tensile strength and 
express the result in terms of the elastic modulus and strain, 
we obtain for the length Z the formula 

I = ^V[2(5]/jr (8.50) 


In (8.50) d is the maximum allowed strain for the ma¬ 
terial and Ajj is the wavelength of sound in the rod at the 
angular frequency of rupture. The implication of (8.50) is 
that the wavelength of the gravitational waves which can 
be radiated by a rod is at least 1,000,000 times the length of 
the rod. Also the moment of inertia is limited to values less 
than about 




(8.31) 


In (8.51) Pm is again the density, and we are considering a 
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fairly slender rod, for which the length is an order larger than 
the lateral dimensions. Employing (8.51) in (8.49) leads to 

PrC 4>X <5® CO-4 (S 52) 


(8.52) shows that contrary to the appearance of (8.49) 
low-frequency operation with large rods gives more radiation 
than high-frequency spinning of small rods. About 
ergs/sec can be radiated by such a 1-m rod. 

A new method for generation of gravitational waves is 
suggested by the field equations. The source of the gravita¬ 
tional field is the stress energy tensor. Time-dependent 
stresses can be produced electrically in a piezoelectric crystal, 
and these give rise to radiation. The weak-field solutions are 
entirely adequate for this kind of discussion since the metric 
is expected to differ from the Lorentz metric by perhaps one 
part in 10^°. Recall the weak-field solutions 





IS r(3"/.'’')retatdeddV 

- 


(8.53) 


In order to apply (8.53) to the problem of radiation by 
a crystal we first assume that acoustic resonance is employed 
and that one-dimensional compressional waves are set up. 
The components of are then given by 


T 0 
0 


PmC^ 


V. 


V 


- cos oit cos 


[l~U{x^ — 1—Aq cos cot) 



— U(- 

1 

1 

-do 

cos cot)'] 

(8.54) 


pM ^ sin cot sin [1 • 

- U[x^ - 

1 - 

- Aq cos oj 

4) 


- U{- 

-x^ — 1 — 

d-o 

cos co^)] 

(8.55) 


— pM ^p'^8 (Ot cos kgX^ 

'[1 —U{x^ 

' — 

1 — Aq cos cot) 


- U{- 

-x^ — 1 — 

do 

cos coi)] 

(8.56) 


In these expressions it is assumed that the waves travel 
in the direction F,, is the particle velocity, Vg is the sound 
velocity, and C/ is a step function defined hy U{x) = 0 for 
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cc C. 0 and U(x) = 1 for a? > 0. is the vibration amplitude 
of the free end. Making use of these expressions in (8.46) and 
employing the Einstein form of the stress energy pseudoten¬ 
sor enables the total radiated power to be calculated. The 
maximum value of determined by the maximum strain 
which is allowed before rupture takes place. 

For quartz the result for acoustic resonance is 


P< 


HGpm^S^v, 




Xl0-« 


ifl' 


(8.57) 


X 10-12 


In (8.57) S is the cross-sectional area, the term with subscript 
CO gives the radiated power at the fundamental frequency, and 
the second term gives the power radiated at twice the funda¬ 
mental frequency. The resonator must be a multiple of a 
half acoustic wavelength long. The first term of (8.57) is 
seen to be independent of This is because for > 1 the 
crystal is essentially an assemblage of electric quadrupoles 
with a given quadrupole driven out of phase with respect to 
its nearest neighbors. The resulting fundamental frequency 
radiation is approximately that of a single quadrupole. Each 
half wave section has an equivalent moment MAQXj27t, 
where M is the mass of a single half wave resonator; Aq and 
are defined by (8.54), (8.55), and (8.56). If a large number 
of separate resonators are located within a region of linear 
dimensions less than a gravitational-wave half wavelength, 
the radiated power will be proportional to the square of the 
total number of crystals. In order to radiate ergs/sec 

at the fundamental frequency, 10® crystals would be needed, 
each one half an acoustic wavelength thick and with a cross- 
sectional area of 50 cm^. A complex phasing arrangement 
would be needed in order to properly drive the array. 

It appears better to suppress the acoustic resonance 
vibrations and create, by the piezoelectric effect or by electro- 
striction, mechanical stress components which do not reverse 



generation of gravitational waves 143 

sign every half acoustic wavelength. In order to see that this 
is possible it is only necessary to study the solution given (8 ) 
in the literature for the component T ^ of the stress in a longi¬ 
tudinally vibrating crystal with as the thickness di^^ctinn, 
with a conducting plating on the faces normal to x . The 
thickness is assumed small, for simplicity, and an externally 
applied electric field parallel to drives the crystal. The 
component Tn is then 


T 


11 




a)(L—-x) , . ojx 

sin -h sm 


sm 


coL' 


V, 


V, 


V, 


sm 


cojL 


V, 


sin coi (8.58) 


In (8.58) dij is the tensor relating stress to electric field, 
E is the electric field in the direction of x^, and L is the total 
length. At resonance, ico/w, = n, and losses would ha.ve to 
be taken into account in (8.68) by a hyperbolic function m 
the denominator. However, off resonance, for example at 
Lcojv, = 7tl2, (8.58) becomes 

sin (ot (8.59) 


T 


11 




cox . (OX 
cos-h sm ■— 


V 


(8.59) is seen to have a component — d^iE^ sin cot, which 
does not reverse sign every half acoustic wavelength. A 
single large crystal, driven in this manner, will then ^ve 
volume-integrated stress components which are large. The 
radiated power would be expected to be 

- -T20;^ ^ ^ ^ 

In (860) P„,a, is the effective tensile strength in dynes 
per square centimeter and A is again the gravitational-wave 
wavelength. The alternative method is to use electrostric- 
tion, either in a solid or a liquid. For a liquid a region coin- 
parable to a half gravitational wavelength on a side would 
have to be used. The space derivatives of the driving electric 
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field would give rise to electrostriction stresses. A result 
similar to (8.60) gives the radiated power. 

Waves 1 m long could be radiated by a crystal with di¬ 
mensions 50 cm on a side. If it is driven just below the 
breaking point, each cyrstal would radiate 10“^® ergs/sec, 
assuming Pniax is its static published value. This amounts to 
100,000 gravitons/sec. Single crystal detectors of the type 
considered earlier may detect a power of about 10~^ ergs/sec 
at these wavelengths. A large gap therefore still exists be¬ 
tween what can be generated and what can be detected in a 
small laboratory. Complex detection and generation arrays 
can narrow this gap. Large amounts of electrical power 
would have to be dissipated in crystals driven to the fracture 
point — perhaps 10® watts in a crystal 50 cm on a side. This 
might well be substantially reduced if low-temperature opera¬ 
tion can be achieved. Also one might hope that low-tem¬ 
perature high-frequency operation might raise the effective 
tensile strengths. All these issues need careful experimental 
investigation. If the numbers employed earlier cannot be 
improved upon, it would require a crystal roughly 100 m 
on a side, and a large detection system, to generate and de¬ 
tect the gravitational radiation. 


8.6 Other Radiation Experiments 

The foregoing discussion was limited to very low energy 
gravitons. A large number of calculations were carried out 
by the author for the purpose of studying the feasibility of 
radiation experiments at higher energies. The results in¬ 
dicate that the situation is probably hopeless at intermediate 
energies corresponding to atomic, molecular, and nuclear 
transitions. As we have seen, the equations of general rela¬ 
tivity have solutions in the weak field approximation which 
are very similar to those of electrodynamics. The lowest- 
order radiation processes in electrodynamics give radiation 
proportional to ^ where e is the electronic charge. An 
electronic transition for which the selection rules permit 
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graviton radiation would in lowest order be calculable by- 
expressions similar to those of electrodynamics except that 
for we need to substitute Gw}, where m is the mass of the 
electron. There are other, less important modifications. 
Another factor of the order of the square of the ratio of the 
Bohr radius to the wavelength would have to be included 
because quadrupole radiation is the lowest nonvanishing 
order. Now, Gm^je^ fn so it follows that graviton 

emission will be at least forty-three orders smaller than 
photon emission. At extremely high energies the situation 
will improve because the mass increases while the charge 
does not. For Gm^je^ 1, energies of the order of 10®® ev 
are needed. Since even 10®® ev is beyond the range of present¬ 
ly conceived machines the outlook for such experiments is 
not promising. 
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CHAPTER 9 


Selected Topics in General Relativity 


There are reasons for thinking that the elementary formations 
which go to make up the atom are held together hy gravitational 
forces. A. Einstein 


9.1 Unified Field Theories 

In this concluding chapter we give a brief account of 
some important aspects of general relativity which are very 
significant, from a philosophical point of view, and may lead 
to developments in elementary particle theory. 

The success in geometrization of the gravitational field 
led Einstein and others to search for a way to geometrize 
electromagnetism. The desirable elimination of forces by 
having particles move along geodesics ceases as soon as 
electrical forces are introduced. This is indicated by Eq. 
(4.54), which has the Lorentz force included in it. 

When the general theory of relativity was formulated it 
was believed that the only forces in nature were those of 
gravity and electromagnetism, and that by geometrization of 
electromagnetism a unified treatment covering all of physics 
would result. The field equations relate a geometrical object 
on the left to a physical object on the right. Einstein felt 
that this was undesirable and that the use of the electromag¬ 
netic stress tensor in this way was at best a provisional solu¬ 
tion. Also Einstein never accepted the postulates of the quan¬ 
tum mechanics. He hoped that a completely deterministic 
description of elementary particle and quantum phenomena 
would emerge from a search for the unified field (4, 5). 

No complete geometrization of electrodynamics was 
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achieved. Indeed the search for a unified field theory would 
seem now to be a much more formidable task since nuclear 
forces and the phenomena of high-energy physics would have 
to be included. 

It is interesting that the 1916 general theory of relativity 
and classical electromagnetism can be expressed in terms of 
the geometrical quantity alone, without ordinarily alter¬ 
ing the physical content. This partial geometrization t was 
accomplished by Rainich (1) in 1927 and has recently been 
enlarged by Wheeler (2, 3) and Misner. Suppose we have 
only gravitation and electromagnetism, without charges. 
The Maxwell stress tensor (4.15) satisfies the relations* 


t No geometrical interpretation is given to the Maxwell tensor. It is 
eliminated from the equations by expressing the Maxwell stress tensor in 

terms follow directly from (4.15). To deduce (9.3) and for 

proof of other equations in this section, the following relations are very 
helpful. A generalized kronecker delta may be defined in terms of the Levi- 
Civita tensor density by writing 

Saf^ 


eafiyK 


Sf < 57 " 

daP Sffp dyP d^P 


, ^5/^ , 

This indicates that the generalized kronecker delta changes sign on inter¬ 
change of any pair of upper or lower indices, and that it vanishes if the same 
index appears more than once in citlicr the uioper or lower groups. 

Most authors take = 1 and - 1- Care is required whea 

employing both Eafiyt and in the same ealculation. 11 wo take the 

of FIIP 

then 

A straightforward reduction shows that to employ and remain consist¬ 

ent with the above expressions we must write: 

These expressions may be employed to cstablisli the fohowing identity, 
valid for any two antisymmetric tensors in a four dimensional space: 

Bua — *Aaa 


(footnote confd) 
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== 0 (9.1) 

^00 > ( 9 - 2 ) 

= (9.3) 

If the gravitational field equations are written in the form 


R ^ 


, „ „ StzG ^ 


/* 


then contraction on the indices and the use of (9.1) leads to 


R = 0 

so the field equations are 

(9.4) 

„ StzG 

= - 34 - Lr 

(9.5) 

The use of (9.3) and (9.5) lead to 


R/R/ = IdJR^^R’^r 

(9.6) 


The geometrical relations (9.4) and (9.6) apply to any 

whose sources are a divergence-free Max¬ 
well tield. Another relation is needed to guarantee that the 
antisymmetric tensor out of which (and ie„,) are con- 

stmcted really satisfies Maxwell’s equations. This additional 
relation is 


refx^vR^'^->‘Ry‘' V—g] 

"^iXfiyR^'^’^Ry^ \/—g~ 

L R^rR'^^ J,,~ 

R^^ R<'^ 


J,/S 


(9.r) 


\ - A V • 1 ) ana (^y.7 ) are the formulation for gravity and 
electroma^etism, in the absence of charges. The expression 

(9.7) requires a somewhat lengthy proof which we outline, 
following Misner and Wheeler. 


str^s “"if *** electromagnetic field tensor, the Maxwell 

stress tensor may be written 





Multiplying this by T**' as 


given by (4.15) then gives (9.3). 
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We shall require the dual of the electromagnetic field 
tensor defined here by 

(9.8) 

An operation e** is defined by 

= F^^ cos a + *Fj,y sin a (9.9) 

We say that (9.9) constitutes a duality rotation on F^y,^ The 
tensor is defined by 

(9.10) 

corresponding to a duality rotation —a. 

By an appropriate choice of a we may obtain a tensor 
which is more simple then F^^. Consider the invariants 

= iiPccfi COS a — *F^^ sin a)^ ^ 

== iFa /5 F°^^ cos 2a sin 2a 

sin 2a -1- cos 2a (9.11b) 

We recall that (9,11a) is — E^, and (9.11b) is 2E • H in a 
Lorentz frame, with E and H the electric and magnetic fields 
associated with 1^^. Select the angle a such that (9.11b) 
vanishes; the other invariant then becomes 

(9.12) 

The angle a is further restricted as to quadrant by requiring 
the sign in the right expressions of (9.12) to be negative. A 
field for which these requirements have been met is caUed 
an extremal field; it is a pure electric field, in the appropriate 

t The dual of the dual of a tensor is the negative of the original tensor, 
in consequence of the fact that g is negative. The dual operation corresponds 
to a duality rotation angle of ?r/2. It is irni)ortant to note that 

Also, from the identity of the previous footnote, 
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Lorentz frame. The extremal field is a convenient quantity 
to manipulate. For the remainder of this section will 
denote the extremal field. An arbitrary field may be ob¬ 
tained from by a duality rotation and a scale factor. It 
may be verified, using the formula ( 4 . 15 ) for the Maxwell 
stress tensor that all the components of are unchanged by 
® rotation. Given an which originates in a Max¬ 

well field, we cannot uniquely specify the Maxwell field, but 
we can determine it up to a constant times a duality rotation. 
It will be convenient to write the Maxwell equations in terms 

Of then to write another expression quadratic in ^ in 
terms of R which then ler ^ ' " ' — “ “ ^ 

tions are expressed as 


0 = 




These are rearranged to give 




s to (9.7). The Maxwell equa- 

cos a -f- ~ 

- — 1 sin a 

dx^/ 

\ ^ 

(9.13a) 

jsinoc+(*^/*p,- 

dcx.\ 

~ — I cos oc 

CX^J 


(9.13b) 

dcL 

dx- ^ 

(9.14a) 

doL 

dx- ^ 

(9.14b) 


n/r 1 • 1 (9.15) 

Soc -2(*f^ J 


dcc^ 


in teams' Tthi'RiccTtetor rthat 
terms of geometrical objects? ^ ^ ^ 
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For convenience a fourth-rank tensor is introduced 
which is constructed out of the Ricci tensor so as to have the 
symmetry of the Riemann tensor. Such a tensor is 
given by 

+ SR,- - {9.17Si) 

and 

- Rr^>^) (9.17b) 

By considering the extremal Maxwell tensor in a Lorentz 
frame it may be verified that 

g- = — iocfiiyS - (9.18a) 

This is a tensor equation; therefore, it is valid in general. 
It may also be shown that 

~ (9.18b) 

By solving (9.18a) and (9.18b) for we obtain 

0 

E^„ys £<rr’'* (9-19) 

Now define the tensor Fafiys 

EaftyS ~ ‘K ^yS fLvEa^pf^^ == ^( ^ySRfi^ ^P^Ra^^ 

(9.20) 

and calculate as 

= ¥ySnv{^a^Rp^ ” R^'^) R^^’y ^/-g (9.21) 

== is,sy,R^^->^R/V~~g 

(9.5) and the relation for in terms of led to the rela¬ 
tion (9.12); — (G/c^) — (R^^Rj^^)K Comparing (9.21) 

and (9.16) we see that (9.16), which is obtained from Max- 
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well’s equations, can be written in terms of the Ricci tensor 
as 


dec v~~g 

dx^ R^r 


(9.22) 


Since (9.22) is the gradient of a scalar, its curl vanishes, so 
we obtain expression (9.7), 

Given the Ricci tensor we may determine the Maxwell 
field (up to a duality rotation) in the following way. Use 
(9.22). Select some reference point and calculate 

rx 

ai{x) = J^a + aQ 

Since the curl of a ^ vanishes, this is independent of the path, 
a^j is arbitrary. This value of a may be used to obtain 
Fnv = ^iivi since (9.19) may be used for the purpose of 
calculating at each point. (9.22) does not exist for a null 
field. Another fundamental difficulty is that allowed geomet¬ 
ric data on an initial hypersurface may correspond to more 
than one Maxwell tensor (45). This non-uniqueness exists, 
for example, when the initial Maxwell fields are localized in 
two separated regions. Duality rotation in one of the regions 
gives different physics but may not alter the geometric data 
for a finite time. 



Wheeler notes that charges may be described by 
admitting regions with non-Euclidean topologies which 
are multiply connected. A schematic diagram of a 
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pair of charges formed in this way is shown in 'Fig. 9.1. Such 
a theory is equivalent to a theory with charge and a Eu¬ 
clidean topology. Misner has noted that the total number 
of lines which thread such a “handle” is a constant of the 
motion. 

Wheeler has also discussed stable solutions of the com¬ 
bined Maxwell-Einstein equations. His“geons” (gravitation¬ 
al electromagnetic entity) are objects endowed with mass 
which are made up of fields, held together by gravitation 
(34). 


9.2 Equations of Motion 


It was noted by Einstein and Grommer (6) in 1927 that 
the equations of motion of a system of masses are contained 
within the gravitational field equations and therefore do not 
need to be postulated separately, as in electrodynamics. 
Bergmann has remarked that constraints on the possible 
motions appear as a result of relations such as the Bianchi 
identities; the nonlinearity then leads to laws of motion such 
as Newton’s laws. An extensive literature exists on the de¬ 
duction of the equations of motion and their connection with 
the problem of radiation and the conservation laws. We 
should like to stress that the proposed methods for detection 
and generation of gravitational radiation, given earlier, in¬ 
volve electrical as well as gravitational forces. Any conclu¬ 
sions derived from consideration of the equations of motion 
and the radiation problem for gravitation alone may not 
apply to the results of Chapter 8. 

The following discussion of the equations of motion is 
due to V. Fock (T). In the zeroth approximation, space is 
considered to be flat, and a Lorentz metric is used. The stress 
energy tensor is then taken to be (Latin, space indices) 

TOO ^ to* ^ pmCv^ (9.23) 


In this approximation, TO^.^ 

ai'oo ' ^xoi 
H" 




and 


cdt 


dx^ 


0 


(9.24) 
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The weak-field solutions (7.13) are 


T’oo 


4G r (r'oo)retarded^^^^^ 


J 


|r — r 1 

2G 




4Gr (r-o,),,tarded^V 
J |r — r'l 


<Pok 


r iT'oo)retarded^''^'^dV 


(9.25) 


dt^j 


For low velocities z;/c <C 1, and in this approximation 
only <Pqq and <??ot need be retained. Using (7.2) and (7.10) it 
then follows that 


^00 ^ 1 + ^9^00? ^ (9.26) 

§0i ^ ^oi (9.27) 

In the next approximation we must consider = 0; this 

gives, using (3.38), 

~ + T- + T>“‘ = 0 (9.28) 

which we write as 
^T'oo 

TIT + 








+ r\, T- + r^.v == 0 


(9.30) 


The last expression, (9.30), will now be seen to be the 
equations of motion in the Newtonian approximation. Again, 
at low velocities, Using this and retaining only 

first-order terms then leads to 


a 

yj(pM’y’)+-g^ 


GpM 

c2 



(^'oo) 


retarded 


d^x' 



= 0 


(9.31) 


pM is the matter density when measured in the local rest 
frame of the matter. The first term on the left is the inertial 
term, the second is a force per unit volume due to the space 
derivatives of the pressure, and the last term is the gradient 
of the gravitational potential. 
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We may deduce the geodesic equation for a particle, 
following Papapetrou (35). First integrate the relation 
= 0 over the entire three space, employing the form 
(3.76). Surface integrals may be dropped because the particle 
is localized, and we obtain 

Now integrate the product xrT= 0 over the three space, 
again dropping surface integrals, with the result 

VtC J V 


^oj = - J (9-28A) 


(9.28B) 

Since these integrands are zero except for a small region 
it is appropriate to represent t>y series 

The prefix zero refers to the value at the point which is 
taken as the position of the particle. Account must now be 
taken of the structure of the particle. Papapetrou defines a 
single pole particle as one which has only some of the inte¬ 
grals J Tf^^(-~g)^(Px not equal to zero. A pole-dipole particle 
has at least some of the integrals f Tf<'''dau>{~~-g)H^x and 
c Xf^^'{—gfcPx not equal to zero. Extension to higlier order 
multipoles is evident. For a pole particle only the first term 
in thl series for is needed and the relations (9.28A), 

(9.28B) give 

r , , dx'^ 

J T«i^{—g)‘<Px = ^^0 

In (9.28D) lot fi — 0, thin givoK 

U‘‘ = imx^ldxo, 




(9.28C) 

f7V'»(- 

—g)^(Px 

(9.281)) 

W 

in tern 

ns of tlie four 

velocity 
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The rest mass ttiq has been set equal to {U^)~^ f T^^{—g)^d^x, 
Substituting (9.28E) back into (9.28D) leads to 

J = moU^U^/Uo (9.28F) 

We may now employ (9.28F) in (9.28C). The result 
after using (9.28E) is 

dUi^ 

= 0 (9.31A) 


which is the geodesic equation. 

Extension of this procedure gives equations of motion 
for a pole-dipole particle. The spin is defined by 

= jda:'^TfO{—g)U^x — J dxfT‘‘’>{—g)*d^x 


and the equations of motion for are 


ds 


+ 


31 

ds 


- UfiUp 




= 0 


(9.31B) 


Again, as in chapter 8, dfds denotes covariant differentiation 
with respect to s. 

Schiff (36) has employed (9.31B) to study the precession 
of a gyroscope which has a non gravitational constraining 
force F, in the earth’s gravitational field. His result to lowest 
brder is that the spin vector I® measured by a co-moving 
observer satisfies the equation 


dt 



X v/2mc^ + {BGMI2c^r^)r 


X v 


-1- (G//cV^)(3[cc) • r]r/r^ — co) 


X I® 


(9.31C) 


In (9.31C) V is the velocity of the gyroscope, J is the moment 
of inertia of the earth, o) is the angular velocity of the earth, 
r is the radial coordinate. The first term on the right gives 
the special relativity Thomas precession, the second and 
third terms are effects of general relativity. For a spinning 
top on earth, a precession of 10“® radians per day would 
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be expected. Observation of such effects on earth or in a 
satellite experiment would provide a test of general relativity. 
Somewhat similar suggestions were put forth (unpublished) 
by R. A. Ferrell (Gravity Research Foundation Essay, New 
Boston, New Hampshire, April 1, 1959) and later by R. A. 
Pugh (Weapons Systems Evaluation Group Memorandum 
11, November 12, 1959). 

It is possible to deduce the equations of motion of isolat¬ 
ed masses in a purely geometrical way (8, 9) using the left 
side of the field equations alone. Higher approximations to 
the equations of motion including effects of radiation have 
been obtained by a number of authors (10, 25, 27). 


9.3 Mach’s Principle 

The Mach principle states that inertia is due entirely to 
the mutual action of matter. Thus a particle is considered to 
require force to accelerate it relative to the other matter in 
the universe. At least part of the centrifugal force appears as 
a mutual effect, in general relativity. Consider the geo¬ 
desic equation for the space indices, keeping terms up to 
{vjcY: 


dV 


daP\^ dx^ dx^ 


Here v is the ordinary velocity of the particle. To first order 
in the Christoffel symbols are 

F^oo H" Ki,ii 

With these substitutions the geodesic equation becomes 


^dV Idx^Y' ["y?/ 
^ dx^ dx^ 


V(t%„„/2) 


d(eKA 

dt 


■ .. ditfi ^ 


(9.82) 
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We employ (9.32) to study the motion of a particle inside 
a hollow spherical shell of radius R and rest mass M, First 
let the shell rotate with angular velocity co with respect to an 
observer within the shell. Let Rg, 0, and (p be the spherical 
coordinates of an element of the shell and let a, Oq, and cpQ 
be the spherical coordinates of the point within the shell at 
which the field components are calculated. will denote 
the four velocity of an element of the shell as seen by an ob¬ 
server within. The rectangular components are 




coRs sin 6 sin (p 


^ (9.33) 

TJZ — A. r/O =- 

’ V[l-"Wsin2 0/c2] 

The mechanical stresses in the shell are assumed to be small, 
so the stress energy tensor is approximately pmC^U/^U^. The 
expressions (7.13) are then 


0 


V 


J |r — r'l 


(9.34) 


We are using 0^^^ for the quantity reserving cp 

for a coordinate. 


r—r'l = [^ 2 +— 2 ajR^[cos 6 cos 0 Q-f sin d sin Oq cos((p — (Po)J] ^ 

(9.35) 

r — r'l"^ may be expanded in the series 


r-r' -1 = 


R 


s 




a 


^Rs‘ 


_|_ (cos Q cos Oq -)- sin 0 sin 6q cos [(p —) 


Rs 




+ 2 ^ - 2 (sin 0 sin cos [cp— 930 ] + cos 0 cos 0 ^)^ - 1 - 


(9.36 ) 


Now pM is the density in the rest frame of an element of 
the shell, while d?x' is a three-dimensional volume element in 
the rest frame of the observer. In order to express the result 
in terms of the rest mass of the shell we note that a volume 
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element (Px" in the shell rest frame covers a region (Px' such 
that 

(Px' 

(Px" = (9.37) 

Utilizing (9.33) through (9.37) and carrying out the integra¬ 
tions then gives to order the values 


0 


00 


1 + 


3c^ 


2a>^a^ co^ sin^ 6 


15c2 




0 


11 


'o)^Rs^ ^ sin^ Oq sin^ cp^ oo^a^' 


Sc^ 


5c^ 


0 


22 


0 


33 


0 

^10 
0 


4>GM 
~^Rs 
4GM 
'^Rs 
4>GM 

0 

AiGMay^a'^ sin^ Oq cos cpQ sin cpQ 


0 


5c2 


m^Rs^ co^a^ sin^ Oq cos^ (Pq 
-1 


15c2 

oj^a^ 

'1^2 


“ 5c^Rs 

AGMcoa sin Oq sin <pQ 
s 


(9.38) 


4GMct)a sin Oq cos (^Jq 


20 


s 


0 . 


30 


0O>: 


23 


0 


0 




agm 


1 


<x> 




From these h^ii, may be calculated as 

— 0 fiv fiv0 (9.39) 

A study of (9.38) and (9.39) indicates that the last two 
terms on the riglit of (9.32) are of an order higher than {vjcy^ 
and the factor daPjds may 1)0 omitted from somci of the re¬ 
maining t(;rms. For this case the geodesic cciuation tlien 
takes the form 
d'^x 


/i2 .. 


(Is 


.2 


VK 


(d. ) a r + V X V X A 


(9.40) 
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where V = and If we now set = x, 

x^ = y, x^ = z and write (9.40) in component form, the 
equations of motion for a particle inside the shell are 


^d^x GM 
^ ds^~ 3c^ 
2 d^y GM 


5 


a>‘=‘x — So) 


^co^y + 8ct> 


dt_ 

dx~ 

dt 


(9.40a) 


2 d^z SGMco^z 

“ dP"" 15c^Rs 

It is clear from (9.40a) that we are observing centrifugal 
force on a body at rest inside the spinning shell. The effect 
is, in general, a small one — a factor ^ GMjc^Rs times the 
value resulting from rotation relative to the other mass in 
the universe. This result was given long ago by Thirring. 

If the shell does not spin but is given a translational 
acceleration in the x^ direction relative to a test particle in¬ 
side, we have, from (9.34) and (9.37), 


^10 — ^10 — cA-^ — 


AiGMvs 


(9.41) 


c^Rsil — Vs^jc^)^ 

Here Vs is the velocity of the shell. (9.40) then gives for a 
particle momentarily at rest inside 

^GM I v^v^ • 


dt^ d^R^iA — Vs^lc^)^ 


1 + 


i) 


(9.40b) 


C2(l — Vs^lc^), 

In (9.40b) is the acceleration of the shell. 

Consider a solid mass M instead of the shell and accel¬ 
erate M relative to a small mass m. If m is momentarily at 
rest its acceleration may be calculated using (9.32) and em¬ 
ploying the potentials of (7.13a). If we retain terms only up 
to vjc the acceleration of m is 


,.2 


ds^ 


GM 


R 


M 


(n- 


V • R^\ 4GM 


cR 


M 


+ 


c^R 


M 


r ^ cR^ 1 


a 


+ 


4GAr /a • R 


i ^ 

'UR^) 


(9.40c) 
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V and a are the velocity and acceleration, respectively, 
of M, and is the radius vector from m to M. (9.40c) con¬ 
tains terms which depend on the direction of the acceleration 
relative to the radius vector from m to M. These are roughly 
vjc times the isotropic effect of the acceleration. 

From a consideration of (9.40) we may conclude that the 
term VK includes the centrifugal force associated with rel¬ 
ative rotation, and the term {daPlds)^dAldt represents forces 
due to relative velocity and relative acceleration. The term 
V X V X A leads to Coriolis acceleration. For the rotating 
shell V X A gives an internal field similar to the magnetic 
field inside a spinning charged sphere. These considerations 
indicate that for the universe as a whole GMIc^R 1 if 
Mach’s principle is a valid one. 

This treatment shows that inertial effects depend at least 
in some degree on the mutual action of matter, according to 
general relativity. This does not satisfy the requirements 
of Mach’s principle because with no matter present other than 
a test particle, inertial effects still appear. Special relativity 
is a limiting class of solutions, if no boundary conditions are 
imposed. Wc may always choose coordinates such that in¬ 
ertia appears as a local property of space and all terms on the 
right side of (9.40) are transformed away. It is evident from 
an equation such as (8.13) that inertia will be unequivocally 
isotropic^ in consequence of the form of the stress tensor for 

t Cocconi and Salpeter (29, 30) have con.sidered the effect of anisotropy of 
inertia on the structure of s|)eetral lines. For the gamma radiation of nuclei 
in crystals, the atomic magnetic field leads to (2,7 + 1) equally spaced com¬ 
ponents of tiic nuclear energy level. With the anisotropy of inertia AM, each 
component is shifted by an amount {AMjM)TF^, where T is the average 
kinetic energy of the nucleon and is a coefficient whose value depends on 
.7, on the magnetic cpiantum number, and on the orientation of the magnetic 
field relative to the direction toward the galactic center. In the Mossbauer 
effect observations, the transition fretiuencies in absorber and emitter are 
compared. Relative motion is produced and the counting rate is observed 
as a function of the relative velocity of emitter and source. Resonance occurs 
successively as lines overlap, provided the selection rules permit excitation 
by the shifted lines. The number of peaks observed will then change if the 
relative spacing of the Zeeman components is altered and if the atomic fields 

(footnote confd) 
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a particle with scalar type rest mass. For these reasons the 
experiment of Hughes, Robinson and Beltran-Lopez (39) 
must be regarded as providing strong support for the present 
formulation of the general theory of relativity, and no sup¬ 
port for a strong form of Mach’s principle. 

Einstein (11) and more recently Wheeler have explored 
the possibility that Mach’s principle is not a consequence of 
the field equations, but is a boundary condition. If we have 
an isolated system we should require that at large distances 
the metric join in the proper way to the metric of the remain¬ 
der of the universe. Satisfying this boundary condition 
would represent the mutual action of the matter in the uni¬ 
verse on the mass in the “isolated” system. Carrying through 
such a program might lead to relations between inertial 
properties and the far-distant nature of the universe. 

Other points of view on Mach’s principle have been 
presented by Dicke and by Sciama (37). 


9.4 Remarks on Cosmology 

The equations of general relativity have solutions which 
appear capable of giving a description of the universe and its 
evolution. Long ago Einstein considered the problem and 
concluded that the equations do not admit a static solution 
or t e universe unless an additional “cosmological” term 
IS introduced into the field equations. Later it was shown 
by Iriedman that nonstatic solutions exist, without addi- 


are aligned. If the atomic magnetic fields are randomly oriented in emitter 

r T of inertia broadens the lines. A Mossbauer 

eriect experiment gave a negative result (38). 

.. elegant and much more sensitive test of anisotropy of inertia was 

eoently ^med out by Hughes, Robinson, and Beltran-Lopez (39). They 

solutiln'^of'r magnetic resonance frequency of the Li' nucleus in a 

the or°L2t^n ofTr ® rotation changes 

usuallv “PParatus relative to our galaxy. This method is un- 

Zton irthe r , 7 ““Sy ‘‘rises from a P, 

l^te vaL of V accuracy with which the abso? 

Their ressonance ftequency can be measured, 

xneir negative result establishes AMjM < 
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tion of a cosmological term. The discovery of a red shift of 
spectral lines, increasing with distance, suggests that the 
universe is indeed nonstatic, in an expansion phase at the 
present time. 

A universe with positive curvature would be closed, 
finite, and unbounded. This is a very attractive idea. Other 
cosmologies have been constructed, notably by Bondi, Gold, 
and Hoyle. We limit our discussion to Friedman’s solution 

The assumption is made that the universe is spatially 
isotropic. This is supported by the observation that the 
density of stars appears to be the same in all dhections. An 
appropriate metric is 

— ds^ = dl^ ~~ (9.42) 

dl^ — g^^dx^dx^ — ■+ 4~ dx^^) (9.43) 


In consequence of symmetry we can expect to have 
the same kind of space dependence at all times, so g^^ should 
be the product of a function ^ of a?® alone and a function a of 
-|- x'^^ + x^^ — T^. Here the coordinates x^, x^, x^, and 
are dimensionless numbers. We need expressions for the g^^ 
in terms of the curvature invariant of the space. Suppose 
we have a curved space and we introduce at a given point a 


pair of vectors and B^; the quantity R^/fys A‘=^A'y is 
a scalar. Another scalar is (gayg/Ji — aiid 

the tensor gay gfin — gasgfiy been constructed to have the 
same symmetry as Rafiys ■ The quantity , given by 


is also a scalar. Straightforward calculation shows that Kji is 
unchanged if A and B arc replaced by any linear combination 
of A and B. is called the Ricmannian curvature. For an 
isotropic space, will also be independent of the directions 
of Af^ and B^^, so at a given point Kji has the same value for 
arbitrary and B^; therefore the Riemann tensor in iso- 


(^oty ^5/? S 
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tropic space is given from (9.44) as 


SaS ^fiy 


) 


(9.45) 


(9.46) 


We employ (9.45) for the three space. The component 

is given by both (9.45) and (3.48). (9.45) gives zero for it 
and we write 

-^1213 = —^11,23/2 + ^^^[^33. 2 ^11,3] = 0 

W^e have made use of the fact that 

B Bt 3 

^11=^22 = ^33 Nc 
(9.43) and (9.46) give 

^ ^ da d^a 

dr) 


2(^1 + 


a 


dx^ dr 


0 


(9.47) 


r dr dr^ 

Expressions for components and may also 

be obtained by comparison of (8.48) and (9.45). These are 

^mbined and the term d^ajdr^ is eliminated using (9.47). 
The result is s v / 




(9.48) 


\dr/ ar dr 

(9.47) and (9.48) express the requirements on a and ^ which 
are imposed by the isotropy of the space. 

Solutions of these equations are of the form 


a 


^2 + (9.49) 

An adequate description will result if the constants are 
chosen so that 

1 4a 


K, 


4CoC 


a = 


Kr, = 


1 H“ ar^ ’ ^ ^2 

The metric therefore assumes the form 
^2 


(9.50) 


-ds^ = 


[dx^^ + _ ^ 0 - 


(9.51) 


[1 + 

^2 

= [1 4- a^2]2 sin2 6 d<f^ + rHe^} — dcd>‘' (9.52) 
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we again note that ^ has the dimensions of length, r is 
dimensionless. The “radius” is given by with 


For a > 0 




dr 


(9.53) 


2 Va 


(9.54) 


From (9.54) we see that a is a factor related to our unit of 
length. Coordinates may be used in which oc = 1, when 
(9.54) applies. The use of (9.45) to calculate the curvature 
scalar gives 

3/?,.^== (9.55) 

== 24ar^-2 j 

We see from (9.54) and (9.56) that for a = 1 we have a 
finite, unbounded universe with positive curvature. With 
(X = 0 we have a “Euclidean” open space and with a nega¬ 
tive we have a curved open space with negative curvature. 
The assumption is now made that the stress energy tensor 
can be assumed, to a good approximation, to have only the 
one component all others vanishing. The field 

equations are then 


R 


00 


itoOO 4" 




^nCrpM 


«2 


R 


11 


1 r. I 
kgxxR = (4a + 


\daf>l 




0 


(9.57) 


(9.58) 


with identical expressions for i ?22 “ 

From subtraction of these equations there results 


A<n^PMG 

+ . 


0 


(9.59) 


Since both p^ and ^ are positive, (9.59) indicates that 
cannot have minima, cannot have inflection points, and can- 
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not be time independent. For ‘‘Euclidean” space, 
(9.57) gives 

\ ^ dcc^) 3c^ 


a = 0; 
(9.60) 


For a 0 (9.58) and (9.59) give 


d^Y Sti^^pmG 
dx^f 3c^ 


(9.61) 


(9.58) has assumed the pressure to be zero. The total 
mass M at any time is given by 


M = 


^ \3 

\T~+~o^/ ^ 


(9.62) 


If the total mass is conserved, ^^pm is independent of 
time. We may therefore set StzG^^pmI^c^ = ^o» where ^^isa 



Figure 9.2 


positive constant, and obtain 



(9.63) 


For a = +1 we have a closed universe with the time de¬ 
pendence of ^ as shown in Fig. 9.2. If a = —1, ^ increases 
monotonically. 

The most attractive possibility, a = +1, can be tested 
by using Eqs. (9.60) and (9.61). 

1 d^ 

^ dx^ 


can be calculated from astronomical data on the expansion 
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rate, and pj\^ may also be estimated from observations. If 

/I cWY 

then a positive value for a is allowed. Present observational 
astronomy data ^ive too small a value of p]\^, to give positive oc. 
Wheeler and Ivanenko liave both remarked that the idea of 
a closed positively curved universe should not be discarded 
before we are certain that our knowledge of pm is sufficientlv 
good to rule it out. Other presently unknown Zuv^s ot 
energy may give values of matter energy density consider¬ 
ably in excess of what is indicated by observations usinsr 
visible light alone. ^ 


9.5 Hamiltonian Formulation (14-17, 19-22, 26) 

Fhe thcoiy wliich has been presented thus far is in 
Lagrangiau form. The functions (6.5) and (6.10) were shown 
to lead to the field equations, and (6.10) permitted us to 
write the Lagrange equations (6.14). 

In quantum theory the Hamiltonian form is the one 
which has been employed most successfully. In mechanics 
the Lagrangian L is a iunetion, of the coordinates and veloc¬ 
ities q^. 1 he Hamiltonian is a function of tlie coordinates 
and momenta and is defined as 


H = p.q. — L 

The momenta are given by 

Pi = 

The equations of motion ar(i 


dq^ 


(9.64) 


(9.65) 


Fur any otlu^’ dynamical variable, suy b, 

a?,. bp. i)i "" i)/f ] dp\dpi dq, ' 


(9.06) 


[b, li\ (9.67) 
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In (9.67) [6, H] is the Poisson bracket, defined by 

du dv du dv 

iu, v] — 

In order for this ‘‘canonical” formalism to be applicable 
the p’s and g’s must be independent variables in the following 
sense. No constraints, that is no a priori relations should 
exist to relate them.^ We have seen that several kinds of con¬ 
straints do exist in general relativity. Expression (7.96) ex¬ 
presses some of them as a set of equations which must be 
satisfied by the initial values of the fields, and also at all 
subsequent times. The term constraint is used also to express 
the fact that there is some arbitrariness in the theory — for 
example the gauge group in electrodynamics or some aspects 
of covariance. Here we are following very closely the dis¬ 
cussion of Dirac (16, 22). 

When the Hamiltonian form is considered, the con¬ 
straints may appear as soon as the momenta are defined by 
(9.65). We may find that the p’s are not independent func¬ 
tions of the ^’s. Suppose there are a number of independent 
relations which appear as a consequence of (9.65), which we 
write 

'P»(p. ?) = 0 (9.68) 

Dirac calls all equations of constraint and equations which 
are valid in consequence of constraints weak equations and 
writes them as q) 0. 

The method of “Lagrangian” multipliers will be employ- 

t Such a requirement corresponds, in the La^j^ran^^ian formalisni, to 
being able to solve for the Qi in terms of and q^. The Lagrainge ecpiations 

d 3L aL 
dt dqj dqj 

may be written in the form 

a^L aL 

The latter set of equations may be solved for q^ only if the determinant with 
elements d^L,jdqidqj does not vanish. If it does vanish, some procedure 
essentially equivalent to the one given here must be used- 
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ed. When the action principle (5J = 0 is written to obtain 
equations of motion, we obtain first 


dl 


J 


dh dt 


. . , ,. oH ^ dH -1 

oPiqi -hPiSqi (9.69) 


Now let the Lagrange multipliers be We add 




U 


m 


^Sq. + ^ dv 
. Hi ^ Sp, 


(9.70) 


to the integrand of (9.69) and then set the result equal to 
zero. For variations which vanish at the end points, the 
equations of motion, 



- - _L. ^ ZJ^ 

SPi ^ “ dp, 



Hi ” Hi 


(9.71) 


then follow. The are functions of p’s and ^’s. The time 
derivative of any dynamical variable, say A, is given^ by 



BA . 
dp, 


lA, H] + u„,[A, 


(9.72) 


Some of the functions represent real arbitrariness in 
the problem. For these, different choices of the arbitrary 
functions correspond to the same physics. Such “constraints” 
are called first class by Dirac and he provides a systematic 
procedure for singling them out. Other types of constraint 
may represent real restrictions, for example on the allowed 
initial values, and these are also taken properly into account 
by Dirac’s procedure. Since 9 ?^ vanishes for all time it follows 
from (9.72) that (letting A = 9 ?,^.) 


L9m'> Pm\ == <> (9.73) 


i Ihe constraint relations must not be employed anywliere in the cal¬ 
culation of Poisson brackets; that is, the jp and q are considered independent 
in this instance. 
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If (9.73) can be solved for the the problem is finished, if 
not the equations (9.73) may give additional relations be¬ 
tween the p/s and g/s such as 

p) ^ 0 (9.74) 

We may again write 

[x^, H] + u„[x^, <p„-] = 0 (9.75) 

These expressions may be employed to solve for u^, obtaining 

“m = U„{q, p) (9.76) 

A study of (9.73) and (9.75) indicates that any 
solution of 

<Pm] ^ VmlZk ^ T’m] ^ ^ (9.77) 

may be added to (9.76). The general solution is therefore 

'^m = (9.78) 

Vam is a set of independent solutions of (9.77). 

Here the Va are arbitrary coefficients, functions of time. 
By use of (9.78) we may eliminate the introducing the 
as variables. These plus the independent p^ and may be 
less than the initial number of p and q variables. (9.78) and 
(9.72) give for the general equations of motion 

A = [A, if'] + t;JA, pj (9.79) 

with 

H’ = H+U„ <p^; <p„ (9.80) 

Formally a constraint is said to be first class if its 
Poisson bracket with the Plamiltonian and all Xk T’m 

vanishes. The other constraints are called second class. 
Other authors have called the equations primary con¬ 
straints and the equations secondary constraints. 

(9.73) indicates that the Poisson bracket of with H' 
vanishes. Also 


H] ^ H] ^0 


(9.81) 
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(9.81) is seen to vanish in consequence of (9.77) and (9.78), 
after multiplying (9.73) by Also (9.73), (9.75), and 

(9.78) show that the Poisson brackets of H' with X/c T’m 
vanish. H' is therefore first class. The Poisson brackets of 
with and X/a vanish in consequence of (9.77). The 
Poisson bracket of with H vanishes because of the vanish¬ 
ing of its Poisson bracket with if', (p^ is therefore first class. 

Equation (9.79) is seen to be made up of the first-class 
functions if' and <Pa- The are arbitrary, not restricted by 
equations of motion- The number of arbitrary functions of 
the time v^, in the general solutions is seen to be the number 
of independent first class (p's. Dirac remarks that in practice 
we know what arbitrary functions there are in the general 
solution, in consequence of the invariance properties of the 
action function. 'Phis allows us to identify the first-class func¬ 
tions without calculating Poisson brackets. As we shall see 
later one may need to do this before the Hamiltonian is known. 

For a field we may imagine all space to be divided into 
cells. (9.65) is still valid, but is being varied only in a 
particular cell of volume zIt. 'The Hamiltonian is 

(9.64A) 

The momentum density is given by 


dLjr . 


(9.65A) 


We may write (i is summed over all variables) 

ff — j (9.64B) 

and introduce the Ilarniltonian density by 

.r* ^ — Lp (9.64C) 

From (9.f)4B) we have (/ is summed from 1 to 3) 

r/ . dLi., BL. \ ,, 

dH = J - dq, - bq, - J d'X 

17 'dL,.- i) iILp . \ 

= J (</.■ -I- a - '5-7.) ) 
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On the far right of (9.64D) we have omitted a surface integral. 
Employing the field equations (4.8) reduces (9.64D) to 

SH = j* (^,-< 571 ,- — :fcfiq,)d^x ( 9 . 64 E) 

From (9.64E) it is apparent that is a function of 
and so we may write 


(9.64F) 


J[( 


(BJe B 
Bq^ dx^ Bq 


r\ BM' 


idM’ dse 

dx^ dn 


-) 


djt. 


d^x 


Again surface integrals have been omitted on the far right 
of (9.64F). Comparing coefficients of dq^, dn^ in (9.64E) and 
(9.64F) then leads to the field equations 

(9.66A) 

B B^\ 

71. 


/a^ 

a dje\ 


” 8a!‘ 8qi,J 


To obtain the time derivative of some quantity given 
as the three volume integral of a density S we write 




dS . dS , dS . dS 

Qi —h q^^ j —)— —|— 


BS' . \ 


' 87t/-^ ' d7t, (9.67A) 

If (o and depend on higher derivatives of q^ and these 
would also have to be included in(9.67A), (9.64F)and (9.66A). 
We now make use of (9.66A) in (9.67A), omitting surface 
integral terms, and obtain 


f 

[{8^ 

a 8^\ 

183^ 

a 8^\ 

J 

LUg. “ 

” 8a!’ 8qiJ 

\ 

Bx^ dTZi^J 


/B<f 

B 

Big \ 

/BJf’ 

B dj^\i 

\d7Zi 

" Bx^ 

8^1, J 

\ 8qi ~ 

" 8x! 8qiJ_ 


[E, H] =J 


d^x 


CVdEdH 

dEdH-^ 


dqi_ 


d^x 


( 9 . 67 B) 
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This expression defines the Poisson brackets for fields The 
functional derivative notation on the extreme right is em¬ 
ployed by some authors. If E is not available as a volume 
integrated density, it may readily be converted into one bv 
Simply using a delta function. 

If there are constraints relating the q, and n,, we may 

proceed as in the deduction of (9.71). The result is that the 
term 

is added to the right side of the expre.ssion for and the term 

\ dq^ doD^ dq^ J 

is added to the right side of the expression for tz,, in (9 66A) 
The equation of motion (9.67B) becomes * 

./ \ dos Bn ■ , / I Btt . Pmi ^ j ~7^ ■ I 


cr/dS' 

^ ""J l(^ 


^ JP \ / o 


djti Ba'J Btz,^ J\Jq~ " a? Yq~ 


B Bje 


J l\dq, B.v^' Bq,^j/\ djt, Bw^ 

{dn, ajr,, J \ dq, dx^ 'dq^_, ^ - Sxi j J d‘ 

(9.67C) 

and the result corresponding to (9.79) is 

h ~ [K, H'] + J qJ_E, q>„^(Ex (9.79A) 

In general relativity the field variables are those needed 
to describe the geometry of the three space at any given time. 
I he field history is then a succession of curved three spaces 
It IS therefore reasonalile to search for a canonical formalisni 
which will end iq) witii only six g,, and their canonical mo¬ 
menta as the re(juired variables. 

Consider now the gravitational Tagrangian density^ 

t In this set'tion we are using units for which mnOje^ = 


1 . 
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— ^ — r^vT^p^r) a /”^ 

= (9,82) 

+ 2(g-“«g^/> — ^/—g 

The Lagrangian is J [V—g] The field variables are the 
gp.vl gfLv,(i are the velocities. The are functions of 

alone, independent of the g^,;,o* The momenta canonical¬ 
ly conjugate to are defined by the statement 

J dsea [ V-g] = J 4V-..od*® (9.83) 

A study of shows that it does not contain any 

terms which are quadratic in g^o,o* follows therefore that 
a function of g^v and ga,v,i, which we denote by //*. 

^ f^gau, go,v,i) (9.84) 

The functions are constraint relations. It can now be 
shown that the equations of motion in no way restrict g^o.o 
and that the constraints (9.84) therefore are “first class. 
If the infinitesimal coordinate transformation 

cc^ 

is carried out, 

gg p,v{^ ) ^gpv ~ g/ip^^,v "i“ gvp^^,p (9.85) 

Let — c)'^pp{x^), where ^p(x^) is an arbitrary func¬ 

tion. Then on the hyperplane = c we have 

bg^^> = 0; ag^^,o == (gfip^^^ + gup^p^)P^ 

4 - pQ 

Since p is arbitrary it follows that an arbitrary change in 
g/ 40,0 does not affect or g,.,.,o* 

It is convenient (following remarks of V. Bargmann) to 
note that some variables of the four-dimensional manifold 
will not change their values on the given spacelike surface 
for coordinate transformations which do not affect the coor¬ 
dinates on the spacelike surface. Such variables are called 
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intrinsic. More precisely, if we carry out the coordinate 
transformations 


~ x^, x^^ x^) 

and Ff^ satisfies the conditions 


F^{c, x'^, x^, x^) ~ c' 
F\c, x^, x^, x^) = x\ 


the variables which are intrinsic do not change their values. 
Thus the are intrinsic. Ihe and are not intrinsic 
and may be varied arbitrarily by choice of F^ or particular 
coordinates. 


We may expect that the intrin.sic field variables will play 
the central role in the Hamiltonian formulation of general 
relativity. 

The 1 oisson bracket relations to be employed here are 
(9.67B). 


[to a/?? — 1 ( 3 ^*' -f- < 3 ^/^ (3a^')('53(a7' — x) 

is the value of g^^ at = c, x'\ x'\ 

Dirac notes that a study ol q indicates that it can be 
modified to make vanish, without affecting the 

equations of motion. Tliis can be done by adding the ex¬ 
pression 

V-g= {(V^ggOOXog^O/gOO}^,, _ {(VZgg00)_,g*0/g00J 

(9.86) 

to iif£,i/_g. g#a7 ciin he written as a surface in- 

tegral, so it will not arfect the equations of motion. The new 
Lagrangian density function ^o^-s/—g is obtained after a 
somewhat tedious reduction as 

+ u -- (^Va, « - h r { gah, v — {gafi, b F gbfi, a 

+ i ^oin, i {- g^^^ ) g^^ H- 2 (g/^« g^^i — )g^^} ] V—g 

(9.87) 



176 


GENERAL RELATIVITY AND GRAVITATIONAL WAVES 


with 

^iiv — gfiv _ ( 9 . 88 ) 

Carrying out the operations listed in (9.83) with oSP* instead 

of gives 

6£(’o*V-g = V-g) (9.89) 

So the momenta Tt'^ are 

i/—g (9.90) 

which may also be written 

r^ah V—g= (gragsb “ ^grsgab)^''^ (9.91) 

Suppose that we have the hyperplane — c and some 
intrinsic field variable ?] such that 97 depends on the and 
on the field variables for other than gravitational fields 
which may be present. Displace the surface = c so that 
each point goes into xt^ + a/^, with af^ small and a function 
of x^ alone. 


dr] == j af^ (Px 

(9.92) 

with a function of x^: 


= / ^ 0 ^'^ 

(9.93) 

Let be a unit normal to aP = c: 


h = = 0; == - 1 

(9.94) 

An expression for P with these requirements is 


0 

0 

il 

(9.95) 

Let be defined by the relation 


— In = = g'‘®(— g®“)“^f/, 

(9.96) 

^0 = ( - 

(9.97) 

Recalling (9.88) enables us to write 


«"g.a = e’''‘g/,« = <5/ - a^o/go® 

(9.98) 

e”g»c = a/ 

(9.99) 

= -g'^/g"" 

(9.100) 
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Using these relations (9.97) becomes 

fo = (- (9.101) 

and (9.98) becomes 

^ = j (9.102) 

We want to have d'rjldcc^ — [r], H] and this is accom¬ 
plished by writing for the Hamiltonian 


m 

H= (9.108) 

provided 

bh^L\ = = 


and will have both matter and gravitational parts 
and depend only on the intrinsic variables. The standard 
definition of the Hamiltonian gives for the gravitational one, 
denoted hj Hq 

Ho = J(’'”gr»,o — ^o*V—S,)'Hx (9.104) 

using (9.87) and (9.90) this is expressed as 

Ha = J^’'’(g„,o + rojsfO)(Px 4- iJ[-g]* [(«”«'*'’ - e'‘‘e‘’>) 

■ {gr.,.ug“"-(g™,.+ g.«,r)g““} {g^^gal,.. - (gaAb +gM,a)g''"}/g®" 

— gi.v. i ga/i. j {(g"- g“'’—g'‘“g‘'^) g" - 52 (g>“ g‘‘^-g/“‘ g^* )g*'^} ] d?x 

(9.105) 


Using (9.91), (8.29) and (9.100) the first integral on the right of 
(9.105) is written 


J w”(g„.o + l\Jg ^) = J{( “- g®“) “M - g) (gra gsi - igrngal>W’^'’ 

-I- g..0«“”K‘‘gr«.. — ‘'i(’*”gr.).8]}«f’» (9.106) 

Dirac gives a procedure for eliminating the velocities and 
reducing the remainder of (9.105) to the canonical form 
(9.108). A study of the S(‘C()nd integrand of (9.105) indicates 
that no term of the form with only a function 
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grs. ^rs,y Originate there. It is therefore necessary to 
calculate only the contribution to (— Since 
also depends only on the intrinsic variables it may 

be calculated assuming g^Q == 0. The determinant g may be 
written in terms of the determinant of as gg^^ = ^g. 
After omitting a surface integral the result is 


+ gmoe^'’i^^^grs,v - 2(7r-g,J,, ]] (9.107) 


with =3^jyrjj;, and the parts of the matter Hamiltonian cor¬ 

responding to the parts and of (9.103). There are 
thus only six degrees of freedom corresponding to the g,.,. left, 
since tt/*® vanish in consequence of the modification of 
which arranged for (9.84) to vanish. Also the relations (9.103) 
plus the requirement H] = 0 leads to and van¬ 
ishing in consequence of the constraints. 

In the weak-field approximation the Hamiltonian which 
includes terms up to (A^^)^ is given by 


H 




i —rs—TS 

171 7C 




~t~ 4 grs, u g', 


rs,u 


4:grr, u gs3,u 


1 2 grs, r guu, s ^grs, u gru, s -f-e^ ml] d^X 

J ([ ^oo] ^ 1 ) (^rs.rs— grr.ss — Ml) d^X 

-/groCajl”* -jeMr)d?X 


(9.107a) 


For a particle of rest mass m the action is 


I 


f 


m ds 


m 


r/_<, 

Jl dT dT) 


where y/^ are the coordinates of the particle. 


L 


V 


m 


( 


dyf^ dy^\^ 
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difldT = 1 and is not a dynamical variable. The momentum 
Pj. is given by 



m<> dyf\~h 

dT V df dT j 


The particle contributes a term 


Hm — PrV^ — Lm 

to the Hamiltonian. 

This may be written in the form 

Hm = —gsO^'^^P^ + (—_j_ 

From (9.103) it follows that 

<3^ml = -f- Pr ^s)^ — y) 

y^Mr ~ — — y) 

A different approach to the Hamiltonian formulation of 
general relativity has been given by Arnowitt, Deser, and 
Misner. Their results are as follows.^ 

A canonical form for general relativity involving only 
two pairs of unconstrained conjugate variables has been 
obtained. The action 


•^ = / [-g]* Rd^x 

expressed in first-order form is 

i-g)iR=- 0 + (“g" R + in^- 

+2g„,»g)‘([-gOO] -1)., go,n (9-108) 
Here 


= (-g)^ - gi„. 

is the matrix inverse to is constructed only of 


t With minor ultorutions this section is essentially a reprint of the 
letter to the editor of Nuoxm cimmto by Arnowitt, Deser, and Misner ( 28 ). 
Reprinted by permission of the authors and Nuovo ciniento. 
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n is the trace of and the bar notation is used for the 
covariant derivative with respect to the metric. Variation 
of this action gives rise to all the usual Einstein field equa¬ 
tions in Palatini form. End-point variations give the generator 

G = j [— dcc^] (Px 

where the stress tensor terms, J dx^, vanish as a conse¬ 
quence of the constraint equations = 0 

(obtained from varying Eq. (9.108) with respect to 
(— g00)-^goz)* 

Now make use of an orthogonal decomposition of and 
according to the general scheme 


/ • • = / • -f- 

J t3 / t3 I 




+ fi,3- + f3,z 


(9.109) 


where l/V^ is the inverse of the flat-space Laplacian operator 
with appropriate boundary conditions. 





(9.109a) 

(9.109b) 


and the transverse traceless part of i.e., is the re¬ 

mainder. 


Such a decomposition is meaningful only when a set of 
coordinate conditions are imposed. The following are 
selected: 

giij = = 0 (9.110) 

General covariance is maintained, since imposition of these 
coordinate conditions is equivalent to using certain invariant 
functionals of the metric as independent variables in place of 
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coordinates.^ Conditions (9.110) imply that the functionals 


t 




(9.111) 


are chosen. These ensure an asymptotically flat metric. 

Inserting the orthogonal decomposition into the genera¬ 
tor we find after some partial integrations that 







]d^x 


(9.112) 


since both the variation d and the time derivative in the ac¬ 
tion commute with the decomposition. Here and j 
may be expressed in terms of and by'solving the 
four constraint equations ^ G/ = 0 for them. We see, there¬ 
fore, that the generator (9.112) is in standard canonical form 
G = pdq — Hdt with the additional momentum terms r,® 
characteristic of a field theory: 





H- TJ>da}^]d^cc 


(9.112a) 


Since and t do not appear explicitly in the expressions for 
gi^ and they do not appear explicitly in the equations 
Gf^^ — 0. Detailed calculation shows that they also do not 


t Thus, as is well known, the apparently noncovariant DeDonder con¬ 
ditions, ([ V -g] gf^^),v = 0 are equivalent to the covariant statement that four 
linearly independent scalar functions are to be used as tlie coordinates a>«, 
where = 0. Sueli an equivalence between coordinate conditions and 

covariant statements holds for this case and indeed for any other choice of 
coordinates which does not depend on the of some initial frame. 

tt Aj^propriate boundary conditions have, of course, been imposed on 
the operator (1/V^) in order to relate (D.HO) and (9.111). A simple proce¬ 
dure making this operator manifestly nonsingular consists in rewriting (9.111) 
as -(l/2V'q7r^, x’ier<xr-^g^ Hubsequently taking the limit a“>() 

at the end of the analysis. 

t That Gfjp ■-= 0 has at least perturbation series solutions for ami 

can easily be checked, as (^an the fact that these coordinate iionditiojis are 
maintained in such an expansion. The Hamiltonian can thus be exhibited as 
an infinite series in powers of the canonical variables beginning with the 
linearized theory’s Jf. 
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appear in the solutions for and ji*' Thus the quanti- 
ties 


= -g^. u{gir> 

iTf = „rsTT^ (9.H3) 

do not depend explicitly on the coordinates. This coordinate 
independence allows us to derive the standard conservation 
laws. From (9.112a) and the fact that the action of Eq. 
(9.108) reduces to 


/ — J ~ dJ^x (9.108a) 

it follows that the two independent pairs of canonical vari¬ 
ables are and These, therefore, obey simple 

Poisson bracket relations 

femn^^(r), (9.114) 

where the transverse traceless d function, is defined 

as in linearized theory and is, of course, independent of the 
metric (31). 

The energy E of the field is defined to be the numerical 
value of the Hamiltonian for the given solution of the equa¬ 
tions of motion. Thus, in the evaluation of 


^ (9.11S) 

be expressed in terms of the canonical vari¬ 
ables and so E can be evaluated as a surface integral (even 
t lough the Hamiltonian cannot be reduced to this form). 
From Eq. (9.109a) we obtain 


E == 




Si3,3) 


(9.116) 


here dS^ = ^e^^dx^dx^ is the two-dimensional surface 
element at infinity in rectangular coordinates. We have 
assumed that the metric becomes asymptotically flat and 
tRe coordinates rectangular at spatial infinity. Coordinate 
transformations preserving these boundary conditions can 
ere e rigorously treated by the linearized theory where 
is a coordinate scalar (31). Thus, it is not necessary to express 
the metric in the canonical coordinate system to evaluate 
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the energy; neither is it necessary to use rectangular coordi¬ 
nates provided we make the usual flat space tensor trans¬ 
formations to the desired (e.g., spherical) coordinates. This 
definition effectively states that the energy of a closed sys¬ 
tem may be obtained from the coefficient of 1 /r in the asymp¬ 
totic expansion of It is also the gravitational mass of 
the system as seen by a distant test particle (due to the 
boundary conditions). The constancy of the energy ensures 
that this coefficient is time-independent. Equation (9.116) 
holds without modification when a point particle is coupled 
to the gravitational field; E now represents the total energy 
of the coupled system. For the simple case of the Schwarz- 
schild metric one of course finds that the energy is the mass 
parameter. 

It is a desirable criterion for the energy and momentum 
of a closed system that they involve only those quantities 
required to specify the state of the system at a given time. 
For general relativity these variables are and when no 
coordinate conditions are employed but do not include go^, 
which serve to describe how the coordinates will be chosen at 
a later time. The expression (9.115) for the energy satisfies 
this condition. Analogous considerations to the above hold 
for the total momentum of the gravitational system where 
again is a coordinate scalar in the asymptotic region. An 
example of an initial value problem where go^ is not deter¬ 
mined is the “many body” spatial metric of Lichnerowicz 
(33): 

gii = + 2 --- aj]^ (9.117a) 

7^^ = 0 (9.117b) 

Here Eq. (9.116) yields unambiguously the value 2 for 
the energy. 

This formalism also allows us to establish criteria for 
the existence of gravitational radiation; these are in exact 
correspondence to tliose of electromagnetic theory. Thus, 
the condition for radiation escaping to infinity is stated in 
terms of the Poynting vector there. The energy flux across 
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a surface element dS^ at infinity is 

= — 2n^\^dS^ (9.118) 

where, again, the right-hand side may be evaluated in any 
asymptotically rectangular coordinate system. More gener¬ 
ally, the criterion for the existence of gravitational waves at 
any point is the presence (in the canonical coordinate frame) 
of a nonvanishing or there. Alternately, such a 

situation means that there exists an excitation of the gravi¬ 
tational field in one of its canonical variables (which are in¬ 
dependent of the source variables). This is identical with the 
electromagnetic definition of a wave, which requires the 
existence of or the canonical variables of the 

electromagnetic field. In the obvious case where no waves 
are expected to exist, namely, the Schwarzschild solution, 
one can verify that g^^^ and vanish everywhere in the 
canonical frame, which justifies its being regarded as a 
“one-particle” solution. This has not been established for 
the “many body” initial conditions treated above. 

9.6 Remarks on Quantization of General Relativity 

It has been noted by Landau, by Klein, and by Pauli 
that quantized general relativity would have to be considered 
in elementary particle theory, at high energies. A cutoff 
momentum P for quantum electrodynamics is suggested, 
with 

In 1/a (9,119) 

Here a is the fine-structure constant and m is the electron 
rest mass. The result (9.119) was given by Kalian and by 
Pauli as a consequence of the requirement of mathematical 
consistency. Landau noted that for energies corresponding 
to the cutoff momentum, the gravitational interaction is 
already of the same order of magnitude as the Coulomb 
interaction. In other words, as the energy increases the mass 
and gravitational interaction increase. The charge, however, 
does not change. 

The quantum principles are very well established. They 
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take into account the ‘‘granular” structure of nature and the 
interaction between measuring devices and the system being 
measured. As P. G. Bergmann has remarked, unquantized 
general relativity cannot peacefully coexist with quantum 
theory. The gravitational field is the metric tensor; in this 
sense it is the most fundamental field in physics. Its quanti¬ 
zation may affect, in a profound way, our notions of space 
and time, particularly in the realm of small dimensions (high 
energy). 

It is interesting that it is possible to construct a quantity 
having the dimensions of length, out of G, and c, it is 

L = iO“33 (9.120) 

This is much too small to be directly related to elemen¬ 
tary particle dimensions now known. If ^ we consider quanti¬ 
zation using Feynman’s formulation we should write for the 
wave function for' the field 


V^later = \ c') ?/yearlier(c')dc' 

where the kernel is given by 

(9.121) 

H 


with H meaning a sum over all field histories. For the gravi¬ 
tational field. 


ilfjlfi 


exp 


% 

Ih 


r J* ^ 


(9.122) 


H 


The Feynman sum should in principle include all possible 
metrics and topologies as acceptable histories, consistent 
with the Lorentz signature. Since 


R 



2 


it follows from (9.122) that if we consider a region of four¬ 
dimensional space-time, L units on a side, the quantum 


t Tliis discussion follows tlic conunents of J. A. Wheeler. 
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fluctuations would be expected to be of the order 

Also we might expect perhaps amelioration of the difficulties 
of quantum field theory (18). There are a number of func¬ 
tions employed there which are singular on the light cone. 
Consider two points which lie on a light cone in flat space. 
If the space is now curved the second point will no longer 
be on the light cone. If a Feynman sum is carried out, over 
all metrics, the difficulties associated with the light cone may 
turn out to have measure zero. This indicates, unfortunately, 
that any quantization scheme which starts with a flat-space 
linear approximation will immediately have all the diver¬ 
gences. To remove them, it may be necessary to quantize 
the nonlinear theory from the very beginning. 

9.7 Spinors in General Relativity 

We discuss briefly the formalism used by Bargmann* to 
discuss spinor fields in general relativity. A field of y ma¬ 
trices is introduced which satisfy the anticommutation law 

yfcYV yvy[I ( 9 . 124 !) 

I is the unit matrix. The components of y are assumed to be 
continuous functions of the coordinates, and to transform 
like a vector under coordinate transformations. Under a 
spinor transformation 

(Spinor )new = (Spinor )oid (9.125) 

the y transform like 

(ya)new ~ *5^ ^(ya)old‘S' (9.126) 

It is more simple to deal only with real representations of the 
y together with spinor transformations 5 whose matrix ele¬ 
ments are real. 

The following formalism is arranged to guarantee that 

t As discussed by D. Brill and J. A. Wheeler, Revs. Modern Phys. 29, 
465 (1957). See also references 42, 43, 44. 
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tensors will transform like tensors, that the effect of a spinor 
transformation will be as given by (9.125), and that the effect 
of a spinor transformation on a spinor tensor will be 

/iip)new “ ‘S' (9.127) 

For covariant differentiation additional 4x4 matrices 
Fa are introduced. These F’s are determined from (9.124) up 
to a multiple of the unit matrix by the relation 

rj,,v — ( 9 . 128 ) 

The CO variant derivative of an object with spinor transfor¬ 
mation properties with respect to is written and has 
the properties 


V^(/4*) = (V^^)* (9.129) 

== 0 

This * means transpose of the complex conjugate. Linear 
combinations of the y matrices used in special relativity can 
be employed. These are written y^ and satisfy 

+ Yvy/. = y.* = y,; y^* = (9.139) 

A system of Vierbein is often used. A locally Lorentz metric 

is introduced at each point by writing 

dx^ = h^t^dxt^ (9.131) 

Here the are the Lorentz metric coordinates and are the 
general coordinates. 

A set of y’s satisfying (9.124) is 


yf. = h/y^ (9.132) 

An expression for the which is a solution of (9.128) is 


r 

-*■ /ti 4 


(f 

t>KX 


L\ dxf^) 






'-^a. 


]IV 


a (9.133) 


Here is arbitrary and = Ky/^y*^ ~~ y^y/*). If the Vier¬ 
bein formalism is used, a similarity transformation of the 
spinors corresponds to a Lorentz transformation of the Vier¬ 
bein. 
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The covariant derivative of a spinor y) is given by 

- r,y> (9.134) 

The covariant derivative of a spinor tensor is 

r, - (9.135) 

The Panli conjugate of tp is written and defined by 

(9.136) 

with p a matrix which is Hermitian and which is selected so 
that pi is also Hermitian. In the Vierbein formalism 
p = iy^. Real quantities such as current density can be 
written 

/^ = yj^iyf^ip (9.137) 

The Dirac equation in general relativity is 

= 0 (9.138) 

with the arbitrary traces of the arranged to include the 
effect of the four potential. (9.138) may be obtained from 
the variational principle. 

dj [yy^y^^i'^ccW) + p('y>^y}]V—§ (Fee — 0 (9.139) 

with and y) varied independently. 

The assumption that the atom is a natural clock may be 
justified from these relations. We may introduce a Lorentz 
frame along the world line of the atom. The flat space (9.138) 
is valid since the atomic dimensions are too small for the 
Riemann tensor to affect the internal motion. To compare 
the radiation with that of a distant atom, a coordinate trans¬ 
formation is carried out so that only the time coordinate is 
changed, and the same time coordinates are employed for 
atoms in different places. From (9.124), y® — y°(~g®^)^. 
The energy is associated with y^^djdx^^ so the energy levels 

become = F^®(-goo)^ in agreement with the red shift 
formula. 
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Exercises 


Just do the exercises diligently^ then you will find out what 
you have understood and what you have not. 

A. Sommerfeld (to W. Heisenberg) 


1 . Discuss the co variant and contra variant components of 
a vector in a plane in which oblique coordinates are 
employed. 

2 . Deduce expressions (3.30) and (3.31), the transforma¬ 
tion laws for the Christoffel symbols of the first and 
second kinds. 

3. Show that the number of independent components of 

ill S' space of N dimensions is {N^ll2){N^ — 1 ). 

4. Discuss the embedding of a curved space of N dimen¬ 
sions in a flat space of a greater number of dimensions. 

5. Derive expression (3.49) for Ra^yd- 

6 - Set up a coordinate system on a torus; calculate all 
components of Ra/Syd' 

7. Evaluate s^fiys 

8 . A contra variant vector is being parallel displaced along 
a small circle of a sphere. It has components = 0 
and A’’ — 2 at the point 6 — 7tl4i, 9 ? = 0 . The path 
along which it is being displaced is 0 = constant. Cal¬ 
culate the values of and as a function of 99 along 
the small circle. 

9. An observer in a uniformly accelerated frame measures 
time using the synchronized clocks fixed in an inertial 
frame with respect to which he is being accelerated. 
Calculate for the accelerated observer. 

10 . Carry out the operations leading to (4.15). 

11 . Derive the expression (5.6). 
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EXERCISES 


12. Derive (5.13a) by starting- with a four potential with 

== — 0 and Aq = f{r). 

13. Write the exact red-shift formula given in expression 
(5.21) as a series expansion, compare it with a series 
obtained using the relativistic Doppler effect in the 
equivalence-principle deduction of Chapter 1. To what 
order in the change in gravitational potential do they 
agree? Does this mean that general relativity is in dis¬ 
agreement with the equivalence principle? 

14. Derive expressions (9.1), (9.2), and (9.3). 

15. Show that (9.44) gives a value of K which is the same for 
any linear combination of A and B, 

16. Given that the Lagrangian is 

^ ~~ 9 ^ 0 ^ + + ^ 3^3 + ^4 

with and functions of q^, and q^. Solve for 

the equations of motion. 

17. Carry out the reductions leading to (9.87), (9.91), (9.95), 
(9.106), (9.107). 

18. A satellite is equipped with a natural clock which is to 
be compared with an identical clock at rest on the 
earth’s surface. If we compare elapsed intervals at the 
end of a given coordinate time, show that (Singer, Phys. 
Rev. 104, 11, (1956)). 


1 

1 

H 

GM 

a —2hlr 

30)2 \ 

ds-^ 

2c^r 

Vl+ hjr 

^TtGpMJ 


is the interval for the earth clock, Sg is the interval for 
the satellite clock, h is the height of the assumed circular 
orbit above the surface of the earth, M is the mass of the 
earth, pm is the mean density of the earth, r and o> refer 
to the radius and angular velocity, respectively, of the 
earth. 

19. Identify the terms in (9.107) which are the three dimen¬ 
sional scalar curvature. 
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